Changeset 2475 for LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F

Timestamp:

Mar 22, 2016, 1:19:13 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Reinstate writehist, writedyn and bilan_dyn in dyn3dmem so that ouputs in the dynamics (dyn_hist* and dynzon files) may be generated when in parallel.
EM

File:

• LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F (modified) (4 diffs)

LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F

@@ -1603 +1603 @@
                ENDIF
 
+              ! Ehouarn: re-compute geopotential for outputs
+c$OMP BARRIER
+c$OMP MASTER
+              CALL geopot_loc(ip1jmp1,teta,pk,pks,phis,phi)
+c$OMP END MASTER
+c$OMP BARRIER
+
 #ifdef CPP_IOIPSL
              IF (ok_dynzon) THEN 
@@ -1638 +1645 @@
 #ifdef CPP_IOIPSL
              if (ok_dyn_ins) then
-                 CALL writehist_loc(itau,vcov,ucov,teta,phi,q,
+                 CALL writehist_loc(itau,vcov,ucov,teta,pk,phi,q,
      &                              masse,ps,phis)
              endif
@@ -1750 +1757 @@
 
 #ifdef CPP_IOIPSL
+              ! Ehouarn: re-compute geopotential for outputs
+c$OMP BARRIER
+c$OMP MASTER
+              CALL geopot_loc(ip1jmp1,teta,pk,pks,phis,phi)
+c$OMP END MASTER
+c$OMP BARRIER
+               
                IF (ok_dynzon) THEN
                CALL bilan_dyn_loc(2,dtvr*iperiod,dtvr*day_step*periodav,
@@ -1774 +1788 @@
 #ifdef CPP_IOIPSL
               if (ok_dyn_ins) then
-                 CALL writehist_loc(itau,vcov,ucov,teta,phi,q,
+                 CALL writehist_loc(itau,vcov,ucov,teta,pk,phi,q,
      &                              masse,ps,phis)
               endif ! of if (ok_dyn_ins)