Changeset 5101 for LMDZ6/branches/Amaury_dev/libf/phydev

Timestamp:

Jul 23, 2024, 8:22:55 AM (6 months ago)

Author:

abarral

Message:

Handle DEBUG_IO in lmdz_cppkeys_wrapper.F90
Transform some files .F -> .[fF]90
[ne compile pas à cause de writefield_u non défini - en attente de réponse Laurent]

Location:

LMDZ6/branches/Amaury_dev/libf/phydev

Files:

• iophy.F90 (modified) (2 diffs)
• phyetat0.F90 (modified) (1 diff)
• phyredem.F90 (modified) (1 diff)
• phys_state_var_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/libf/phydev/iophy.F90

@@ -103 +103 @@
     
 #ifndef CPP_IOIPSL_NO_OUTPUT
-    call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
+    CALL flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
                       'APPLE',phys_domain_id)
 #endif
@@ -157 +157 @@
 !$OMP MASTER    
     if (is_sequential) then
-      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
+      CALL histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
                    1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
     else
-      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
+      CALL histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
                    1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
     endif

LMDZ6/branches/Amaury_dev/libf/phydev/phyetat0.F90

@@ -7 +7 @@
 
   USE dimphy, only: klon
-  USE iostart, ONLY : open_startphy,get_field,close_startphy
-  USE iophy, ONLY : init_iophy_new
-  USE geometry_mod, ONLY : longitude_deg, latitude_deg
+  USE iostart, ONLY: open_startphy,get_field,close_startphy
+  USE iophy, ONLY: init_iophy_new
+  USE geometry_mod, ONLY: longitude_deg, latitude_deg
 
   IMPLICIT NONE

LMDZ6/branches/Amaury_dev/libf/phydev/phyredem.F90

@@ -4 +4 @@
 SUBROUTINE phyredem (fichnom)
 
-  USE geometry_mod, ONLY : longitude_deg, latitude_deg
+  USE geometry_mod, ONLY: longitude_deg, latitude_deg
   USE iostart, ONLY: open_restartphy, close_restartphy, enddef_restartphy, put_field, put_var
 

LMDZ6/branches/Amaury_dev/libf/phydev/phys_state_var_mod.F90

@@ -7 +7 @@
 !======================================================================
 
-!USE dimphy, only : klon
+!USE dimphy, ONLY: klon
  
 
@@ -17 +17 @@
 !======================================================================
   SUBROUTINE phys_state_var_init()
-!  use dimphy, only : klon
+!  use dimphy, ONLY: klon
 
 !  if (.not.allocated(rlat)) then
@@ -29 +29 @@
 !======================================================================
   SUBROUTINE phys_state_var_end
-!  use dimphy, only : klon
+!  use dimphy, ONLY: klon
 
 !  deallocate(rlat,rlon)

LMDZ6/branches/Amaury_dev/libf/phydev/physiq_mod.F90

@@ -13 +13 @@
               d_u, d_v, d_t, d_qx, d_ps)
 
-      USE dimphy, only : klon,klev
-      USE infotrac_phy, only : nqtot
-      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 phys_state_var_mod, only : phys_state_var_init
+      USE dimphy, ONLY: klon,klev
+      USE infotrac_phy, ONLY: nqtot
+      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 phys_state_var_mod, ONLY: phys_state_var_init
       USE mod_grid_phy_lmdz, ONLY: nbp_lon,nbp_lat
 
@@ -33 +33 @@
       integer,intent(in) :: nlon ! number of atmospheric colums
       integer,intent(in) :: nlev ! number of vertical levels (should be =klev)
-      logical,intent(in) :: debut ! signals first call to physics
-      logical,intent(in) :: lafin ! signals last call to physics
+      logical,intent(in) :: debut ! signals first CALL to physics
+      logical,intent(in) :: lafin ! signals last CALL to physics
       real,intent(in) :: pdtphys ! physics time step (s)
       real,intent(in) :: paprs(klon,klev+1) ! interlayer pressure (Pa)
@@ -75 +75 @@
 
 ! initializations
-if (debut) then ! Things to do only for the first call to physics 
+if (debut) then ! Things to do only for the first CALL to physics
 ! load initial conditions for physics (including the grid)
-  call phys_state_var_init() ! some initializations, required before calling phyetat0
-  call phyetat0("startphy.nc")
+  CALL phys_state_var_init() ! some initializations, required before calling phyetat0
+  CALL phyetat0("startphy.nc")
 
 ! Initialize outputs:
@@ -84 +84 @@
 !$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)
+  ! NB: getin() is not threadsafe; only one thread should CALL it.
+  CALL getin("iwrite_phys",iwrite_phys_omp)
 !$OMP END MASTER
 !$OMP BARRIER
@@ -93 +93 @@
   ! 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)
+  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)
+  CALL histbeg_phy("histins.nc",itau0,zjulian,dtime,nhori,nid_hist)
 #endif
 
@@ -105 +105 @@
 ! IOIPSL
   ! define vertical coordinate
-  call histvert(nid_hist,"presnivs","Vertical levels","Pa",klev, &
+  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', &
+  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', &
+  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', &
+  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', &
+  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)
+  CALL histend(nid_hist)
 #endif
 
@@ -166 +166 @@
 #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))
+  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
@@ -192 +192 @@
 ! if lastcall, then it is time to write "restartphy.nc" file
 if (lafin) then
-  call phyredem("restartphy.nc")
+  CALL phyredem("restartphy.nc")
 endif