Changeset 4727 for LMDZ6/branches/LMDZ_ECRad/libf/dyn3d

Timestamp:

Oct 19, 2023, 4:02:57 PM (21 months ago)

Author:

idelkadi

Message:

Merged trunk changes -r4488:4726 LMDZ_ECRad branch

Location:

LMDZ6/branches/LMDZ_ECRad

Files:

• . (modified) (1 prop)
• libf/dyn3d/abort_gcm.F (modified) (2 diffs)
• libf/dyn3d/conf_gcm.F90 (modified) (2 diffs)
• libf/dyn3d/dynetat0.F90 (modified) (1 diff)
• libf/dyn3d/gcm.F90 (modified) (2 diffs)
• libf/dyn3d/replay3d.F90 (modified) (1 diff)

LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/abort_gcm.F

@@ -12 +12 @@
       USE ioipsl_getincom
 #endif
-
-#ifdef CPP_XIOS
-    ! ug Pour les sorties XIOS
+ ! ug Pour les sorties XIOS
       USE wxios
-#endif
 
 #include "iniprint.h"
@@ -34 +31 @@
       write(lunout,*) 'in abort_gcm'
 
-#ifdef CPP_XIOS
-    !Fermeture propre de XIOS
-      CALL wxios_close()
-#endif
+      IF (using_xios) THEN
+!Fermeture propre de XIOS
+        CALL wxios_close()
+      ENDIF
 
 #ifdef CPP_IOIPSL

LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/conf_gcm.F90

@@ -15 +15 @@
   USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &
                           iflag_top_bound, mode_top_bound, tau_top_bound, &
-                          ngroup
+                          ngroup, maxlatfilter
   USE logic_mod, ONLY: fxyhypb, iflag_phys, ok_etat0, ok_gradsfile, &
                        ok_guide, ok_limit, ok_strato, purmats, read_start, &
@@ -310 +310 @@
   CALL getin('dissip_deltaz',dissip_deltaz )
   CALL getin('dissip_zref',dissip_zref )
+
+  ! maxlatfilter
+  maxlatfilter=-1.0
+  CALL getin('maxlatfilter',maxlatfilter)
+  if (maxlatfilter > 90) &
+       call abort_gcm("conf_gcm", 'maxlatfilter should be <=90', 1)
+
 
   ! ngroup

LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/dynetat0.F90

@@ -161 +161 @@
       ELSE
          CALL msg('Tracer <'//TRIM(var)//'> is missing => initialized to its parent isotope concentration.', modname)
-         q(:,:,:,iq) = q(:,:,:,iqIsoPha(iName,iPhase))
+         ! Camille 9 mars 2023: attention!! seuls les tags qui correspondent à
+         ! izone=izone_init (définie dans isotrac_mod) sont initialisés comme
+         ! les parents. Sinon, c'est nul.
+         ! j'ai fait ça en attendant, mais il faudrait initialiser proprement en
+         ! remplacant 1 par izone_init dans la ligne qui suit.
+         IF(tracers(iq)%iso_iZone == 1) THEN
+           q(:,:,:,iq) = q(:,:,:,iqIsoPha(iName,iPhase))
+         ELSE
+           q(:,:,:,iq) = 0.
+         END IF
       END IF
     !--------------------------------------------------------------------------------------------------------------------------

LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/gcm.F90

@@ -14 +14 @@
 
 
-#ifdef CPP_XIOS
-  ! ug Pour les sorties XIOS
+! ug Pour les sorties XIOS
   USE wxios
-#endif
 
   USE filtreg_mod
@@ -156 +154 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-#ifdef CPP_XIOS
-  CALL wxios_init("LMDZ")
-#endif
+  IF (using_xios) THEN
+    CALL wxios_init("LMDZ")
+  ENDIF
 
 

LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/replay3d.F90

@@ -182 +182 @@
 ! Initialisation de la parametrisation
 !---------------------------------------------------------------------
-      call get_ini_module
+      call call_ini_replay
 
 !---------------------------------------------------------------------