Changeset 3978 for LMDZ6/trunk

Timestamp:

Aug 30, 2021, 8:37:06 PM (3 years ago)

Author:

fhourdin

Message:

Nouvelle option des programme iotd*
Possibilite de controler les variables de sorties via
iotd.def

Location:

LMDZ6/trunk/libf/phylmd

Files:

• iotd.h (modified) (1 diff)
• iotd_ecrit.F90 (modified) (3 diffs)
• iotd_ini.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/iotd.h

@@ -13 +13 @@
       INTEGER dim_coord(4)
       real iotd_ts,iotd_t0
+      integer :: n_names_iotd_def
+      character*20, dimension(200) :: names_iotd_def
+      character*20 :: un_nom
 
-      common/ecritd_c/imax,jmax,lmax,nid,dim_coord,iotd_t0,iotd_ts
-!$OMP THREADPRIVATE(/ecritd_c/)
+      common/iotd_ca/imax,jmax,lmax,nid,dim_coord,iotd_t0,iotd_ts
+      common/iotd_cb/n_names_iotd_def,names_iotd_def
+!$OMP THREADPRIVATE(/iotd_ca/)
+!$OMP THREADPRIVATE(/iotd_cb/)

LMDZ6/trunk/libf/phylmd/iotd_ecrit.F90

@@ -63 +63 @@
       
 
+
+       if (n_names_iotd_def>0 .and..not.any(names_iotd_def==nom)) return
 !***************************************************************
 ! Initialisation of 'firstnom' and create/open the "diagfi.nc" NetCDF file
@@ -94 +96 @@
 !!          print*,'ntime ',ntime
            date=iotd_t0+ntime*iotd_ts
-           print*,'iotd_ecrit ',iotd_ts,ntime, date
+           !print*,'iotd_ecrit ',iotd_ts,ntime, date
 !          date= float (zitau +1)/float (day_step)
 
@@ -160 +162 @@
           ierr = NF_REDEF (nid)
           ierr = NF_DEF_VAR(nid,nom,NF_FLOAT,ndim,dim_cc,varid)
-          print*,'DEF ',nom,nid,varid
+          !print*,'DEF ',nom,nid,varid
           ierr = NF_ENDDEF(nid)
       else
          ierr= NF_INQ_VARID(nid,nom,varid)
-          print*,'INQ ',nom,nid,varid
+          !print*,'INQ ',nom,nid,varid
 ! Commandes pour recuperer automatiquement les coordonnees
 !             ierr= NF_INQ_DIMID(nid,"longitude",id(1))

LMDZ6/trunk/libf/phylmd/iotd_ini.F90

@@ -39 +39 @@
       CHARACTER*11 date0b
 
-      INTEGER ierr
+      INTEGER :: ierr
 
       integer :: nvarid
@@ -60 +60 @@
 
 
-!-----------------------------------------------------------------------
+      !-----------------------------------------------------------------------
+      ! Possibilité de spécifier une liste de variables à sortir
+      ! dans iotd.def
+      ! Si iotd.def existe et est non vide,
+      ! seules les variables faisant à la fois l'objet d'un call iotd_ecrit
+      ! et étant spécifiées dans iotd.def sont sorties.
+      ! Sinon, toutes les variables faisant l'objet d'un call iotd_ecrit
+      ! sont sorties
+      !-----------------------------------------------------------------------
+      n_names_iotd_def=0
+      open(99,file='iotd.def',form='formatted',status='old',iostat=ierr)
+         if ( ierr.eq.0 ) then
+            ierr=0
+            do while (ierr==0)
+               read(99,*,iostat=ierr) un_nom
+               if (ierr==0) then
+                  n_names_iotd_def=n_names_iotd_def+1
+                  names_iotd_def(n_names_iotd_def)=un_nom
+               endif
+            enddo
+         endif
+         print*,n_names_iotd_def,names_iotd_def(1:n_names_iotd_def)
+      close(99)