Ignore:
Timestamp:
Nov 19, 2025, 12:40:01 PM (6 weeks ago)
Author:
emillour
Message:

Mars PCM:
Minor cleaning: don't systematically read/write "totcloudfrac" (water ice
cloud fraction) in startfi/restartfi files. Do it only if CLFvarying is .true.
EM

Location:
trunk/LMDZ.MARS/libf/phymars
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

    r3903 r3964  
    3838use tabfi_mod, only: tabfi
    3939use callkeys_mod, only: startphy_file, rdstorm, hdo
     40use callkeys_mod, only: CLFvarying, CLFfixval
    4041
    4142implicit none
     
    483484            minval(dtau), maxval(dtau)
    484485
    485 ! Sub-grid cloud fraction
    486 if (startphy_file) then
     486! Sub-grid water ice cloud fraction
     487totcloudfrac(:)=1.0 ! default value; no cloud fraction
     488if (startphy_file) then
     489  if (CLFvarying) then
    487490   call get_field("totcloudfrac",totcloudfrac,found,indextime)
    488491   if (.not.found) then
    489      write(*,*) "phyetat0: <totcloudfrac> not in file WARNING put to 1"
    490      totcloudfrac(:) = 1.0 !valeur par defaut (CLFfixval par defaut)
    491    endif
    492 else
    493    totcloudfrac(:)=1.0
     492     write(*,*) "phyetat0: <totcloudfrac> not in file; set to ",CLFfixval
     493   endif
     494   totcloudfrac(:)=CLFfixval
     495  endif ! of if (CLFvarying)
    494496endif ! if (startphy_file)
    495497write(*,*) "phyetat0: total cloud fraction <totcloudfrac> range:", &
  • trunk/LMDZ.MARS/libf/phymars/phyredem.F90

    r3808 r3964  
    207207  use paleoclimate_mod,    only: paleoclimate
    208208  use callkeys_mod,        only: calltherm, dustinjection, calllott_nonoro
     209  use callkeys_mod, only: CLFvarying
    209210  use paleoclimate_mod,    only: paleoclimate, h2o_ice_depth, lag_co2_ice, d_coef, zdqsdif_ssi_tot
    210211
     
    274275
    275276  ! Sub-grid cloud fraction
    276   call put_field("totcloudfrac","Total cloud fraction",totcloudfrac,time)
     277  if (CLFvarying) then
     278    call put_field("totcloudfrac","Total cloud fraction",totcloudfrac,time)
     279  endif
    277280 
    278281  ! Dust conversion factor
Note: See TracChangeset for help on using the changeset viewer.