Changeset 138 for codes/icosagcm/trunk/src/etat0_dcmip4.f90

Timestamp:

02/16/13 17:03:57 (11 years ago)

Author:

dubos

Message:

Transport now working again - tested with dcmip4.1.0

File:

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

codes/icosagcm/trunk/src/etat0_dcmip4.f90

@@ -174 +174 @@
     
       q(:,:,1)=theta(:,:)
+      IF(nqtot>2) q(:,:,3)=1.
 
       DO l=1,llm
@@ -200 +201 @@
                                                             -cos(latc)*sin(lat)*cos(lon-lonc))/(sqrt(1-K**2)))
              q(ij,l,2)=ABS(g/preff*(-1./radius*duodeta*dthetaodlat-(2*sin(lat)*omega+vort)*dthetaodeta))
+             IF(nqtot>3) q(ij,l,4)=cos(lon)*cos(lat)
            ENDDO
          ENDDO