Changeset 30

Timestamp:

07/26/12 18:00:18 (12 years ago)

Author:

sdubey

Message:

bug fix for advection, with transfer

File:

• codes/icosagcm/trunk/src/advect_tracer.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/advect_tracer.f90

@@ -61 +61 @@
     INTEGER,SAVE :: iadvtr=0
     LOGICAL,SAVE:: first=.TRUE. 
-    !------------------------------------------------------sarvesh
+!------------------------------------------------------sarvesh
+       CALL transfert_request(f_ps,req_i1)
+       CALL transfert_request(f_u,req_e1)
+       CALL transfert_request(f_u,req_e1)
+       CALL transfert_request(f_q,req_i1)
+        
     IF ( first ) THEN 
        CALL allocate_field(f_rhodz,field_t,type_real,llm) 
@@ -140 +145 @@
           uc = f_uc(ind) 
           uc = uc/real(iapp_tracvl) 
-          massflxc = f_massflx(ind)
+          massflxc = f_massflxc(ind)
           massflxc = massflxc*dt
           ! now massflx is time-integrated
@@ -207 +212 @@
        CALL advtrac(massflx,wg) ! compute vertical mass fluxes
     END DO
-
-    !   CALL transfert_request(f_wg,req_i1) 
 
     DO ind=1,ndomain