Changeset 175

Timestamp:

10/16/13 11:58:42 (11 years ago)

Author:

ymipsl

Message:

Bug fix in dissipation.
=> better conservation of the mass

YM

File:

• codes/icosagcm/trunk/src/dissip_gcm.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/dissip_gcm.f90

@@ -691 +691 @@
       CALL swap_geometry(ind)
       dtheta_rhodz=f_dtheta_rhodz(ind) 
+      mass=f_mass(ind)
     
       DO  l = ll_begin, ll_end
@@ -781 +782 @@
     DO l=llb,lle
 !$SIMD
-      DO ij=jj_begin,ij_end
+      DO ij=ij_begin,ij_end
 
           divgrad_i(ij,l)=1./Ai(ij)*(ne(ij,right)*grad_e(ij+u_right,l)*le(ij+u_right)  +  &