Changeset 1216 for trunk/LMDZ.GENERIC/libf/phystd/lect_start_archive.F

Timestamp:

Apr 3, 2014, 9:09:47 AM (11 years ago)

Author:

emillour

Message:

Generic model:
Major cleanup, in order to ease the use of LMDZ.GENERIC with (parallel) dynamics
in LMDZ.COMMON: (NB: this will break LMDZ.UNIVERSAL, which should be thrashed
in the near future)

• Updated makegcm_* scripts (and makdim) and added the "-full" (to enforce full recomputation of the model) option
• In dyn3d: converted control.h to module control_mod.F90 and converted iniadvtrac.F to module infotrac.F90
• Added module mod_const_mpi.F90 in dyn3d (not used in serial mode)
• Rearanged input/outputs routines everywhere to handle serial/MPI cases. physdem.F => phyredem.F90 , phyetat0.F => phyetat0.F90 ; all read/write routines for startfi files are gathered in module iostart.F90
• added parallelism related routines init_phys_lmdz.F90, comgeomphy.F90, dimphy.F90, iniphysiq.F90, mod_grid_phy_lmdz.F90, mod_phys_lmdz_mpi_data.F90, mod_phys_lmdz_mpi_transfert.F90, mod_phys_lmdz_omp_data.F90, mod_phys_lmdz_omp_transfert.F90, mod_phys_lmdz_para.F90, mod_phys_lmdz_transfert_para.F90 in phymars and mod_const_mpi.F90 in dyn3d (for compliance with parallelism)
• added created generic routines 'planetwide_maxval' and 'planetwide_minval', in module "planetwide_mod", that enable obtaining the max and min of a field over the whole planet. This should be further imroved with computation of means (possibly area weighed), etc.

EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/lect_start_archive.F (modified) (16 diffs)

trunk/LMDZ.GENERIC/libf/phystd/lect_start_archive.F

@@ -3 +3 @@
      &     q,qsurf,surfith,nid)
 
-      USE surfdat_h
-      USE comsoil_h
-      USE tracer_h
+!      USE surfdat_h
+      USE comsoil_h, ONLY: nsoilmx, layer, mlayer, volcapa, inertiedat
+      USE tracer_h, ONLY: igcm_co2_ice
+      USE infotrac, ONLY: tname, nqtot
+!      USE control_mod
 
 c=======================================================================
@@ -26 +28 @@
 #include "dimensions.h"
 #include "dimphys.h"
-#include "planete.h"
+!#include "planete.h"
 #include "paramet.h"
 #include "comconst.h"
 #include "comvert.h"
 #include "comgeom2.h"
-#include "control.h"
-#include "logic.h"
+!#include "control.h"
+!#include "logic.h"
 #include "ener.h"
 #include "temps.h"
 #include "netcdf.inc"
-#include"advtrac.h"
+!#include"advtrac.h"
 c=======================================================================
 c   Declarations
@@ -79 +81 @@
       REAL vcov(iip1,jjm,llm),ucov(iip1,jjp1,llm) ! vents covariants
       REAL h(iip1,jjp1,llm),ps(iip1,jjp1)
-      REAL q(iip1,jjp1,llm,nqmx),qtot(iip1,jjp1,llm)
+      REAL q(iip1,jjp1,llm,nqtot),qtot(iip1,jjp1,llm)
 
 c autre variables dynamique nouvelle grille
@@ -97 +99 @@
       REAL co2ice(ngridmx) ! CO2 ice layer
       REAL emis(ngridmx)
-      REAL q2(ngridmx,nlayermx+1),qsurf(ngridmx,nqmx)
+      REAL q2(ngridmx,nlayermx+1),qsurf(ngridmx,nqtot)
 c     REAL phisfi(ngridmx)
 
@@ -119 +121 @@
       real co2iceS(iip1,jjp1)
       real emisS(iip1,jjp1)
-      REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqmx)
+      REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
 
       real ptotal, co2icetotal
@@ -283 +285 @@
       allocate(psold(imold+1,jmold+1))
       allocate(phisold(imold+1,jmold+1))
-      allocate(qold(imold+1,jmold+1,lmold,nqmx))
+      allocate(qold(imold+1,jmold+1,lmold,nqtot))
       allocate(co2iceold(imold+1,jmold+1))
       allocate(tsurfold(imold+1,jmold+1))
@@ -296 +298 @@
       allocate(surfithold(imold+1,jmold+1))
       allocate(mlayerold(nsoilold))
-      allocate(qsurfold(imold+1,jmold+1,nqmx))
+      allocate(qsurfold(imold+1,jmold+1,nqtot))
 
       allocate(var (imold+1,jmold+1,llm))
@@ -703 +705 @@
 
 ! Surface tracers:      
-      do iq=1,nqmx
+      do iq=1,nqtot
         ! initialize all surface tracers to zero
         call initial0((jmold+1)*(imold+1), qsurfold(1,1,iq))
@@ -709 +711 @@
 
 
-!      print*,'tnom=',tnom
+!      print*,'tname=',tname
 !      print*,'nid',nid
 !      print*,'nvarid',nvarid
 !      stop
 
-      DO iq=1,nqmx
-          txt=trim(tnom(iq))//"_surf"
+      DO iq=1,nqtot
+          txt=trim(tname(iq))//"_surf"
           if (txt.eq."h2o_vap") then
             ! There is no surface tracer for h2o_vap;
@@ -748 +750 @@
         ENDIF
 
-      ENDDO ! of DO iq=1,nqmx
+      ENDDO ! of DO iq=1,nqtot
 
 
@@ -894 +896 @@
 
 ! Tracers:      
-      do iq=1,nqmx
+      do iq=1,nqtot
          call initial0((jmold+1)*(imold+1)*lmold,qold(1,1,1,iq) )
       enddo
 
-      DO iq=1,nqmx
-        txt=tnom(iq)
+      DO iq=1,nqtot
+        txt=tname(iq)
         write(*,*)"lect_start_archive: loading tracer ",trim(txt)
         ierr = NF_INQ_VARID (nid,txt,nvarid)
@@ -925 +927 @@
         ENDIF
 
-      ENDDO ! of DO iq=1,nqmx
+      ENDDO ! of DO iq=1,nqtot
 
 
@@ -1254 +1256 @@
 c     write(49,*) 'ucov',vcov
 
-      if (nqmx .gt. 0) then
+      if (nqtot .gt. 0) then
 c traceurs surface
-      do iq = 1, nqmx
+      do iq = 1, nqtot
             call interp_horiz(qsurfold(1,1,iq) ,qsurfs(1,1,iq),
      &                  imold,jmold,iim,jjm,1,
@@ -1262 +1264 @@
       enddo
 
-      call gr_dyn_fi (nqmx,iim+1,jjm+1,ngridmx,qsurfs,qsurf)
+      call gr_dyn_fi (nqtot,iim+1,jjm+1,ngridmx,qsurfs,qsurf)
 
 c traceurs 3D
-      do  iq = 1, nqmx
+      do  iq = 1, nqtot
             call interp_vert(qold(1,1,1,iq),var,lmold,llm,
      &        apsold,bpsold,aps,bps,psold,(imold+1)*(jmold+1))
@@ -1305 +1307 @@
 
 c     Periodicite :
-      do  iq = 1, nqmx
+      do  iq = 1, nqtot
          do l=1, llm
             do j = 1, jjp1
@@ -1316 +1318 @@
 ! no need to transfer "co2ice" any more; it is in qsurf(igcm_co2_ice)
 
-      endif !! if nqmx .ne. 0
+      endif !! if nqtot .ne. 0
 
 c-----------------------------------------------------------------------