Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmd/surf_landice_mod.F90

Timestamp:

Jul 24, 2024, 4:23:34 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/surf_landice_mod.F90 (modified) (9 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/surf_landice_mod.F90

@@ -66 +66 @@
     !****************************************************************************************
     INTEGER, INTENT(IN) :: itime, knon
-    INTEGER, DIMENSION(klon), INTENT(in) :: knindex
-    REAL, INTENT(in) :: dtime
+    INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
+    REAL, INTENT(IN) :: dtime
     REAL, DIMENSION(klon), INTENT(IN) :: swnet ! net shortwave radiance
     REAL, DIMENSION(klon), INTENT(IN) :: lwnet ! net longwave radiance
@@ -135 +135 @@
 #ifdef ISO
     REAL, DIMENSION(ntiso,klon), INTENT(OUT)     :: xtevap     
-!    real, DIMENSION(niso,klon) :: xtrun_off_lic_0_diag ! est une variable globale de
+!    REAL, DIMENSION(niso,klon) :: xtrun_off_lic_0_diag ! est une variable globale de
 !    fonte_neige
 #endif
@@ -163 +163 @@
     REAL, DIMENSION(niso,klon) :: xtsnow_prec,xtsol_prec
     REAL, DIMENSION(klon) :: snow_prec,qsol_prec
-!    real, DIMENSION(klon) :: run_off_lic_0_diag
+!    REAL, DIMENSION(klon) :: run_off_lic_0_diag
 #endif
 
@@ -433 +433 @@
       !****************************************************************************************
 
-      if (z0m_landice > 0.) THEN
+      IF (z0m_landice > 0.) THEN
         z0m(1:knon) = z0m_landice
         z0h(1:knon) = z0h_landice
@@ -449 +449 @@
         coefd = log(z03) - coefc * ta3
         do j = 1, knon
-          if (temp_air(j) < ta1) THEN
+          IF (temp_air(j) < ta1) THEN
             z0m(j) = z01
-          else if (temp_air(j)>=ta1 .and. temp_air(j)<ta2) THEN
+          ELSE IF (temp_air(j)>=ta1 .AND. temp_air(j)<ta2) THEN
             z0m(j) = exp(coefa * temp_air(j) + coefb)
-          else if (temp_air(j)>=ta2 .and. temp_air(j)<ta3) THEN
+          ELSE IF (temp_air(j)>=ta2 .AND. temp_air(j)<ta3) THEN
             ! if st > 0, melting induce smooth surface
             z0m(j) = exp(coefc * temp_air(j) + coefd)
@@ -473 +473 @@
       ! this layer
 
-      if (ok_bs) THEN
+      IF (ok_bs) THEN
         fluxbs(:) = 0.
         do j = 1, knon
@@ -484 +484 @@
         ! 1st step: erosion of fresh snow accumulated during the time step
         do j = 1, knon
-          if (precip_snow(j) > 0.) THEN
+          IF (precip_snow(j) > 0.) THEN
             rhos(j) = rhofresh_bs
             ! blowing snow flux formula used in MAR
@@ -507 +507 @@
             fluxbs_1(j) = max(-maxerosion, fluxbs_1(j))
 
-            if (precip_snow(j) > abs(fluxbs_1(j))) THEN
+            IF (precip_snow(j) > abs(fluxbs_1(j))) THEN
               ok_remaining_freshsnow(j) = .TRUE.
               bsweight_fresh(j) = 1.
@@ -530 +530 @@
         ! which depends on surface snow density
         do j = 1, knon
-          if (ok_remaining_freshsnow(j)) THEN
+          IF (ok_remaining_freshsnow(j)) THEN
             fluxbs_2(j) = 0.
           else