Changeset 4540

Timestamp:

May 17, 2023, 2:17:37 AM (12 months ago)

Author:

fhourdin

Message:

Modification sorties physiqex

File:

• LMDZ6/trunk/libf/phylmd/output_physiqex_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/output_physiqex_mod.F90

@@ -40 +40 @@
 !$OMP THREADPRIVATE(itau)
 
+integer, save :: ioex=1
 
 
@@ -46 +47 @@
 if(debut)then
 
-   CALL init_iophy_new(latitude_deg,longitude_deg)
+   call getin_p("iwrite_phys",iwrite_phys)
+   call getin_p("ioex",ioex)
 
-   dtime=pdtphys
-   itau=0
-   call histbeg_phy("histins.nc",itau,zjulian,dtime,nhori,nid_hist)
-   print*,'NNNNNNN ',nid_hist,debut
+   if ( ioex == 1 ) then
+      CALL iophys_ini(pdtphys)
+   else if ( ioex == 2 ) then
 
-   call getin_p("iwrite_phys",iwrite_phys)
-   print*,'NNNNNNN OK0'
-   t_ops=pdtphys*iwrite_phys ! frequency of the IOIPSL operation
-   t_wrt=pdtphys*iwrite_phys ! frequency of the outputs in the file
-   print*,'NNNNNNN OK1'
+      CALL init_iophy_new(latitude_deg,longitude_deg)
+      dtime=pdtphys
+      itau=0
+      call histbeg_phy("histins.nc",itau,zjulian,dtime,nhori,nid_hist)
+      print*,'NNNNNNN ',nid_hist,debut
+      print*,'NNNNNNN OK0'
+      t_ops=pdtphys*iwrite_phys ! frequency of the IOIPSL operation
+      t_wrt=pdtphys*iwrite_phys ! frequency of the outputs in the file
+      print*,'NNNNNNN OK1'
 
    !$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,'Temp','Atmospheric temperature','K', &
-                  nbp_lon,jj_nb,nhori,klev,1,klev,zvertid,32, &
-                  'inst(X)',t_ops,t_wrt)
-     print*,'NNNNNNN OK2a',nid_hist,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)
-     print*,'NNNNNNN OK2b',nid_hist,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)
-     print*,'NNNNNNN OK2c',nid_hist,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
-     print*,'NNNNNNN OK2',nid_hist,t_ops,t_wrt
-     call histend(nid_hist)
-     print*,'NNNNNNN OK3'
+       ! 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,'Temp','Atmospheric temperature','K', &
+                    nbp_lon,jj_nb,nhori,klev,1,klev,zvertid,32, &
+                    'inst(X)',t_ops,t_wrt)
+       print*,'NNNNNNN OK2a',nid_hist,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)
+       print*,'NNNNNNN OK2b',nid_hist,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)
+       print*,'NNNNNNN OK2c',nid_hist,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
+       print*,'NNNNNNN OK2',nid_hist,t_ops,t_wrt
+       call histend(nid_hist)
+       print*,'NNNNNNN OK3'
 #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()
+      !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
+      !$OMP END MASTER
+      !$OMP BARRIER
+   endif
 
 endif
@@ -112 +118 @@
 #ifndef CPP_IOIPSL_NO_OUTPUT 
 if (modulo(itau,iwrite_phys)==0) then
-  call histwrite_phy(nid_hist,.false.,"Temp",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))
-!$OMP MASTER
+  if ( ioex == 1 ) then
+     call iophys_ecrit('temp',klev,'Temperature','K',t)
+     call iophys_ecrit('u',klev,'zonal wind','m/s',t)
+     call iophys_ecrit('v',klev,'meridinal wind','m/s',t)
+     call iophys_ecrit('ps',1,'Surface pressure','Pa',paprs(:,1))
+  else if ( ioex == 2 ) then
+     call histwrite_phy(nid_hist,.false.,"Temp",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))
+     !$OMP MASTER
      CALL histsync(nid_hist)
-!$OMP END MASTER
+     !$OMP END MASTER
+  endif
 endif
 #endif