Changeset 3111
- Timestamp:
- Nov 2, 2023, 6:24:46 PM (14 months ago)
- Location:
- trunk/LMDZ.MARS/libf/phymars
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/conf_phys.F
r3098 r3111 21 21 use surfdat_h, only: albedo_h2o_cap,albedo_h2o_frost, 22 22 & frost_albedo_threshold, inert_h2o_ice, 23 & frost_metam_threshold 23 & frost_metam_threshold,old_wsublimation_scheme 24 24 use time_phylmdz_mod, only: ecritphy,day_step,iphysiq,ecritstart, 25 25 & daysec,dtphys … … 890 890 & frost_albedo_threshold 891 891 892 ! TMP: old_wsublimation_scheme 893 write(*,*) "Old water sublimation scheme?" 894 old_wsublimation_scheme = .true. 895 call getin_p("old_wsublimation_scheme",old_wsublimation_scheme) 896 write(*,*) "old_wsublimation_scheme",old_wsublimation_scheme 897 892 898 ! call Titus crocus line -- DEFAULT IS NONE 893 899 write(*,*) "Titus crocus line ?" -
trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90
r3098 r3111 20 20 real, save, allocatable :: d_coef(:,:) ! Diffusion coeficent 21 21 LOGICAL,SAVE :: lag_layer ! does lag layer is present? 22 !$OMP THREADPRIVATE(h2o_ice_depth, lag_co2_ice,albedo_perenialco2)22 !$OMP THREADPRIVATE(h2o_ice_depth,d_coef,lag_co2_ice,albedo_perenialco2) 23 23 24 24 CONTAINS -
trunk/LMDZ.MARS/libf/phymars/surfdat_h.F90
r2999 r3111 31 31 real,save :: z0_default ! default (constant over planet) surface roughness (m) 32 32 33 LOGICAL, SAVE :: old_wsublimation_scheme ! TEMPORARY : TO USE THE OLD WATER SUBLIMATION SCHEME (i.e., using Cd instead of Ch), true by default 34 33 35 !$OMP THREADPRIVATE(albedo_h2o_cap,albedo_h2o_frost,inert_h2o_ice, & 34 36 !$OMP frost_albedo_threshold,frost_metam_threshold,TESice_Ncoef, & 35 37 !$OMP TESice_Scoef,iceradius,dtemisice, & 36 !$OMP zmea,zstd,zsig,zgam,zthe,hmons,summit,base,z0,z0_default ) 38 !$OMP zmea,zstd,zsig,zgam,zthe,hmons,summit,base,z0,z0_default, & 39 !$OMP old_wsublimation_scheme) 37 40 38 41 !! mountain top dust flows -
trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
r3107 r3111 20 20 & igcm_hdo_vap, igcm_hdo_ice, 21 21 & igcm_stormdust_mass, igcm_stormdust_number 22 use surfdat_h, only: watercaptag, frost_albedo_threshold, dryness 22 use surfdat_h, only: watercaptag, frost_albedo_threshold, dryness, 23 & old_wsublimation_scheme 23 24 USE comcstfi_h, ONLY: cpp, r, rcp, g, pi 24 25 use watersat_mod, only: watersat … … 952 953 c the subtimestep. 953 954 saved_h2o_vap(:)= zq(:,1,igcm_h2o_vap) 954 955 955 DO ig=1,ngrid 956 956 subtimestep = ptimestep/nsubtimestep(ig) … … 963 963 zb(1:ngrid,2:nlay)=zkh(1:ngrid,2:nlay)*zb0(1:ngrid,2:nlay) 964 964 & /float(nsubtimestep(ig)) 965 zb(1:ngrid,1)=zcdv(1:ngrid)*zb0(1:ngrid,1) 965 if(old_wsublimation_scheme) then 966 zb(1:ngrid,1)=zcdv(1:ngrid)*zb0(1:ngrid,1) 966 967 & /float(nsubtimestep(ig)) 968 else 969 zb(1:ngrid,1)=zcdh(1:ngrid)*zb0(1:ngrid,1) 970 & /float(nsubtimestep(ig)) 971 endif 967 972 zb(1:ngrid,1)=dryness(1:ngrid)*zb(1:ngrid,1) 968 973 … … 986 991 zd(ig,1)=zb(ig,1)*z1(ig) 987 992 zq1temp(ig)=zc(ig,1)+ zd(ig,1)*qsat(ig) 988 989 zdqsdif(ig)=rho(ig)*dryness(ig)*zcdv(ig) 990 & *(zq1temp(ig)-qsat(ig)) 993 if(old_wsublimation_scheme) then 994 zdqsdif(ig)=rho(ig)*dryness(ig)*zcdv(ig) 995 & *(zq1temp(ig)-qsat(ig)) 996 else 997 zdqsdif(ig)=rho(ig)*dryness(ig)*zcdh(ig) 998 & *(zq1temp(ig)-qsat(ig)) 999 endif 991 1000 c write(*,*)'subliming more than available frost: qsurf!' 992 1001
Note: See TracChangeset
for help on using the changeset viewer.