Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/phydev/physiq_mod.F90

Timestamp:

Jul 23, 2024, 3:29:36 PM (8 weeks ago)

Author:

abarral

Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:

• LMDZ6/branches/Amaury_dev/libf/phydev/physiq_mod.F90 (modified) (5 diffs)

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

@@ -56 +56 @@
 real :: temp_newton(klon,klev)
 integer :: k
-logical, save :: first=.true.
+logical, save :: first=.TRUE.
 !$OMP THREADPRIVATE(first)
 
@@ -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
@@ -137 +137 @@
 !$OMP END MASTER
 !$OMP BARRIER
-endif ! of if (debut)
+END IF ! of if (debut)
 
 ! increment local time counter itau
@@ -154 +154 @@
 d_v(1:klon,1)=-v(1:klon,1)/86400.
 ! newtonian relaxation towards temp_newton()
-do k=1,klev
+DO k=1,klev
   temp_newton(1:klon,k)=280.+cos(latitude(1:klon))*40.-pphi(1:klon,k)/rg*6.e-3
   d_t(1:klon,k)=(temp_newton(1:klon,k)-t(1:klon,k))/1.e5
-enddo
+END DO
 
 
-print*,'PHYDEV: itau=',itau
+PRINT*,'PHYDEV: itau=',itau
 
 ! 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
+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))
+END IF
 #endif
 
@@ -191 +191 @@
 
 ! if lastcall, then it is time to write "restartphy.nc" file
-if (lafin) then
+IF (lafin) then
   CALL phyredem("restartphy.nc")
-endif
+END IF
 
-end subroutine physiq
+END SUBROUTINE  physiq
 
 END MODULE physiq_mod