Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.f90

Timestamp:

Jul 24, 2024, 2:54:37 PM (2 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.f90 (modified) (5 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.f90

@@ -94 +94 @@
   ! +  ==================
 
-  integer :: ikl   ,isl   ,jsl   ,ist      !
-  integer :: ist__s,ist__w                 ! Soil/Water  Body Identifier
-  integer :: islsgn                        ! Soil/Snow Surfac.Identifier
-  real :: eps__3                        ! Arbitrary    Low Number
-  real :: etaMid,psiMid                 ! Layer Interface's Humidity
-  real :: mu_eta                        !     Soil thermal Conductivity
-  real :: mu_exp                        ! arg Soil thermal Conductivity
-  real :: mu_min                        ! Min Soil thermal Conductivity
-  real :: mu_max                        ! Max Soil thermal Conductivity
-  real :: mu_sno(knonv),mu_aux          !     Snow thermal Conductivity
-  real :: mu__dz(knonv,-nsol:nsno+1)    ! mu_(eta,sno)   / dz
-  real :: dtC_sv(knonv,-nsol:nsno)      ! dt      / C
-  real :: IRs__D(knonv)                 ! UpwardIR Previous Iter.Contr.
-  real :: dIRsdT(knonv)                 ! UpwardIR           T Derivat.
-  real :: f_HSHL(knonv)                 ! Factor common to HS and HL
-  real :: dRidTs(knonv)                 ! d(Rib)/d(Ts)
-  real :: HS___D(knonv)                 ! Sensible Heat Flux Atm.Contr.
-  real :: f___HL(knonv)                 !
-  real :: HL___D(knonv)                 ! Latent   Heat Flux Atm.Contr.
+  INTEGER :: ikl   ,isl   ,jsl   ,ist      !
+  INTEGER :: ist__s,ist__w                 ! Soil/Water  Body Identifier
+  INTEGER :: islsgn                        ! Soil/Snow Surfac.Identifier
+  REAL :: eps__3                        ! Arbitrary    Low Number
+  REAL :: etaMid,psiMid                 ! Layer Interface's Humidity
+  REAL :: mu_eta                        !     Soil thermal Conductivity
+  REAL :: mu_exp                        ! arg Soil thermal Conductivity
+  REAL :: mu_min                        ! Min Soil thermal Conductivity
+  REAL :: mu_max                        ! Max Soil thermal Conductivity
+  REAL :: mu_sno(knonv),mu_aux          !     Snow thermal Conductivity
+  REAL :: mu__dz(knonv,-nsol:nsno+1)    ! mu_(eta,sno)   / dz
+  REAL :: dtC_sv(knonv,-nsol:nsno)      ! dt      / C
+  REAL :: IRs__D(knonv)                 ! UpwardIR Previous Iter.Contr.
+  REAL :: dIRsdT(knonv)                 ! UpwardIR           T Derivat.
+  REAL :: f_HSHL(knonv)                 ! Factor common to HS and HL
+  REAL :: dRidTs(knonv)                 ! d(Rib)/d(Ts)
+  REAL :: HS___D(knonv)                 ! Sensible Heat Flux Atm.Contr.
+  REAL :: f___HL(knonv)                 !
+  REAL :: HL___D(knonv)                 ! Latent   Heat Flux Atm.Contr.
   REAL :: TSurf0(knonv),dTSurf          ! Previous Surface Temperature
-  real :: qsatsg(knonv) !,den_qs,arg_qs ! Soil   Saturat. Spec. Humidity
-  real :: dqs_dT(knonv)                 ! d(qsatsg)/dTv
-  real :: Psi(   knonv)                 ! 1st Soil Layer Water Potential
-  real :: RHuSol(knonv)                 ! Soil Surface Relative Humidity
-  real :: etaSol                        ! Soil Surface          Humidity
-  real :: d__eta                        ! Soil Surface Humidity Increm.
-  real :: Elem_A,Elem_C                 !   Diagonal Coefficients
-  real :: Diag_A(knonv,-nsol:nsno)      ! A Diagonal
-  real :: Diag_B(knonv,-nsol:nsno)      ! B Diagonal
-  real :: Diag_C(knonv,-nsol:nsno)      ! C Diagonal
-  real :: Term_D(knonv,-nsol:nsno)      !   Independant Term
-  real :: Aux__P(knonv,-nsol:nsno)      ! P Auxiliary Variable
-  real :: Aux__Q(knonv,-nsol:nsno)      ! Q Auxiliary Variable
-  real :: Ts_Min,Ts_Max                 ! Temperature Limits
+  REAL :: qsatsg(knonv) !,den_qs,arg_qs ! Soil   Saturat. Spec. Humidity
+  REAL :: dqs_dT(knonv)                 ! d(qsatsg)/dTv
+  REAL :: Psi(   knonv)                 ! 1st Soil Layer Water Potential
+  REAL :: RHuSol(knonv)                 ! Soil Surface Relative Humidity
+  REAL :: etaSol                        ! Soil Surface          Humidity
+  REAL :: d__eta                        ! Soil Surface Humidity Increm.
+  REAL :: Elem_A,Elem_C                 !   Diagonal Coefficients
+  REAL :: Diag_A(knonv,-nsol:nsno)      ! A Diagonal
+  REAL :: Diag_B(knonv,-nsol:nsno)      ! B Diagonal
+  REAL :: Diag_C(knonv,-nsol:nsno)      ! C Diagonal
+  REAL :: Term_D(knonv,-nsol:nsno)      !   Independant Term
+  REAL :: Aux__P(knonv,-nsol:nsno)      ! P Auxiliary Variable
+  REAL :: Aux__Q(knonv,-nsol:nsno)      ! Q Auxiliary Variable
+  REAL :: Ts_Min,Ts_Max                 ! Temperature Limits
   ! #e1 real      Exist0                        ! Existing Layer Switch
-  real :: psat_wat, psat_ice, sp        ! computation of qsat
-
-  integer :: nt_srf,it_srf,itEuBk          ! HL: Surface Scheme
+  REAL :: psat_wat, psat_ice, sp        ! computation of qsat
+
+  INTEGER :: nt_srf,it_srf,itEuBk          ! HL: Surface Scheme
   parameter(nt_srf=10)                     ! 10 before
-  real :: agpsrf,xgpsrf,dt_srf,dt_ver   !
-  real :: etaBAK(knonv)                 !
-  real :: etaNEW(knonv)                 !
-  real :: etEuBk(knonv)                 !
-  real :: fac_dt(knonv),faceta(knonv)   !
-  real :: PsiArg(knonv),SHuSol(knonv)   !
+  REAL :: agpsrf,xgpsrf,dt_srf,dt_ver   !
+  REAL :: etaBAK(knonv)                 !
+  REAL :: etaNEW(knonv)                 !
+  REAL :: etEuBk(knonv)                 !
+  REAL :: fac_dt(knonv),faceta(knonv)   !
+  REAL :: PsiArg(knonv),SHuSol(knonv)   !
 
 
@@ -436 +436 @@
             * exp (6827.*(1.         /273.16-1./TsisSV(ikl,isl)))
 
-      if(TsisSV(ikl,isl)<=273.16) then
+      IF(TsisSV(ikl,isl)<=273.16) THEN
         qsatsg(ikl) = 0.622 * psat_ice / (sp - 0.378 * psat_ice)
       else
@@ -466 +466 @@
         END DO
 
-    if(ist__s==1) then ! to reduce computer time
+    IF(ist__s==1) then ! to reduce computer time
 
     DO it_srf=1,nt_srf                                     !
@@ -615 +615 @@
       TsisSV(ikl,isl)   = Aux__Q(ikl,isl)  *TsisSV(ikl,isl+1) &
             +TsisSV(ikl,isl)
-      if(isl==0.and.isnoSV(ikl)==0) then
-
+      IF(isl==0.and.isnoSV(ikl)==0) THEN
        TsisSV(ikl,isl)  = min(TaT_SV(ikl)+30,TsisSV(ikl,isl))
        TsisSV(ikl,isl)  = max(TaT_SV(ikl)-30,TsisSV(ikl,isl))
@@ -673 +672 @@
 
 
-end subroutine sisvat_tso
+END SUBROUTINE sisvat_tso