Changeset 3957 for LMDZ6/branches/LMDZ-tracers/libf/dyn3d

Timestamp:

Jul 11, 2021, 11:39:01 PM (3 years ago)

Author:

dcugnet

Message:

In readTracFiles: the separator between the tracer name and its phase is no longer hardcoded and equal to "-",
but is a parameter ("phases_sep") which default value is "_".
Few more fixes.

File:

• LMDZ6/branches/LMDZ-tracers/libf/dyn3d/dynetat0.F90 (moved) (moved from LMDZ6/branches/LMDZ-tracers/libf/dyn3d/dynetat0.f90) (11 diffs)

LMDZ6/branches/LMDZ-tracers/libf/dyn3d/dynetat0.F90

@@ -6 +6 @@
 ! Purpose: Initial state reading.
 !-------------------------------------------------------------------------------
+  USE readTracFiles_mod, ONLY: known_phases, old_phases, nphases, phases_sep
   USE infotrac,    ONLY: nqtot, niso, tracers, iTraPha, tnat, alpha_ideal, tra
   USE netcdf,      ONLY: NF90_OPEN,  NF90_NOWRITE, NF90_INQ_VARID, NF90_NoErr, &
@@ -17 +18 @@
   USE serre_mod, ONLY: clon, clat, grossismx, grossismy
   USE temps_mod, ONLY: annee_ref, day_ini, day_ref, itau_dyn, start_time
-  USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
+  USE ener_mod,  ONLY: etot0, ptot0, ztot0, stot0, ang0
 
   IMPLICIT NONE
@@ -37 +38 @@
 !===============================================================================
 ! Local variables:
-  CHARACTER(LEN=256) :: sdum, var, modname, oldH2O(3), newH2O(3)
+  CHARACTER(LEN=256) :: sdum, var, modname, oldH2O
   INTEGER, PARAMETER :: length=100
-  INTEGER :: iq, fID, vID, idecal, ix!, iml, jml, lml, nqt
+  INTEGER :: iq, fID, vID, idecal, ix, ip
   REAL    :: time, tab_cntrl(length)               !--- RUN PARAMS TABLE
- TYPE(tra), POINTER :: tr
+  TYPE(tra), POINTER :: tr
 !-------------------------------------------------------------------------------
   modname="dynetat0"
-  oldH2O=['H2Ov ','H2Ol ','H2Oi ']
-  newH2O=['H2O-g','H2O-l','H2O-s']
 
 !--- Initial state file opening
@@ -107 +106 @@
   CALL check_dim(jm,jjm,'jm','jm')
   CALL check_dim(lllm,llm,'lm','lllm')
-  CALL get_var1("rlonu",rlonu)
-  CALL get_var1("rlatu",rlatu)
-  CALL get_var1("rlonv",rlonv)
-  CALL get_var1("rlatv",rlatv)
-  CALL get_var2("cu"   ,cu)
-  CALL get_var2("cv"   ,cv)
-  CALL get_var2("aire" ,aire)
+  CALL get_var1("rlonu", rlonu)
+  CALL get_var1("rlatu", rlatu)
+  CALL get_var1("rlonv", rlonv)
+  CALL get_var1("rlatv", rlatv)
+  CALL get_var2("cu",       cu)
+  CALL get_var2("cv",       cv)
+  CALL get_var2("aire",   aire)
   var="temps"
   IF(NF90_INQ_VARID(fID,var,vID)/=NF90_NoErr) THEN
@@ -121 +120 @@
   END IF
   CALL err(NF90_GET_VAR(fID,vID,time),"get",var)
-  CALL get_var2("phisinit",phis)
-  CALL get_var3("ucov",ucov)
-  CALL get_var3("vcov",vcov)
-  CALL get_var3("teta",teta)
-  CALL get_var3("masse",masse)
-  CALL get_var2("ps",ps)
+  CALL get_var2("phisinit", phis)
+  CALL get_var3("ucov",     ucov)
+  CALL get_var3("vcov",     vcov)
+  CALL get_var3("teta",     teta)
+  CALL get_var3("masse",   masse)
+  CALL get_var2("ps",         ps)
 
 !--- Tracers
@@ -132 +131 @@
     tr => tracers(iq)
     var = tr%name
-    IF(NF90_INQ_VARID(fID,var,vID)==NF90_NoErr) THEN
-      CALL err(NF90_GET_VAR(fID,vID,q(:,:,:,iq)),"get",var); CYCLE
+    ix = strIdx([('H2O'//phases_sep//known_phases(ip:ip), ip=1, nphases)], var)
+    IF(NF90_INQ_VARID(fID, var, vID) == NF90_NoErr) THEN
+      CALL get_var3(var, q(:,:,:,iq))
 #ifdef INCA
-    ELSE IF(var == "O3") THEN          !--- INCA and O3 missing: take OX instead
-      WRITE(lunout,*) 'Tracer O3 is missing - it is initialized to OX'
-        IF(NF90_INQ_VARID(fID,"OX",vID) == NF90_NoErr) THEN
-           CALL err(NF90_GET_VAR(fID,vID,q(:,:,:,iq)),"get",var); CYCLE
-        END IF
+    ELSE IF(NF90_INQ_VARID(fID, 'OX', vID) == NF90_NoErr .AND. var == 'O3') THEN
+      WRITE(lunout,*)TRIM(modname)//': Tracer <O3> is missing => initialized to OX'
+      CALL get_var3('OX', q(:,:,:,iq))
 #endif
-    ELSE                               !--- Old file, water: H2Ov/l/i instead of H2O-g/-l/-s
-      ix = strIdx(newH2O, var)         !--- Current tracer is water (new name) ?
-      IF(ix /= 0) THEN                 !--- Then read the field, using the old name.
-        IF(NF90_INQ_VARID(fID,oldH2O(ix),vID) == NF90_NoErr) THEN
-           CALL err(NF90_GET_VAR(fID,vID,q(:,:,:,iq)),"get",var); CYCLE
-        END IF
+    ELSE IF(ix /= 0) THEN              !--- Old file, water: H2Ov/l/i instead of H2O_g/_l/_s
+      oldH2O = 'H2O'//old_phases(ix:ix)
+      IF(NF90_INQ_VARID(fID, oldH2O, vID) == NF90_NoErr) THEN
+        WRITE(lunout,*)TRIM(modname)//': Tracer <'//TRIM(var)//'> is missing => initialized to '//TRIM(oldH2O)
+        CALL get_var3(oldH2O, q(:,:,:,iq))
+      END IF
+    ELSE
+      WRITE(lunout,*)TRIM(modname)//': Tracer <'//TRIM(var)//'> is missing => initialized to zero'
+      q(:,:,:,iq)=0.
+      !--- CRisi: for isotopes, theoretical initialization using very simplified Rayleigh distillation law
+      IF(niso > 0 .AND. tr%iso_num > 0) THEN
+        IF(tr%iso_zon == 0) q(:,:,:,iq) = q(:,:,:,tr%iprnt)         *        tnat(tr%iso_num)  &
+                                        *(q(:,:,:,tr%iprnt)/30.e-3)**(alpha_ideal(tr%iso_num)-1)
+        IF(tr%iso_zon == 1) q(:,:,:,iq) = q(:,:,:,iTraPha(tr%iso_num,tr%iso_pha))
       END IF
     END IF
-    WRITE(lunout,*)TRIM(modname)//": Tracer <"//TRIM(var)//"> is missing"
-    WRITE(lunout,*)"         It is hence initialized to zero"
-    q(:,:,:,iq)=0.
-    !--- CRisi: for isotops, theoretical initialization using very simplified
-    !           Rayleigh distillation las.
-    IF(niso > 0 .AND. tr%iso_num > 0) THEN
-      IF(tr%iso_zon == 0) q(:,:,:,iq) = q(:,:,:,tr%iprnt)         *        tnat(tr%iso_num)  &
-                                      *(q(:,:,:,tr%iprnt)/30.e-3)**(alpha_ideal(tr%iso_num)-1)
-      IF(tr%iso_zon == 1) q(:,:,:,iq) = q(:,:,:,iTraPha(tr%iso_num,tr%iso_pha))
-    END IF
   END DO
-
   CALL err(NF90_CLOSE(fID),"close",fichnom)
   day_ini=day_ini+INT(time)
@@ -176 +171 @@
     s1='value of '//TRIM(str1)//' ='
     s2=' read in starting file differs from parametrized '//TRIM(str2)//' ='
-    WRITE(sdum,'(10x,a,i4,2x,a,i4)'),TRIM(ADJUSTL(s1)),n1,TRIM(ADJUSTL(s2)),n2
+    WRITE(sdum,'(10x,a,i4,2x,a,i4)'),TRIM(s1),n1,TRIM(s2),n2
     CALL ABORT_gcm(TRIM(modname),TRIM(sdum),1)
   END IF
@@ -182 +177 @@
 
 
-SUBROUTINE get_var1(var,v)
+SUBROUTINE get_var1(var, v)
   CHARACTER(LEN=*), INTENT(IN)  :: var
   REAL,             INTENT(OUT) :: v(:)
@@ -190 +185 @@
 
 
-SUBROUTINE get_var2(var,v)
+SUBROUTINE get_var2(var, v)
   CHARACTER(LEN=*), INTENT(IN)  :: var
   REAL,             INTENT(OUT) :: v(:,:)
@@ -198 +193 @@
 
 
-SUBROUTINE get_var3(var,v)
+SUBROUTINE get_var3(var, v)
   CHARACTER(LEN=*), INTENT(IN)  :: var
   REAL,             INTENT(OUT) :: v(:,:,:)
@@ -206 +201 @@
 
 
-SUBROUTINE err(ierr,typ,nam)
+SUBROUTINE err(ierr, typ, nam)
   INTEGER,          INTENT(IN) :: ierr   !--- NetCDF ERROR CODE
   CHARACTER(LEN=*), INTENT(IN) :: typ    !--- TYPE OF OPERATION