Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer

Timestamp:

Jul 23, 2024, 3:29:36 PM (4 months 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]

Location:

LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer

Files:

• nucleation_tstep_mod.F90 (modified) (3 diffs)
• stratH2O_methox.F90 (modified) (1 diff)
• sulfate_aer_mod.F90 (modified) (4 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/nucleation_tstep_mod.F90

@@ -135 +135 @@
 
 
-  !    Fortran 90 subroutine binapara
+  !    Fortran 90 SUBROUTINE binapara
 
   !    Calculates parametrized values of nucleation rate,
@@ -299 +299 @@
      kinetic_n,kinetic_i)
 
-  !    Fortran 90 subroutine newbinapara
+  !    Fortran 90 SUBROUTINE newbinapara
 
   !    Calculates parametrized values for neutral and ion-induced sulfuric acid-water particle formation rate
@@ -614 +614 @@
      kinrhotresi=EXP(kinrhotresi) !1/cm3
      
-     IF (kinrhotresi<rhoali) kinetic_i=.true.
+     IF (kinrhotresi<rhoali) kinetic_i=.TRUE.
      
      IF (kinetic_i) THEN    

LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratH2O_methox.F90

@@ -30 +30 @@
   LOGICAL, INTENT(IN) :: debutphy   ! flag for first physiq step
 ! Output variables
-! tendency buffer used in add_phys_tend subroutine (in physiq_mod)
+! tendency buffer used in add_phys_tend SUBROUTINE (in physiq_mod)
   REAL, INTENT(INOUT), DIMENSION(klon,klev)  :: dq_ch4mmr
   

LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90

@@ -72 +72 @@
 !        include conversion from molec/cm3 to Pa
 !        ph2so4=solh2so4(temp,xa)*(1.38065e-16*temp)/10.
-!        print*,' ph2so4=',ph2so4,temp,R2SO4(ilon,ilev)
+!        PRINT*,' ph2so4=',ph2so4,temp,R2SO4(ilon,ilev)
 !        good match with Martin, et Ayers, not with Gmitro (the famous 0.086)
 
@@ -135 +135 @@
       (R2SO4B(ilon,ilev,IK)*1.e-2))**third
 
-!             print*,'R,Rwet(m),kelvin,h2so4(%),ro=',RRSI(ik),radwet,fkelvin, &
+!             PRINT*,'R,Rwet(m),kelvin,h2so4(%),ro=',RRSI(ik),radwet,fkelvin, &
 !              &  R2SO4B(ilon,ilev,IK),DENSO4B(ilon,ilev,IK)
-!             print*,' equil.h2so4(molec/cm3), & 
+!             PRINT*,' equil.h2so4(molec/cm3), &
 !              & sigma',solh2so4(temp,xa),surftension(temp,xa)
 
@@ -830 +830 @@
          psh2o_out=100.*exp(psh2o_out)
       end if
-!      print*,psh2o_out
+!      PRINT*,psh2o_out
       
        return
@@ -993 +993 @@
 !        h2so4 weight precent
          wph2so4_out = 9800.*sulfmolal/(98.*sulfmolal+1000.)
-!         print*,rh,pph2o,psh2o(T),vpice(T)
-!         print*,T,aw,sulfmolal,wph2so4_out
+!         PRINT*,rh,pph2o,psh2o(T),vpice(T)
+!         PRINT*,T,aw,sulfmolal,wph2so4_out
          wph2so4_out = max(wph2so4_out,15.)
          wph2so4_out = min(wph2so4_out,99.999)