Ignore:
Timestamp:
Jul 24, 2024, 1:17:08 PM (7 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

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
Files:
8 edited

Legend:

Unmodified
Added
Removed

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

    r5105 r5113  
    336336
    337337  data      Z0mBS0 /   0.5e-6/            ! MINimum Snow Roughness Length
    338                                           ! ! for Momentum if Blowing Snow
    339                                           ! ! Gallee et al. 2001 BLM 99 (19)
     338                                          ! for Momentum if Blowing Snow
     339                                          ! Gallee et al. 2001 BLM 99 (19)
    340340  data      Z0m_S0/    0.00005/           ! MINimum Snow Roughness Length
    341                                           ! ! MegaDunes    included
     341                                          ! MegaDunes    included
    342342  data      Z0m_S1/    0.030  /           ! MAXimum Snow Roughness Length
    343                                           ! !        (Sastrugis)
     343                                          !        (Sastrugis)
    344344  ! #GL data      Z0_GIM/    0.0013/            ! Ice Min Z0 = 0.0013 m (Broeke)
    345                                           ! ! Old Ice Z0 = 0.0500 m (Bruce)
    346                                           ! !              0.0500 m (Smeets)
    347                                           ! !              0.1200 m (Broeke)
     345                                          ! Old Ice Z0 = 0.0500 m (Bruce)
     346                                          !              0.0500 m (Smeets)
     347                                          !              0.1200 m (Broeke)
    348348  data      Z0_ICE/    0.0010/            ! Sea-Ice Z0 = 0.0010 m (Andreas)
    349                                           ! !    (Ice Station Weddel -- ISW)
     349                                          !    (Ice Station Weddel -- ISW)
    350350  ! for aerolian erosion
    351351  data      SblPom/ 1.27/   ! Lower Boundary Height Parameter
    … …  
    436436
    437437    rCd10n =  1. / 26.5 ! Vt / u*t = 26.5
    438                  ! ! Budd et al. 1965, Antarct. Res. Series Fig.13
    439                  ! ! ratio developped during assumed neutral conditions
     438                 ! Budd et al. 1965, Antarct. Res. Series Fig.13
     439                 ! ratio developped during assumed neutral conditions
    440440
    441441
    … …  
    554554        ! sss__G      =  0.27417     * gravit
    555555
    556   ! ! ______________               _____
    557   ! ! Newton-Raphson (! Iteration, BEGIN)
    558   ! ! ~~~~~~~~~~~~~~               ~~~~~
     556  ! ______________               _____
     557  ! Newton-Raphson (! Iteration, BEGIN)
     558  ! ~~~~~~~~~~~~~~               ~~~~~
    559559        ! DO iit=1,nit
    560560        ! sss__K      =  gravit      * r_Turb * A_Turb *za__SV(ikl)
    … …  
    580580   !      us__SV(ikl)= max(us__SV(ikl),epsi  )
    581581   !      rCDmSV(ikl)=     us__SV(ikl)/VVa_OK
    582   ! ! #AE     sss__F     =     vonkar     /rCDmSV(ikl)
     582  ! #AE     sss__F     =     vonkar     /rCDmSV(ikl)
    583583   !      ENDDO
    584584
    585   ! ! ______________               ___
    586   ! ! Newton-Raphson (! Iteration, END  )
    587   ! ! ~~~~~~~~~~~~~~               ~~~
     585  ! ______________               ___
     586  ! Newton-Raphson (! Iteration, END  )
     587  ! ~~~~~~~~~~~~~~               ~~~
    588588
    589589   !      us_127      =  exp(    SblPom *log(us__SV(ikl)))
    590590   !      us_227      =  us_127         *    us__SV(ikl)
    591591
    592   ! ! Momentum            Turbulent Scale  u*: 0-Limit in case of no Blow. Snow
    593   ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
     592  ! Momentum            Turbulent Scale  u*: 0-Limit in case of no Blow. Snow
     593  ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    594594   !      dusuth      =  us__SV(ikl) - usthSV(ikl)       ! u* - uth*
    595595   !      signus      =  max(sign(unun,dusuth),zero)     ! 1 <=> u* - uth* > 0
    … …  
    664664        if (is_ok_density_kotlyakov) then
    665665            tt_tmp = TaT_SV(ikl)-TfSnow
    666             ! !vv_tmp = VV10SV(ikl)
     666            !vv_tmp = VV10SV(ikl)
    667667            vv_tmp=VV__SV(ikl) ! Etienne: use wind speed at first model level instead of 10m wind
    668668  ! + ...         [ A compromise between
    … …  
    10041004
    10051005        e_pRad = 2.5   *  coszSV(ikl)       ! exponential argument,
    1006                                             ! ! V/nIR radiation partitioning,
    1007                                             ! ! DR97, 2, eqn (2.53) & (2.54)
     1006                                            ! V/nIR radiation partitioning,
     1007                                            ! DR97, 2, eqn (2.53) & (2.54)
    10081008        e1pRad = 1.-exp(-e_pRad)            ! exponential, V/nIR Rad. Part.
    10091009        exdRad= 1.
    … …  
    11861186
    11871187      TBr_sv(ikl) =sqrt(sqrt(IRu_SV(ikl)/StefBo))      ! Brightness
    1188                                                        ! ! Temperature
     1188                                                       ! Temperature
    11891189      uts_SV(ikl) =     (HSv_sv(ikl) +HSs_sv(ikl)) & ! u*T*
    11901190            /(rhT_SV(ikl) *cp)          !
    … …  
    13181318  ! +  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    13191319      A_Fact      =  1.0000        ! Andreas et al., 2004, p.4
    1320                                    ! ! ams.confex.com/ams/pdfpapers/68601.pdf
     1320                                   ! ams.confex.com/ams/pdfpapers/68601.pdf
    13211321
    13221322  ! Parameterization of z0 dependance on Temperature (C. Amory, 2017)
    … …  
    13411341        Z0_obs = exp(coefa*TaT_SV(ikl) + coefb)
    13421342      else if (TaT_SV(ikl)>=ta2 .and. TaT_SV(ikl)<ta3) then
    1343         ! ! if st > 0, melting induce smooth surface
     1343        ! if st > 0, melting induce smooth surface
    13441344        Z0_obs = exp(coefc*TaT_SV(ikl) + coefd)
    13451345      else

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

    r5105 r5113  
    5050  ! +  =======
    5151
    52   ! !DO isn = nsno,2,-1
     52  !DO isn = nsno,2,-1
    5353  DO ikl = 1,knonv
    5454
    … …  
    7070          /max(epsi,ro__SV(ikl,isn))
    7171
    72     ! ! Densification of the uppermost snow layer if erosion:
     72    ! Densification of the uppermost snow layer if erosion:
    7373    if((dzweqo-dzweqn)>0                    .and. &
    7474          dzsnSV(ikl,isn)>0                    .and. &
    7575          ro__SV(ikl,max(1,isnoSV(ikl)))<roBdSV) then
    7676
    77     ! !characteristic time scale for drifting snow compaction set to 24h
    78     ! !linear densification rate [kg/m3/s] over 24h
     77    !characteristic time scale for drifting snow compaction set to 24h
     78    !linear densification rate [kg/m3/s] over 24h
    7979    densif         = (450. - frsno) / (3600*24)
    8080
    81     ! !Attenuation of compaction rate from 450 to 500 kg/m3
     81    !Attenuation of compaction rate from 450 to 500 kg/m3
    8282    Fac         = 1-((ro__SV(ikl,max(1,isnoSV(ikl))) &
    8383          -roBdSV)/(500.-roBdSV))
    … …  
    103103
    104104  END DO
    105   ! !END DO
     105  !END DO
    106106
    107107

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

    r5105 r5113  
    7373  ! #vm real     EqSn01(knonv)                 ! Snow Energy, after Phase Change
    7474  ! #vm real     EqSn02(knonv)                 ! Snow Energy, after Phase Change
    75                                          ! !              .AND. Last Melting
     75                                         !              .AND. Last Melting
    7676  ! #e1 real     EqSn_1(knonv)                 ! Snow Energy, after Phase Change
    77                                          ! !              .AND. Mass Redistr.
     77                                         !              .AND. Mass Redistr.
    7878  ! Snow/Ice (Mass) Budget
    7979  ! ~~~~~~~~~~~~~~~~~~~~~~
    … …  
    160160  data      dz_Min/2.5e-3/                ! Minim. Snow Layer Thickness
    161161  data      SGDmax/0.003/                 ! Maxim. Snow Grain Diameter  [m]
    162                                           ! ! (Rowe et al. 1995, JGR p.16268)
     162                                          ! (Rowe et al. 1995, JGR p.16268)
    163163
    164164  ! +--Energy Budget (IN)
    … …  
    206206    noSnow(ikl)   = 0                   ! Nb of Layers Updater
    207207    ispiSV(ikl)   = 0                   ! Pore Hole Close OFF Index
    208                                         ! ! (assumed to be the Top of
    209                                         ! ! the surimposed Ice Layer)
     208                                        ! (assumed to be the Top of
     209                                        ! the surimposed Ice Layer)
    210210    zn5_SV(ikl)   = 0.
    211211    rusnSV0(ikl)  = 0.
    … …  
    493493
    494494      if(isnoSV(ikl)<=1 .OR. opt_runoff_ac) rusnew = 0.
    495       ! !if(ivgtSV(ikl)>=1) rusnew = 0.
     495      !if(ivgtSV(ikl)>=1) rusnew = 0.
    496496
    497497  ! #EU                        rusnew = 0.

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

    r5105 r5113  
    253253
    254254      dalbeW =(0.64 - csegal  )*0.0625  ! Warren 1982, RevGeo, fig.12b
    255                                         ! ! 0.0625 = 5% * 1/0.8,   p.81
    256                                         ! ! 0.64   = cos(50)
     255                                        ! 0.0625 = 5% * 1/0.8,   p.81
     256                                        ! 0.64   = cos(50)
    257257      dalbed =     dalbeW      *       min(1,isnoSV(ikl))
    258258  !-------------------------------------------------------------------------
    … …  
    276276      albSn6(4:6)=albSn2
    277277
    278       ! !snow albedo corection if wetsnow
     278      !snow albedo corection if wetsnow
    279279  ! #GL     albSn1 =  albSn1*max(0.9,(1.-1.5*eta_SV(ikl,isn)))
    280280  ! #GL     albSn2 =  albSn2*max(0.9,(1.-1.5*eta_SV(ikl,isn)))
    … …  
    410410  ! +--Integrated Snow/Ice/Soil Albedo: Clouds Correction! Greuell & all., 1994
    411411  ! +  --------------------------------------------------! Glob.&t Planet.Change
    412                                                    ! ! (9):91-114
     412                                                   ! (9):91-114
    413413        alb1sv(ikl) = alb1sv(ikl) + 0.05 *(cld_SV(ikl)-0.5)*SIcenH &
    414414              + dalbed      *    (1.-cld_SV(ikl))
    … …  
    460460    !AO_CK 20/02/2020
    461461
    462     ! ! No check if coupling update since MAR and NEMO albedo are too different
     462    ! No check if coupling update since MAR and NEMO albedo are too different
    463463    ! and since MAR albedo is computed on properties that are not in NEMO
    464     ! ! prescription for each time step with NEMO values
     464    ! prescription for each time step with NEMO values
    465465
    466466  ! #AO      if (LSmask(ikl) .eq. 0 .and. coupling_ao .eq. .TRUE.) then
    … …  
    511511      DO ikl=1,knonv
    512512        sEX_sv(ikl,isn) = 1.0
    513        ! !sEX_sv(ikl,isn) = 0.95 ! if MAR is too warm in summer
     513       !sEX_sv(ikl,isn) = 0.95 ! if MAR is too warm in summer
    514514      END DO
    515515    END DO

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

    r5105 r5113  
    265265    IRs__D(ig) = dIRsdT(ig)* TsisSV(ig, isl) * 0.75                  !:
    266266  END DO
    267   ! !hj
     267  !hj
    268268  !!-----------------------------------------------------------------------
    269269  !! 3)
    … …  
    370370      zx_qs = r2es * FOEEW(Tsf_SV(ig), zdelta) / ps__SV(ig)
    371371      zx_qs = MIN(0.5, zx_qs)
    372       ! !write(*,*)'zcor',retv*zx_qs
     372      !write(*,*)'zcor',retv*zx_qs
    373373      zcor = 1. / (1. - retv * zx_qs)
    374374      zx_qs = zx_qs * zcor

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

    r5105 r5113  
    372372  ! +  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
    373373
    374     ! !EV!isl=  min(isnoSV(1)+1,nsno)
     374    !EV!isl=  min(isnoSV(1)+1,nsno)
    375375
    376376    DO ikl=1,knonv
    … …  
    426426      !  sp = (pst_SV(ikl) + ptopSV) * 10.
    427427
    428       ! !sp=ps__SV(ikl)
    429       ! ! Etienne: in the formula herebelow sp should be in hPa, not
    430       ! ! in Pa so I divide by 100.
     428      !sp=ps__SV(ikl)
     429      ! Etienne: in the formula herebelow sp should be in hPa, not
     430      ! in Pa so I divide by 100.
    431431      sp=ps__SV(ikl)/100.
    432432      psat_ice = 6.1070 * exp(6150. *(1./273.16 - &
    … …  
    476476  !                                                              ! by Atm.Conten
    477477  ! #??.        *max(0,sign(1.,qsatsg(ikl)-QaT_SV(ikl))))        ! NO Limitation
    478                                                            ! ! of Downw.Flux
     478                                                           ! of Downw.Flux
    479479        END DO                                             !
    480480      DO itEuBk=1,2                                        !
    … …  
    522522      ist__w=     1-ist__s                                 ! 1 if    H2O
    523523      d__eta     =  eta_SV(ikl,isl)-etaNEW(ikl)            !
    524       ! ! latent heat flux computation
     524      ! latent heat flux computation
    525525      HL___D(ikl)=( ist__s *ro_Wat *dz_dSV(0) & ! Soil Contrib.
    526526            *(etaNEW(ikl)    -etaBAK(ikl)) / dt__SV & !
    … …  
    623623  ! #EU      TsisSV(ikl,isl)  = max(TaT_SV(ikl)-15.,TsisSV(ikl,isl))
    624624
    625       ! !XF 18/11/2018 to avoid ST reaching 70°C!!
    626       ! !It is an error compensation but does not work over tundra
     625      !XF 18/11/2018 to avoid ST reaching 70°C!!
     626      !It is an error compensation but does not work over tundra
    627627
    628628      endif

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

    r5105 r5113  
    140140  DO ikl=1,knonv
    141141  DO   isn=1,isnoSV(ikl)-3 ! no agregation of 3 first snowlayers
    142                            ! ! XF 04/07/2019
     142                           ! XF 04/07/2019
    143143
    144144      isno_n    =             isnoSV(ikl)-isn+1     ! Snow Normal.Profile
    … …  
    526526  DO ikl=1,knonv
    527527    DO isn=1,isnoSV(ikl)-3 ! no agregation of 3 first snowlayers
    528                            ! ! XF 04/07/2019
     528                           ! XF 04/07/2019
    529529
    530530      isno_n    =             isnoSV(ikl)-isn+1     ! Snow Normal.Profile

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90

    r5110 r5113  
    783783        min_period = 1800. ! en secondes
    784784        dalph_soil = 2.    ! rapport entre les epaisseurs de 2 couches succ.
    785         ! !$OMP MASTER
     785        !$OMP MASTER
    786786        !     IF (is_mpi_root) THEN
    787787        !        OPEN(99,file='soil.def',status='old',form='formatted',iostat=ierr)
    … …  
    795795        !        END IF
    796796        !     ENDIF
    797         ! !$OMP END MASTER
     797        !$OMP END MASTER
    798798        !     CALL bcast(min_period)
    799799        !     CALL bcast(dalph_soil)
Note: See TracChangeset for help on using the changeset viewer.