Replace drift heating to engbalg outside ngsp loop

holm10 · Nov 9, 2024 · 20a1bde · 20a1bde
1 parent a1354c4
commit 20a1bde
Showing 1 changed file with 15 additions and 15 deletions.
diff --git a/bbb/oderhs.m b/bbb/oderhs.m
@@ -4668,21 +4668,6 @@ ccc             call volave(nx, ny, j2, j5, i2, i5, ixp1(0,0), ixm1(0,0),
      .                  + sead(ix,iy)
      .                  + seadh(ix,iy)
 
-            upgcc = 0.5*(up(ix,iy,iigsp)+up(ix1,iy,iigsp))
-            vycc = (cfnidhgy**0.5)*0.5*(vy(ix,iy,iigsp)+vy(ix1,iy,iigsp))
-            v2cc = (cfnidhg2**0.5)*0.5*(v2(ix,iy,iigsp)+v2(ix1,iy,iigsp))
-
-c           Atom kinetic energy source from mol. drift heating
-            reseg(ix,iy,1) = reseg(ix,iy,1) 
-     .          + cfnidh*cfnidhdis*0.5*mg(1)* (upgcc**2 + vycc**2 + v2cc**2)
-     .          * psordis(ix,iy,2)
-
-c           Ion energy source from mol. drift heating
-            resei(ix,iy) = resei(ix,iy)
-     .          + cftiexclg * cfneut * cfneutsor_ei * cnsor * cfnidhdis
-     .          * 0.5*mg(1)*(upgcc**2 + vycc**2 + v2cc**2) * psordis(ix,iy,2) 
-
-
             else
                resei(ix,iy) = resei(ix,iy) + w0(ix,iy)
      .             + cfneut * cfneutsor_ei * ctsor*1.25e-1*mi(1)*
@@ -7846,6 +7831,21 @@ call fd2tra (nx,ny,floxge(0:nx+1,0:ny+1,igsp),
 c        Special case for igsp = 1
 c        Check inertial neutral model
 
+            upgcc = 0.5*(up(ix,iy,iigsp)+up(ix1,iy,iigsp))
+            vycc = (cfnidhgy**0.5)*0.5*(vy(ix,iy,iigsp)+vy(ix1,iy,iigsp))
+            v2cc = (cfnidhg2**0.5)*0.5*(v2(ix,iy,iigsp)+v2(ix1,iy,iigsp))
+
+c           Atom kinetic energy source from mol. drift heating
+            reseg(ix,iy,1) = reseg(ix,iy,1) 
+     .          + cfnidh*cfnidhdis*0.5*mg(1)* (upgcc**2 + vycc**2 + v2cc**2)
+     .          * psordis(ix,iy,2)
+
+c           Ion energy source from mol. drift heating
+            resei(ix,iy) = resei(ix,iy)
+     .          + cftiexclg * cfneut * cfneutsor_ei * cnsor * cfnidhdis
+     .          * 0.5*mg(1)*(upgcc**2 + vycc**2 + v2cc**2) * psordis(ix,iy,2) 
+
+
 
 c        Special case for molecules