Changeset 2475 for LMDZ5/trunk/libf/dyn3dmem/initdynav_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/initdynav_loc.F (modified) (10 diffs)

LMDZ5/trunk/libf/dyn3dmem/initdynav_loc.F

@@ -154 +154 @@
 
       ddid=(/ 1,2 /)
-      dsg=(/ iip1,jjp1 /)
+      dsg=(/ iip1,jjm /)
       dsl=(/ iip1,jjn /)
       dpf=(/ 1,jjb /)
@@ -171 +171 @@
      
 ! Grille U
+
+      do jj = 1, jjp1
+        do ii = 1, iip1
+          rlong(ii,jj) = rlonu(ii) * 180. / pi
+          rlat(ii,jj) = rlatu(jj) * 180. / pi
+        enddo
+      enddo
 
       jjb=jj_begin
@@ -209 +216 @@
 C  Vents U
 C
+      jjn=jj_nb
       call histdef(histuaveid, 'u', 'vent u moyen ',
-     .             'm/s', iip1, jjp1, uhoriid, llm, 1, llm, zvertiidu,
+     .             'm/s', iip1, jjn, uhoriid, llm, 1, llm, zvertiidu,
      .             32, 'ave(X)', t_ops, t_wrt)
 
@@ -216 +224 @@
 C  Vents V
 C
+      if (pole_sud) jjn=jj_nb-1
       call histdef(histvaveid, 'v', 'vent v moyen',
-     .             'm/s', iip1, jjm, vhoriid, llm, 1, llm, zvertiidv,
+     .             'm/s', iip1, jjn, vhoriid, llm, 1, llm, zvertiidv,
      .             32, 'ave(X)', t_ops, t_wrt)
 
@@ -223 +232 @@
 C  Temperature
 C
+      jjn=jj_nb
       call histdef(histaveid, 'temp', 'temperature moyenne', 'K',
-     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
+     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
      .             32, 'ave(X)', t_ops, t_wrt)
 C
@@ -230 +240 @@
 C
       call histdef(histaveid, 'theta', 'temperature potentielle', 'K',
-     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
+     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
      .             32, 'ave(X)', t_ops, t_wrt)
 
@@ -238 +248 @@
 C
       call histdef(histaveid, 'phi', 'geopotentiel moyen', '-',
-     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
+     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
      .             32, 'ave(X)', t_ops, t_wrt)
 C
@@ -251 +261 @@
 C  Masse
 C
-      call histdef(histaveid, 'masse', 'masse', 'kg',
-     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
+      call histdef(histaveid, 'masse', 'masse moyenne', 'kg',
+     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
      .             32, 'ave(X)', t_ops, t_wrt)
 C
@@ -258 +268 @@
 C
       call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa',
-     .             iip1, jjp1, thoriid, 1, 1, 1, -99,
-     .             32, 'ave(X)', t_ops, t_wrt)
-C
-C  Pression au sol
+     .             iip1, jjn, thoriid, 1, 1, 1, -99,
+     .             32, 'ave(X)', t_ops, t_wrt)
+C
+C  Geopotentiel au sol
 C
 !      call histdef(histaveid, 'phis', 'geopotentiel au sol', '-',
@@ -273 +283 @@
       call histend(histvaveid)
 #else
-      write(lunout,*)'initdynav_p: Needs IOIPSL to function'
+      write(lunout,*)'initdynav_loc: Needs IOIPSL to function'
 #endif
 ! #endif of #ifdef CPP_IOIPSL
-      return
       end