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

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/phylmd/phyetat0.F90 (modified) (2 diffs)

LMDZ6/branches/LMDZ-tracers/libf/phylmd/phyetat0.F90

@@ -575 +575 @@
 USE dimphy, only: klon 
 USE print_control_mod, ONLY: lunout
-
+USE readTracFiles_mod, ONLY: known_phases, old_phases, nphases, phases_sep
+USE strings_mod,       ONLY: strIdx
 IMPLICIT NONE
 
@@ -588 +589 @@
 ! Local variables
 LOGICAL found
+INTEGER :: ix, ip
+CHARACTER(LEN=256) :: oldname
 
    CALL get_field(name, field, found)
+   ix = strIdx([('trs_H2O'//phases_sep//known_phases(ip:ip), ip=1, nphases)], name)
+   IF (.NOT. found .AND. ix /= 0) THEN           !--- Old file, water: H2Ov/l/i instead of H2O_g/_l/_s
+      oldname = 'trs_H2O'//old_phases(ix:ix)
+      CALL get_field(oldname, field, found)
+   END IF
    IF (.NOT. found) THEN
      WRITE(lunout,*) "phyetat0: Le champ <",TRIM(name),"> est absent"
      WRITE(lunout,*) "Depart legerement fausse. Mais je continue"
      field(:,:)=default
+   ELSE IF(ix /= 0) THEN
+     WRITE(lunout,*)"phyetat0: Le champ <",TRIM(name),"> est absent => initialisation a <"//TRIM(oldname)//">"
    ENDIF
    WRITE(lunout,*) TRIM(name), descr, MINval(field),MAXval(field)