Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/inlandsis.f90

Timestamp:

Jul 24, 2024, 1:17:08 PM (2 months ago)

Author:

abarral

Message:

Rename modules in misc from *_mod > lmdz_*
Put cbrt.f90, ch*.f90, pch*.f90 in new lmdz_libmath_pch.f90

File:

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

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

@@ -336 +336 @@
 
   data      Z0mBS0 /   0.5e-6/            ! MINimum Snow Roughness Length
-                                          ! ! for Momentum if Blowing Snow
-                                          ! ! Gallee et al. 2001 BLM 99 (19)
+                                          ! for Momentum if Blowing Snow
+                                          ! Gallee et al. 2001 BLM 99 (19)
   data      Z0m_S0/    0.00005/           ! MINimum Snow Roughness Length
-                                          ! ! MegaDunes    included
+                                          ! MegaDunes    included
   data      Z0m_S1/    0.030  /           ! MAXimum Snow Roughness Length
-                                          ! !        (Sastrugis)
+                                          !        (Sastrugis)
   ! #GL data      Z0_GIM/    0.0013/            ! Ice Min Z0 = 0.0013 m (Broeke)
-                                          ! ! Old Ice Z0 = 0.0500 m (Bruce)
-                                          ! !              0.0500 m (Smeets)
-                                          ! !              0.1200 m (Broeke)
+                                          ! Old Ice Z0 = 0.0500 m (Bruce)
+                                          !              0.0500 m (Smeets)
+                                          !              0.1200 m (Broeke)
   data      Z0_ICE/    0.0010/            ! Sea-Ice Z0 = 0.0010 m (Andreas)
-                                          ! !    (Ice Station Weddel -- ISW)
+                                          !    (Ice Station Weddel -- ISW)
   ! for aerolian erosion
   data      SblPom/ 1.27/   ! Lower Boundary Height Parameter
@@ -436 +436 @@
 
     rCd10n =  1. / 26.5 ! Vt / u*t = 26.5
-                 ! ! Budd et al. 1965, Antarct. Res. Series Fig.13
-                 ! ! ratio developped during assumed neutral conditions
+                 ! Budd et al. 1965, Antarct. Res. Series Fig.13
+                 ! ratio developped during assumed neutral conditions
 
 
@@ -554 +554 @@
         ! sss__G      =  0.27417     * gravit
 
-  ! !  ______________               _____
-  ! !  Newton-Raphson (! Iteration, BEGIN)
-  ! !  ~~~~~~~~~~~~~~               ~~~~~
+  !  ______________               _____
+  !  Newton-Raphson (! Iteration, BEGIN)
+  !  ~~~~~~~~~~~~~~               ~~~~~
         ! DO iit=1,nit
         ! sss__K      =  gravit      * r_Turb * A_Turb *za__SV(ikl)
@@ -580 +580 @@
    !      us__SV(ikl)= max(us__SV(ikl),epsi  )
    !      rCDmSV(ikl)=     us__SV(ikl)/VVa_OK
-  ! ! #AE     sss__F     =     vonkar     /rCDmSV(ikl)
+  ! #AE     sss__F     =     vonkar     /rCDmSV(ikl)
    !      ENDDO
 
-  ! !  ______________               ___
-  ! !  Newton-Raphson (! Iteration, END  )
-  ! !  ~~~~~~~~~~~~~~               ~~~
+  !  ______________               ___
+  !  Newton-Raphson (! Iteration, END  )
+  !  ~~~~~~~~~~~~~~               ~~~
 
    !      us_127      =  exp(    SblPom *log(us__SV(ikl)))
    !      us_227      =  us_127         *    us__SV(ikl)
 
-  ! !  Momentum            Turbulent Scale  u*: 0-Limit in case of no Blow. Snow
-  ! !  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+  !  Momentum            Turbulent Scale  u*: 0-Limit in case of no Blow. Snow
+  !  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    !      dusuth      =  us__SV(ikl) - usthSV(ikl)       ! u* - uth*
    !      signus      =  max(sign(unun,dusuth),zero)     ! 1 <=> u* - uth* > 0
@@ -664 +664 @@
         if (is_ok_density_kotlyakov) then
             tt_tmp = TaT_SV(ikl)-TfSnow
-            ! !vv_tmp = VV10SV(ikl)
+            !vv_tmp = VV10SV(ikl)
             vv_tmp=VV__SV(ikl) ! Etienne: use wind speed at first model level instead of 10m wind
   ! + ...         [ A compromise between
@@ -1004 +1004 @@
 
         e_pRad = 2.5   *  coszSV(ikl)       ! exponential argument,
-                                            ! ! V/nIR radiation partitioning,
-                                            ! ! DR97, 2, eqn (2.53) & (2.54)
+                                            ! V/nIR radiation partitioning,
+                                            ! DR97, 2, eqn (2.53) & (2.54)
         e1pRad = 1.-exp(-e_pRad)            ! exponential, V/nIR Rad. Part.
         exdRad= 1.
@@ -1186 +1186 @@
 
       TBr_sv(ikl) =sqrt(sqrt(IRu_SV(ikl)/StefBo))      ! Brightness
-                                                       ! ! Temperature
+                                                       ! Temperature
       uts_SV(ikl) =     (HSv_sv(ikl) +HSs_sv(ikl)) & ! u*T*
             /(rhT_SV(ikl) *cp)          !
@@ -1318 +1318 @@
   ! +  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
       A_Fact      =  1.0000        ! Andreas et al., 2004, p.4
-                                   ! ! ams.confex.com/ams/pdfpapers/68601.pdf
+                                   ! ams.confex.com/ams/pdfpapers/68601.pdf
 
   ! Parameterization of z0 dependance on Temperature (C. Amory, 2017)
@@ -1341 +1341 @@
         Z0_obs = exp(coefa*TaT_SV(ikl) + coefb)
       else if (TaT_SV(ikl)>=ta2 .and. TaT_SV(ikl)<ta3) then
-        ! ! if st > 0, melting induce smooth surface
+        ! if st > 0, melting induce smooth surface
         Z0_obs = exp(coefc*TaT_SV(ikl) + coefd)
       else