Changeset 4146 for LMDZ6/trunk/libf/dyn3dmem/gcm.F90

Timestamp:

May 10, 2022, 6:50:43 PM (2 years ago)

Author:

Laurent Fairhead

Message:

Source code and xml files needed for XIOS output in the LMDZ LonLat? dynamical core.
One flag controls everything: ok_dyn_xios. Parameters controlling outputs need to be put in the xml
files, nothing is set in the code
Work to follow on the dynzon file and the modipsl/libIGCM framework.
ok_dyn_xios = n should not change anything

File:

• LMDZ6/trunk/libf/dyn3dmem/gcm.F90 (modified) (7 diffs, 1 prop)

LMDZ6/trunk/libf/dyn3dmem/gcm.F90

@@ -1 @@
-! $Id: $
+!
+! $Id$
+!
 
 PROGRAM gcm
@@ -26 +28 @@
                        itau_dyn,itau_phy,day_ini,jD_ref,jH_ref,day_end, &
                        dt,hour_ini,itaufin
+#ifdef CPP_XIOS
+  USE mod_xios_dyn3dmem, ONLY: xios_dyn3dmem_init
+#endif
 
   IMPLICIT NONE
@@ -109 +114 @@
   INTEGER :: an, mois, jour
   REAL :: heure
+  ! needed for xios interface
+  character (len=10) :: xios_cal_type
+  INTEGER :: anref, moisref, jourref
+  REAL :: heureref
+  
 
 
@@ -174 +184 @@
      call ioconf_calendar('360d')
      write(lunout,*)'CALENDRIER CHOISI: Terrestre a 360 jours/an'
+     xios_cal_type='d360'
   else if (calend == 'earth_365d') then
      call ioconf_calendar('noleap')
      write(lunout,*)'CALENDRIER CHOISI: Terrestre a 365 jours/an'
+     xios_cal_type='noleap'
   else if (calend == 'gregorian') then
      call ioconf_calendar('gregorian')
      write(lunout,*)'CALENDRIER CHOISI: Terrestre bissextile'
+     xios_cal_type='gregorian'
   else
      abort_message = 'Mauvais choix de calendrier'
@@ -331 +344 @@
   write(lunout,*)'annee_ref, mois, day_ref, heure, jD_ref'
   write(lunout,*)annee_ref, mois, day_ref, heure, jD_ref
-  call ju2ymds(jD_ref+jH_ref,an, mois, jour, heure)
+  call ju2ymds(jD_ref+jH_ref,anref, moisref, jourref, heureref)
   write(lunout,*)'jD_ref+jH_ref,an, mois, jour, heure'
-  write(lunout,*)jD_ref+jH_ref,an, mois, jour, heure
+  write(lunout,*)jD_ref+jH_ref,anref, moisref, jourref, heureref
 #else
   ! Ehouarn: we still need to define JD_ref and JH_ref
@@ -421 +434 @@
   ecripar = .TRUE.
 
+#define CPP_IOIPSL
 #ifdef CPP_IOIPSL
   time_step = zdtvr
@@ -440 +454 @@
   dtav = iperiod*dtvr/daysec
 #endif
+#undef CPP_IOIPSL
+
+! setting up DYN3D/XIOS inerface
+#ifdef CPP_XIOS
+   if (ok_dyn_xios) then
+     CALL xios_dyn3dmem_init(xios_cal_type, anref, moisref, jourref,heureref, an,   &
+          mois, jour, heure, zdtvr)
+  endif
+#endif
+
   ! #endif of #ifdef CPP_IOIPSL
   !