Changeset 4368 for LMDZ6/branches/Ocean_skin/libf/dyn3dmem/conf_gcm.F90

Timestamp:

Dec 6, 2022, 12:01:16 AM (18 months ago)

Author:

lguez

Message:

Sync latest trunk changes to Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/dyn3dmem/conf_gcm.F90 (modified) (19 diffs)

LMDZ6/branches/Ocean_skin/libf/dyn3dmem/conf_gcm.F90

@@ -6 +6 @@
   USE control_mod
 #ifdef CPP_IOIPSL
-  use IOIPSL
+  USE IOIPSL
 #else
   ! if not using IOIPSL, we still need to use (a local version of) getin
-  use ioipsl_getincom
+  USE ioipsl_getincom
 #endif
-  use misc_mod
-  use mod_filtre_fft, ONLY : use_filtre_fft
-  use mod_filtre_fft_loc, ONLY : use_filtre_fft_loc=>use_filtre_fft
-  use mod_hallo, ONLY : use_mpi_alloc
-  USE infotrac, ONLY : type_trac
-  use assert_m, only: assert
+  USE misc_mod
+  USE mod_filtre_fft, ONLY: use_filtre_fft
+  USE mod_filtre_fft_loc, ONLY: use_filtre_fft_loc=>use_filtre_fft
+  USE mod_hallo, ONLY: use_mpi_alloc
+  USE infotrac, ONLY: type_trac
+  USE assert_m, ONLY: assert
   USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &
                           iflag_top_bound, mode_top_bound, tau_top_bound, &
@@ -55 +55 @@
   LOGICAL  fxyhypbb, ysinuss
   INTEGER i
-  character(len=*),parameter :: modname="conf_gcm"
-  character (len=80) :: abort_message
+  CHARACTER(len=*), PARAMETER :: modname="conf_gcm"
+  CHARACTER(len=80) :: abort_message
 #ifdef CPP_OMP
-  integer,external :: OMP_GET_NUM_THREADS
+  INTEGER, EXTERNAL :: OMP_GET_NUM_THREADS
 #endif
 
@@ -102 +102 @@
 
   adjust=.false.
-  call getin('adjust',adjust)
+  CALL getin('adjust',adjust)
 
 #ifdef CPP_OMP
   ! adjust=y not implemented in case of OpenMP threads...
   !$OMP PARALLEL
-  if ((OMP_GET_NUM_THREADS()>1).and.adjust) then
+  IF ((OMP_GET_NUM_THREADS()>1).and.adjust) then
      write(lunout,*)'conf_gcm: Error, adjust should be set to n' &
           ,' when running with OpenMP threads'
      abort_message = 'Wrong value for adjust'
-     call abort_gcm(modname,abort_message,1)
-  endif
+     CALL abort_gcm(modname,abort_message,1)
+  ENDIF
   !$OMP END PARALLEL         
 #endif
 
   itaumax=0
-  call getin('itaumax',itaumax);
-  if (itaumax<=0) itaumax=HUGE(itaumax)
+  CALL getin('itaumax',itaumax);
+  IF (itaumax<=0) itaumax=HUGE(itaumax)
 
   !Config  Key  = prt_level
@@ -146 +146 @@
 ! initialize year_len for aquaplanets and 1D
   CALL getin('calend', calend)
-     if (calend == 'earth_360d') then
-        year_len=360
-      else if (calend == 'earth_365d') then
-        year_len=365
-      else if (calend == 'earth_366d') then
-        year_len=366
-      else 
-        year_len=1
-      endif
-  
+  IF (calend == 'earth_360d') THEN
+    year_len=360
+  ELSE IF (calend == 'earth_365d') THEN
+    year_len=365
+  ELSE IF (calend == 'earth_366d') THEN
+    year_len=366
+  ELSE
+    year_len=1
+  ENDIF
 
   !Config  Key  = dayref
@@ -345 +344 @@
   ngroup=3
   CALL getin('ngroup',ngroup)
-  if (mod(iim, 2**ngroup) /= 0) &
+  IF (mod(iim, 2**ngroup) /= 0) &
        call abort_gcm("conf_gcm", 'iim must be multiple of 2**ngroup', 1)
-  if (2**ngroup > jjm + 1) &
+  IF (2**ngroup > jjm + 1) &
        call abort_gcm("conf_gcm", '2**ngroup must be <= jjm + 1', 1)
 
@@ -357 +356 @@
   iflag_top_bound=1
   CALL getin('iflag_top_bound',iflag_top_bound)
+  IF (iflag_top_bound < 0 .or. iflag_top_bound > 2) &
+       call abort_gcm("conf_gcm", 'iflag_top_bound must be 0, 1 or 2', 1)
 
   ! mode_top_bound : fields towards which sponge relaxation will be done:
@@ -392 +393 @@
   CALL getin('ok_guide',ok_guide)
 
-  if (ok_guide .and. adjust) call abort_gcm("conf_gcm", &
+  IF (ok_guide .and. adjust) call abort_gcm("conf_gcm", &
        "adjust does not work with ok_guide", 1)
 
@@ -632 +633 @@
         WRITE(lunout,*)  &
              'only the file phystoke.nc will still be created ' 
-     END IF
+     ENDIF
 
      !Config  Key  = type_trac
@@ -644 +645 @@
      CALL getin('type_trac',type_trac)
 
-     !Config  Key  = config_inca
-     !Config  Desc = Choix de configuration de INCA
-     !Config  Def  = none
-     !Config  Help = Choix de configuration de INCA :
-     !Config         'none' = sans INCA
-     !Config         'chem' = INCA avec calcul de chemie
-     !Config         'aero' = INCA avec calcul des aerosols 
-     config_inca = 'none'
-     CALL getin('config_inca',config_inca)
-
      !Config  Key  = ok_dynzon 
      !Config  Desc = calcul et sortie des transports 
@@ -677 +668 @@
      ok_dyn_ave = .FALSE. 
      CALL getin('ok_dyn_ave',ok_dyn_ave) 
+
+     !Config  Key  = ok_dyn_xios
+     !Config  Desc = sorties moyennes dans la dynamique
+     !Config  Def  = n 
+     !Config  Help = 
+     !Config          
+     ok_dyn_xios = .FALSE. 
+     CALL getin('ok_dyn_xios',ok_dyn_xios) 
 
      write(lunout,*)' #########################################'
@@ -716 +715 @@
      write(lunout,*)' offline = ', offline
      write(lunout,*)' type_trac = ', type_trac
-     write(lunout,*)' config_inca = ', config_inca
      write(lunout,*)' ok_dynzon = ', ok_dynzon 
      write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins 
      write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave 
+     write(lunout,*)' ok_dyn_xios = ', ok_dyn_xios 
   else
      !Config  Key  = clon
@@ -754 +753 @@
 
      IF( grossismx.LT.1. )  THEN
-        write(lunout,*) &
-             'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
+        write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismx < 1 . *** '
         STOP
      ELSE
@@ -768 +766 @@
      ENDIF
 
-     write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
+     write(lunout,*) 'conf_gcm: alphax alphay ',alphax,alphay
 
      !    alphax et alphay sont les anciennes formulat. des grossissements
@@ -835 +833 @@
         WRITE(lunout,*)  &
              'only the file phystoke.nc will still be created ' 
-     END IF
+     ENDIF
 
      !Config  Key  = type_trac
@@ -847 +845 @@
      CALL getin('type_trac',type_trac)
 
-     !Config  Key  = config_inca
-     !Config  Desc = Choix de configuration de INCA
-     !Config  Def  = none
-     !Config  Help = Choix de configuration de INCA :
-     !Config         'none' = sans INCA
-     !Config         'chem' = INCA avec calcul de chemie
-     !Config         'aero' = INCA avec calcul des aerosols 
-     config_inca = 'none'
-     CALL getin('config_inca',config_inca)
-
      !Config  Key  = ok_dynzon 
      !Config  Desc = sortie des transports zonaux dans la dynamique
@@ -880 +868 @@
      ok_dyn_ave = .FALSE. 
      CALL getin('ok_dyn_ave',ok_dyn_ave) 
+
+     !Config  Key  = ok_dyn_xios
+     !Config  Desc = sorties moyennes dans la dynamique
+     !Config  Def  = n 
+     !Config  Help = 
+     !Config          
+     ok_dyn_xios = .FALSE. 
+     CALL getin('ok_dyn_xios',ok_dyn_xios) 
 
      !Config  Key  = use_filtre_fft
@@ -953 +949 @@
 
      write(lunout,*)' #########################################'
-     write(lunout,*)' Configuration des parametres de cel0' &
-          //'_limit: '
+     write(lunout,*)' Configuration des parametres de cel0_limit: '
      write(lunout,*)' planet_type = ', planet_type
      write(lunout,*)' calend = ', calend
@@ -990 +985 @@
      write(lunout,*)' offline = ', offline
      write(lunout,*)' type_trac = ', type_trac
-     write(lunout,*)' config_inca = ', config_inca
      write(lunout,*)' ok_dynzon = ', ok_dynzon
      write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins 
      write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave 
+     write(lunout,*)' ok_dyn_xios = ', ok_dyn_xios 
      write(lunout,*)' use_filtre_fft = ', use_filtre_fft
      write(lunout,*)' use_mpi_alloc = ', use_mpi_alloc
@@ -1000 +995 @@
      write(lunout,*)' ok_limit = ', ok_limit
      write(lunout,*)' ok_etat0 = ', ok_etat0
+     write(lunout,*)' ok_guide = ', ok_guide
      write(lunout,*)' read_orop = ', read_orop
-  end IF test_etatinit
+  ENDIF test_etatinit
 
 END SUBROUTINE conf_gcm