Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/fonte_neige_mod.F90

Timestamp:

Jul 22, 2024, 9:29:09 PM (4 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/fonte_neige_mod.F90 (modified) (22 diffs)

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

@@ -1 @@
-!
+
 ! $Header$
-!
+
 MODULE fonte_neige_mod
-!
+
 ! This module will treat the process of snow, melting, accumulating, calving, in 
 ! case of simplified soil model.
-!
+
 !****************************************************************************************
   USE dimphy, ONLY : klon
@@ -52 +52 @@
 
 CONTAINS
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE fonte_neige_init(restart_runoff)
 
@@ -60 +60 @@
 ! The variable run_off_lic_0 is initialized to the field read from
 ! restart file. The other variables are initialized to zero.
-!
+
 !****************************************************************************************
 ! Input argument
@@ -73 +73 @@
 !****************************************************************************************
 ! Allocate run-off at landice and initilize with field read from restart
-!
+
 !****************************************************************************************
 
@@ -85 +85 @@
 !****************************************************************************************
 ! Allocate other variables and initilize to zero
-!
+
 !****************************************************************************************
     ALLOCATE(run_off_ter(klon), stat = error)
@@ -131 +131 @@
 !****************************************************************************************
 ! Read tau_calv
-!
+
 !****************************************************************************************
     CALL conf_interface(tau_calv)
@@ -150 +150 @@
     USE isotopes_verif_mod
 #endif
-!
+
 !****************************************************************************************
 ! Input argument
@@ -164 +164 @@
 !****************************************************************************************
 ! Allocate run-off at landice and initilize with field read from restart
-!
+
 !****************************************************************************************
 
@@ -185 +185 @@
 !****************************************************************************************
 ! Allocate other variables and initilize to zero
-!
+
 !****************************************************************************************
 
@@ -226 +226 @@
 #endif
 
-!
-!****************************************************************************************
-!
+!****************************************************************************************
+
   SUBROUTINE fonte_neige( knon, nisurf, knindex, dtime, &
        tsurf, precip_rain, precip_snow, &
@@ -258 +257 @@
 !   precip_rain  precipitations liquides
 !   precip_snow  precipitations solides
-!
+
 ! input/output:
 !   snow         champs hauteur de neige
@@ -264 +263 @@
 !   tsurf_new    temperature au sol
 !   evap
-!
+
   INCLUDE "YOETHF.h"
   INCLUDE "YOMCST.h"
@@ -329 +328 @@
 ! Start calculation
 ! - Initialization
-!
+
 !****************************************************************************************
     coeff_rel = dtime/(tau_calv * rday)
@@ -338 +337 @@
 ! - Increment snow due to precipitation and evaporation
 ! - Calculate the water balance due to precipitation and evaporation (bil_eau_s)
-!
+
 !****************************************************************************************
     WHERE (precip_snow > 0.) 
@@ -374 +373 @@
 ! - Calculate calving and decrement snow, if there are to much snow
 ! - Update temperature at surface
-!
+
 !****************************************************************************************
 
@@ -451 +450 @@
 ! Save ffonte, fqfonte and fqcalving in global arrays for each 
 ! sub-surface separately
-!
+
 !****************************************************************************************
     DO i = 1, knon
@@ -466 +465 @@
 
   END SUBROUTINE fonte_neige
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE fonte_neige_final(restart_runoff &
 #ifdef ISO      
@@ -474 +473 @@
 #endif   
                           )
-!
+
 ! This subroutine returns run_off_lic_0 for later writing to restart file.
-!
+
 #ifdef ISO
     USE infotrac_phy, ONLY: niso
@@ -484 +483 @@
 #endif
 #endif
-!
+
 !****************************************************************************************
     REAL, DIMENSION(klon), INTENT(OUT) :: restart_runoff
@@ -540 +539 @@
 
   END SUBROUTINE fonte_neige_final
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE fonte_neige_get_vars(pctsrf, fqcalving_out, &
               fqfonte_out, ffonte_out, run_off_lic_out &
@@ -553 +552 @@
 ! Cumulate ffonte, fqfonte and fqcalving respectively for
 ! all type of surfaces according to their fraction.
-!
+
 ! This routine is called from physiq.F before histwrite.
 !****************************************************************************************
@@ -608 +607 @@
 
   END SUBROUTINE fonte_neige_get_vars
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
 !#ifdef ISO
 !  subroutine fonte_neige_export_xtrun_off_lic_0(knon,xtrun_off_lic_0_diag)
 !    use infotrac_phy, ONLY: niso
-!
+
 !    ! inputs
 !    INTEGER, INTENT(IN)                      :: knon
 !    real, INTENT(IN), DIMENSION(niso,klon)   :: xtrun_off_lic_0_diag
-!
+
 !    xtrun_off_lic_0(:,:)=xtrun_off_lic_0_diag(:,:)
-!
+
 !  end subroutine fonte_neige_export_xtrun_off_lic_0
 !#endif