Changeset 4539

Timestamp:

May 17, 2023, 12:53:42 AM (12 months ago)

Author:

fhourdin

Message:

Modification for the new physiqex programme

Location:

LMDZ6/trunk/libf/phylmd

Files:

• output_physiqex_mod.F90 (added)
• physiqex_mod.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/physiqex_mod.F90

@@ -17 +17 @@
       USE geometry_mod, only : latitude
 !      USE comcstphy, only : rg
-      USE iophy, only : histbeg_phy,histwrite_phy
-      USE ioipsl, only : getin,histvert,histdef,histend,ymds2ju
-      USE mod_phys_lmdz_para, only : jj_nb
+      USE ioipsl, only : ymds2ju
       USE phys_state_var_mod, only : phys_state_var_init
-      USE mod_grid_phy_lmdz, ONLY: nbp_lon,nbp_lat
       USE phyetat0_mod, only: phyetat0
-      USE iophy, ONLY : init_iophy_new
-      USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
-
-
-#ifdef CPP_XIOS
-      USE xios, ONLY: xios_update_calendar
-      USE wxios, only: wxios_add_vaxis, wxios_set_cal, wxios_closedef
-      USE iophy, ONLY: histwrite_phy
-#endif
+      USE output_physiqex_mod, ONLY: output_physiqex
 
       IMPLICIT none
@@ -69 +58 @@
 
 
-integer,save :: itau=0 ! counter to count number of calls to physics
-!$OMP THREADPRIVATE(itau)
 real :: temp_newton(klon,klev)
 integer :: k
@@ -82 +69 @@
 integer :: itau0
 real :: zjulian
-real :: dtime
-integer :: nhori ! horizontal coordinate ID
-integer,save :: nid_hist ! output file ID
-!$OMP THREADPRIVATE(nid_hist)
-integer :: zvertid ! vertical coordinate ID
-integer,save :: iwrite_phys ! output every iwrite_phys physics step
-!$OMP THREADPRIVATE(iwrite_phys)
-integer,save :: iwrite_phys_omp ! intermediate variable to read iwrite_phys
-                                ! (must be shared by all threads)
-real :: t_ops ! frequency of the IOIPSL operations (eg average over...)
-real :: t_wrt ! frequency of the IOIPSL outputs
-
 
 
@@ -108 +83 @@
   call phys_state_var_init(1) ! some initializations, required before calling phyetat0
   call phyetat0("startphy.nc", clesphy0, tabcntr0)
-  CALL init_iophy_new(latitude_deg,longitude_deg)
 
 ! Initialize outputs:
   itau0=0
-!$OMP MASTER
-  iwrite_phys_omp=1 !default: output every physics timestep
-  ! NB: getin() is not threadsafe; only one thread should call it.
-  call getin("iwrite_phys",iwrite_phys_omp)
-!$OMP END MASTER
-!$OMP BARRIER
-  iwrite_phys=iwrite_phys_omp
-  t_ops=pdtphys*iwrite_phys ! frequency of the IOIPSL operation
-  t_wrt=pdtphys*iwrite_phys ! frequency of the outputs in the file
   ! compute zjulian for annee0=1979 and month=1 dayref=1 and hour=0.0
   !CALL ymds2ju(annee0, month, dayref, hour, zjulian)
   call ymds2ju(1979, 1, 1, 0.0, zjulian)
-  dtime=pdtphys
 
 #ifndef CPP_IOIPSL_NO_OUTPUT
   ! Initialize IOIPSL output file
-  call histbeg_phy("histins.nc",itau0,zjulian,dtime,nhori,nid_hist)
 #endif
 
@@ -137 +100 @@
 !------------------------------------------------------------
 
-! increment local time counter itau
-itau=itau+1
 
 ! set all tendencies to zero
@@ -166 +127 @@
 !------------------------------------------------------------
 
-print*,'PHYDEV: itau=',itau
 
-
-if(debut)then
-!$OMP MASTER
-#ifndef CPP_IOIPSL_NO_OUTPUT 
-! IOIPSL
-  ! define vertical coordinate
-  call histvert(nid_hist,"presnivs","Vertical levels","Pa",klev, &
-                presnivs,zvertid,'down')
-  ! define variables which will be written in "histins.nc" file
-  call histdef(nid_hist,'temperature','Atmospheric temperature','K', &
-               nbp_lon,jj_nb,nhori,klev,1,klev,zvertid,32, &
-               'inst(X)',t_ops,t_wrt)
-  call histdef(nid_hist,'u','Eastward Zonal Wind','m/s', &
-               nbp_lon,jj_nb,nhori,klev,1,klev,zvertid,32, &
-               'inst(X)',t_ops,t_wrt)
-  call histdef(nid_hist,'v','Northward Meridional Wind','m/s', &
-               nbp_lon,jj_nb,nhori,klev,1,klev,zvertid,32, &
-               'inst(X)',t_ops,t_wrt)
-  call histdef(nid_hist,'ps','Surface Pressure','Pa', &
-               nbp_lon,jj_nb,nhori,1,1,1,zvertid,32, &
-               'inst(X)',t_ops,t_wrt)
-  ! end definition sequence
-  call histend(nid_hist)
-#endif
-
-#ifdef CPP_XIOS
-!XIOS
-    ! Declare available vertical axes to be used in output files:    
-    CALL wxios_add_vaxis("presnivs", klev, presnivs)
-
-    ! Declare calendar and time step
-    CALL wxios_set_cal(dtime,"earth_360d",1,1,1,0.0,1,1,1,0.0)
-    
-    !Finalize the context:
-    CALL wxios_closedef()
-#endif
-!$OMP END MASTER
-!$OMP BARRIER
-endif
-! write some outputs:
-! IOIPSL
-#ifndef CPP_IOIPSL_NO_OUTPUT 
-if (modulo(itau,iwrite_phys)==0) then
-  call histwrite_phy(nid_hist,.false.,"temperature",itau,t)
-  call histwrite_phy(nid_hist,.false.,"u",itau,u)
-  call histwrite_phy(nid_hist,.false.,"v",itau,v)
-  call histwrite_phy(nid_hist,.false.,"ps",itau,paprs(:,1))
-endif
-#endif
-
-!XIOS
-#ifdef CPP_XIOS
-!$OMP MASTER
-    !Increment XIOS time
-    CALL xios_update_calendar(itau)
-!$OMP END MASTER
-!$OMP BARRIER
-
-    !Send fields to XIOS: (NB these fields must also be defined as
-    ! <field id="..." /> in iodef.xml to be correctly used
-    CALL histwrite_phy("temperature",t)
-    CALL histwrite_phy("temp_newton",temp_newton)
-    CALL histwrite_phy("u",u)
-    CALL histwrite_phy("v",v)
-    CALL histwrite_phy("ps",paprs(:,1))
-#endif
+call output_physiqex(debut,zjulian,pdtphys,presnivs,paprs,u,v,t)
 
 ! if lastcall, then it is time to write "restartphy.nc" file
@@ -240 +135 @@
 endif
 
+
 end subroutine physiqex