Changeset 2630
- Timestamp:
- Mar 1, 2022, 9:16:34 PM (3 years ago)
- Location:
- trunk/LMDZ.MARS
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/README
r2628 r2630 3625 3625 to the new terminology "mountain top dust flows", accordingly to ticket #71. 3626 3626 The new flag read in callphys.def is "topflows". 3627 3628 == 01/03/2022 == AB+LL 3629 Fix ticket #93. start2archive now correctly stores the variable albedo from 3630 startfi.nc in the output file start_archive.nc, and newstart reads the variable 3631 albedo in start_archive.nc before outputing it in restartfi.nc -
trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F
r2466 r2630 1 1 SUBROUTINE lect_start_archive(ngrid,nlayer,nqtot, 2 & date,tsurf,tsoil, emis,q2,2 & date,tsurf,tsoil,albedo,emis,q2, 3 3 & t,ucov,vcov,ps,co2ice,h,phisold_newgrid, 4 4 & q,qsurf,tauscaling,totcloudfrac,surfith,nid, … … 47 47 REAL,INTENT(OUT) :: tsoil(ngrid,nsoilmx) ! soil temperature 48 48 REAL,INTENT(OUT) :: co2ice(ngrid) ! CO2 ice layer 49 REAL,INTENT(OUT) :: emis(ngrid) 49 REAL,INTENT(OUT) :: albedo(ngrid,2) ! surface albedo 50 REAL,INTENT(OUT) :: emis(ngrid) ! ground emissivity 50 51 REAL,INTENT(OUT) :: q2(ngrid,nlayer+1),qsurf(ngrid,nqtot) 51 52 REAL,INTENT(OUT) :: tauscaling(ngrid) ! dust conversion factor … … 135 136 real totcloudfracS(iip1,jjp1) 136 137 real watercapS(iip1,jjp1) 138 real albedoS(iip1,jjp1) 137 139 138 140 real ptotal, co2icetotal … … 165 167 real, dimension(:,:), allocatable :: totcloudfracold 166 168 real, dimension(:,:), allocatable :: watercapold 169 real, dimension(:,:), allocatable :: albedoold 167 170 168 171 real tab_cntrl(100) … … 367 370 allocate(totcloudfracold(imold+1,jmold+1)) 368 371 allocate(watercapold(imold+1,jmold+1)) 372 allocate(albedoold(imold+1,jmold+1)) 369 373 370 374 allocate(var (imold+1,jmold+1,llm)) … … 692 696 693 697 c----------------------------------------------------------------------- 694 c 5.1 Lecture des champs 2D (co2ice, emis,ps,tsurf,Tg[10], q2surf, tauscaling,totcloudfrac,watercap)698 c 5.1 Lecture des champs 2D (co2ice,emis,ps,tsurf,Tg[10], q2surf, tauscaling,totcloudfrac,watercap,albedo) 695 699 c----------------------------------------------------------------------- 696 700 … … 824 828 IF (ierr .NE. NF_NOERR) THEN 825 829 PRINT*, "lect_start_archive: Failed loading <watercap>" 830 PRINT*, NF_STRERROR(ierr) 831 CALL abort 832 ENDIF 833 ENDIF 834 c 835 ierr = NF_INQ_VARID (nid, "albedo", nvarid) 836 IF (ierr .NE. NF_NOERR) THEN 837 PRINT*, "lect_start_archive: <albedo> not in file" 838 albedoold(:,:) = 0. 839 ELSE 840 #ifdef NC_DOUBLE 841 ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,count, 842 & albedoold) 843 #else 844 ierr = NF_GET_VARA_REAL(nid,nvarid,start,count, 845 & albedoold) 846 #endif 847 IF (ierr .NE. NF_NOERR) THEN 848 PRINT*, "lect_start_archive: Failed loading <albedo>" 826 849 PRINT*, NF_STRERROR(ierr) 827 850 CALL abort … … 1120 1143 & jjm,1,rlonuold,rlatvold,rlonu,rlatv) 1121 1144 call gr_dyn_fi (1,iim+1,jjm+1,ngrid,watercaps,watercap) 1145 1146 c Surface albedo 1147 call interp_horiz (albedoold,albedoS,imold,jmold,iim, 1148 & jjm,1,rlonuold,rlatvold,rlonu,rlatv) 1149 call gr_dyn_fi (1,iim+1,jjm+1,ngrid,albedoS,albedo(:,1)) 1150 1151 albedo(:,2)=albedo(:,1) 1122 1152 1123 1153 c write(46,*) 'emis',emis … … 1482 1512 deallocate(totcloudfracold) 1483 1513 deallocate(watercapold) 1514 deallocate(albedoold) 1484 1515 deallocate(var,varp1) 1485 1516 -
trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F
r2565 r2630 417 417 write(*,*) 'Reading file START_ARCHIVE' 418 418 CALL lect_start_archive(ngridmx,llm,nqtot, 419 & date,tsurf,tsoil, emis,q2,419 & date,tsurf,tsoil,albedo,emis,q2, 420 420 & t,ucov,vcov,ps,co2ice,teta,phisold_newgrid,q,qsurf, 421 421 & tauscaling,totcloudfrac,surfith,nid,watercap) -
trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F
r2565 r2630 89 89 REAL,ALLOCATABLE :: qsurfS(:,:) 90 90 REAL emisS(ip1jmp1) 91 REAL albedoS(ip1jmp1) 91 92 92 93 c Variables intermediaires : vent naturel, mais pas coord scalaire … … 246 247 c tsurf --> tsurfS 247 248 c co2ice --> co2iceS 249 c watercap --> watercapS 248 250 c tsoil --> tsoilS 251 c inertiedat --> ithS 249 252 c emis --> emisS 253 c albedo --> albedoS 250 254 c q2 --> q2S 251 255 c qsurf --> qsurfS … … 262 266 call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiedat,ithS) 263 267 call gr_fi_dyn(1,ngridmx,iip1,jjp1,emis,emisS) 268 call gr_fi_dyn(1,ngridmx,iip1,jjp1,albedo(:,1),albedoS) 264 269 call gr_fi_dyn(llm+1,ngridmx,iip1,jjp1,q2,q2S) 265 270 call gr_fi_dyn(nqtot,ngridmx,iip1,jjp1,qsurf,qsurfS) … … 343 348 344 349 c----------------------------------------------------------------------- 345 c Ecriture des champs (co2ice,emis, ps,Tsurf,T,u,v,q2,q,qsurf)350 c Ecriture des champs (co2ice,emis,albedo,ps,Tsurf,T,u,v,q2,q,qsurf) 346 351 c----------------------------------------------------------------------- 347 352 c ATTENTION: q2 a une couche de plus!!!! … … 361 366 & 'sub grid cloud fraction',' ',2,totcloudfracS) 362 367 call write_archive(nid,ntime,'emis','grd emis',' ',2,emisS) 368 call write_archive(nid,ntime,'albedo','surface albedo',' ', 369 & 2,albedoS) 363 370 call write_archive(nid,ntime,'ps','Psurf','Pa',2,ps) 364 371 call write_archive(nid,ntime,'tsurf','surf T','K',2,tsurfS) -
trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90
r2578 r2630 102 102 REAL :: timestart ! to pick which initial state to start from 103 103 REAL :: surfemis ! constant emissivity when no startfi 104 REAL :: surfalbedo ! constant emissivitywhen no startfi104 REAL :: surfalbedo ! constant albedo when no startfi 105 105 CHARACTER(len=5) :: modname="phyetat0" 106 106
Note: See TracChangeset
for help on using the changeset viewer.