Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_ts2.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_ts2.f90 (modified) (8 diffs)

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

@@ -95 +95 @@
   ! +  ==================
 
-  integer :: ig, jk, isl
-  real :: mu
-  real :: Tsrf(klonv)               ! surface temperature as extrapolated from soil
-  real :: mug(klonv)                 !hj coef top layers
-  real :: ztherm_i(klonv), zdz2(klonv, -nsol:nsno), z1s
-  real :: pfluxgrd(klonv), pcapcal(klonv), cal(klonv)
-  real :: beta(klonv), dif_grnd(klonv)
-  real :: C_coef(klonv, -nsol:nsno), D_coef(klonv, -nsol:nsno)
+  INTEGER :: ig, jk, isl
+  REAL :: mu
+  REAL :: Tsrf(klonv)               ! surface temperature as extrapolated from soil
+  REAL :: mug(klonv)                 !hj coef top layers
+  REAL :: ztherm_i(klonv), zdz2(klonv, -nsol:nsno), z1s
+  REAL :: pfluxgrd(klonv), pcapcal(klonv), cal(klonv)
+  REAL :: beta(klonv), dif_grnd(klonv)
+  REAL :: C_coef(klonv, -nsol:nsno), D_coef(klonv, -nsol:nsno)
 
   REAL, DIMENSION(klonv) :: zx_mh, zx_nh, zx_oh
@@ -123 +123 @@
 
 
-  ! write(*,*)'T check'
+  ! WRITE(*,*)'T check'
 
   ! DO  ig = 1,knonv
@@ -196 +196 @@
   ! IF (knonv>0) THEN
   !  DO ig=1,8
-  !    write(*,*)ig,'sisvat: Tsis ',TsisSV(ig,isnoSV(ig))
-  !    write(*,*)'max-1            ',TsisSV(ig,isnoSV(ig)-1)
-  !    write(*,*)'max-2            ',TsisSV(ig,isnoSV(ig)-2)
-  !    write(*,*)'0                ',TsisSV(ig,0)
-  !!        write(*,*)min(max(isnoSV(ig),0),1),max(1-isnoSV(ig),0)
+  !    WRITE(*,*)ig,'sisvat: Tsis ',TsisSV(ig,isnoSV(ig))
+  !    WRITE(*,*)'max-1            ',TsisSV(ig,isnoSV(ig)-1)
+  !    WRITE(*,*)'max-2            ',TsisSV(ig,isnoSV(ig)-2)
+  !    WRITE(*,*)'0                ',TsisSV(ig,0)
+  !!        WRITE(*,*)min(max(isnoSV(ig),0),1),max(1-isnoSV(ig),0)
   !  ENDDO
   ! END IF
@@ -216 +216 @@
 
     IF (mug(ig)  <= 0.05) THEN
-      write(*, *)'Attention mu low', mug(ig)
+      WRITE(*, *)'Attention mu low', mug(ig)
     ENDIF
     IF (mug(ig)  >= 0.98) THEN
-      write(*, *)'Attention mu high', mug(ig)
+      WRITE(*, *)'Attention mu high', mug(ig)
     ENDIF
 
@@ -243 +243 @@
     ENDDO
   ENDDO
-  ! write(*,*)ig,'Tsis',TsisSV(ig,0)
+  ! WRITE(*,*)ig,'Tsis',TsisSV(ig,0)
 
   ! IF (indice == is_sic) THEN
@@ -329 +329 @@
   DO ig = 1, knonv
     IF (ps__SV(ig)<1.) THEN
-      ! write(*,*)'ig',ig,'ps',ps__SV(ig)
+      ! WRITE(*,*)'ig',ig,'ps',ps__SV(ig)
       ps__SV(ig) = max(ps__SV(ig), 1.e-8)
     ENDIF
     IF (p1l_SV(ig)<1.) THEN
-      ! write(*,*)'ig',ig,'p1l',p1l_SV(ig)
+      ! WRITE(*,*)'ig',ig,'p1l',p1l_SV(ig)
       p1l_SV(ig) = max(p1l_SV(ig), 1.e-8)
     ENDIF
     IF (TaT_SV(ig)<180.) THEN
-      ! write(*,*)'ig',ig,'TaT',TaT_SV(ig)
+      ! WRITE(*,*)'ig',ig,'TaT',TaT_SV(ig)
       TaT_SV(ig) = max(TaT_SV(ig), 180.)
     ENDIF
     IF (QaT_SV(ig)<1.e-8) THEN
-      ! write(*,*)'ig',ig,'QaT',QaT_SV(ig)
+      ! WRITE(*,*)'ig',ig,'QaT',QaT_SV(ig)
       QaT_SV(ig) = max(QaT_SV(ig), 1.e-8)
     ENDIF
     IF (Tsf_SV(ig)<100.) THEN
-      ! write(*,*)'ig',ig,'Tsf',Tsf_SV(ig)
+      ! WRITE(*,*)'ig',ig,'Tsf',Tsf_SV(ig)
       Tsf_SV(ig) = max(Tsf_SV(ig), 180.)
     ENDIF
     IF (Tsf_SV(ig)>500.) THEN
-      ! write(*,*)'ig',ig,'Tsf',Tsf_SV(ig)
+      ! WRITE(*,*)'ig',ig,'Tsf',Tsf_SV(ig)
       Tsf_SV(ig) = min(Tsf_SV(ig), 400.)
     ENDIF
     ! IF (Tsrf(ig).LT.1.) THEN
-    !!          write(*,*)'ig',ig,'Tsrf',Tsrf(ig)
+    !!          WRITE(*,*)'ig',ig,'Tsrf',Tsrf(ig)
     !   Tsrf(ig)=max(Tsrf(ig),TaT_SV(ig)-20.)
     ! ENDIF
     IF (cdH_SV(ig)<1.e-10) THEN
-      ! IF (ig.le.3)   write(*,*)'ig',ig,'cdH',cdH_SV(ig)
+      ! IF (ig.le.3)   WRITE(*,*)'ig',ig,'cdH',cdH_SV(ig)
       cdH_SV(ig) = .5
     ENDIF
@@ -370 +370 @@
       zx_qs = r2es * FOEEW(Tsf_SV(ig), zdelta) / ps__SV(ig)
       zx_qs = MIN(0.5, zx_qs)
-      !write(*,*)'zcor',retv*zx_qs
+      !WRITE(*,*)'zcor',retv*zx_qs
       zcor = 1. / (1. - retv * zx_qs)
       zx_qs = zx_qs * zcor
@@ -470 +470 @@
   ENDDO
 
-end subroutine sisvat_ts2
+END SUBROUTINE sisvat_ts2