Changeset 3453 for trunk/LMDZ.PLUTO


Ignore:
Timestamp:
Oct 9, 2024, 2:06:37 PM (6 weeks ago)
Author:
afalco
Message:

Pluto PCM: cleaned outputs.
AF

Location:
trunk/LMDZ.PLUTO/libf/phypluto
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

    r3452 r3453  
    5757                              carbox, methane,condmetsurf,condcosurf,&
    5858                              oldplutovdifc,oldplutocorrk,oldplutosedim, &
    59                               aerohaze,haze_proffix,source_haze, tsurfmax, &
    60                               albmin_ch4, &
     59                              aerohaze,haze_proffix,haze_radproffix,&
     60                              source_haze, tsurfmax, albmin_ch4, &
    6161                              season, sedimentation,generic_condensation, &
    6262                              specOLR, &
    … …  
    22552255      call writediagfi(ngrid,"tsurf","Surface temperature","K",2,tsurf)
    22562256      call writediagfi(ngrid,"ps","Surface pressure","Pa",2,ps)
     2257      call writediagfi(ngrid,"emis","Emissivity","",2,emis)
    22572258
    22582259      !! Pluto outputs
    2259       ! call writediagfi(ngrid,"rice_ch4","ch4 ice mass mean radius","m",3,rice_ch4)
    22602260      call writediagfi(ngrid,"dist_star","dist_star","AU",0,dist_star)
    22612261
    22622262      if (fast) then
    22632263         call writediagfi(ngrid,"globave","surf press","Pa",0,globave)
    2264          !AF: TODO which outputs?
     2264         call writediagfi(ngrid,"fluxrad","fluxrad","W m-2",2,fluxrad)
     2265         call writediagfi(ngrid,"fluxgrd","fluxgrd","W m-2",2,fluxgrd)
     2266         call writediagfi(ngrid,"capcal","capcal","W.s m-2 K-1",2,capcal)
     2267         ! call writediagfi(ngrid,"dplanck","dplanck","W.s m-2 K-1",2,dplanck)
     2268         call writediagfi(ngrid,"tsoil","tsoil","K",3,tsoil)
    22652269      else
    22662270         if (check_physics_outputs) then
    … …  
    22692273         endif
    22702274
    2271        call writediagfi(ngrid,"emis","Emissivity","",2,emis)
    2272        call writediagfi(ngrid,"temp","temperature","K",3,zt)
    2273        call writediagfi(ngrid,"teta","potential temperature","K",3,zh)
    2274        call writediagfi(ngrid,"u","Zonal wind","m.s-1",3,zu)
    2275        call writediagfi(ngrid,"v","Meridional wind","m.s-1",3,zv)
    2276        call writediagfi(ngrid,"w","Vertical wind","m.s-1",3,pw)
    2277        call writediagfi(ngrid,"p","Pressure","Pa",3,pplay)
     2275         call writediagfi(ngrid,"temp","temperature","K",3,zt)
     2276         call writediagfi(ngrid,"teta","potential temperature","K",3,zh)
     2277         call writediagfi(ngrid,"u","Zonal wind","m.s-1",3,zu)
     2278         call writediagfi(ngrid,"v","Meridional wind","m.s-1",3,zv)
     2279         call writediagfi(ngrid,"w","Vertical wind","m.s-1",3,pw)
     2280         call writediagfi(ngrid,"p","Pressure","Pa",3,pplay)
    22782281      endif
    22792282
    … …  
    22812284!        call writediagsoil(ngrid,"tsurf","Surface temperature","K",2,tsurf)
    22822285!        call writediagsoil(ngrid,"temp","temperature","K",3,tsoil)
    2283 
    2284 !       ! Oceanic layers !AF24: removed
    2285       ! ! Thermal plume model  !AF24: removed
    2286 !        GW non-oro outputs  !AF24: removed
    22872286
    22882287      ! Total energy balance diagnostics
    … …  
    22962295         call writediagfi(ngrid,"shad","rings"," ", 2, fract)
    22972296
    2298 !           call writediagfi(ngrid,"ASRcs","absorbed stellar rad (cs).","W m-2",2,fluxabs_sw1)
    2299 !           call writediagfi(ngrid,"OLRcs","outgoing longwave rad (cs).","W m-2",2,fluxtop_lw1)
    23002297!           call writediagfi(ngrid,"fluxsurfsw","sw surface flux.","W m-2",2,fluxsurf_sw)
    23012298!           call writediagfi(ngrid,"fluxsurflw","lw back radiation.","W m-2",2,fluxsurf_lw)
    2302 !           call writediagfi(ngrid,"fluxsurfswcs","sw surface flux (cs).","W m-2",2,fluxsurf_sw1)
    2303 !           call writediagfi(ngrid,"fluxsurflwcs","lw back radiation (cs).","W m-2",2,fluxsurf_lw1)
    23042299
    23052300         call writediagfi(ngrid,"GND","heat flux from ground","W m-2",2,fluxgrd)
    … …  
    23112306
    23122307      if(enertest) then
    2313 
    23142308         if (calldifv) then
    23152309
    … …  
    23282322         endif
    23292323
    2330 !          if(watercond) then  !AF24: removed
    2331 
    23322324         if (generic_condensation) then
    2333 
    23342325            call writediagfi(ngrid,"genericconddE","heat from generic condensation","W m-2",2,genericconddE)
    23352326            call writediagfi(ngrid,"dt_generic_condensation","heating from generic condensation","K s-1",3,dt_generic_condensation)
    2336 
    23372327         endif
    23382328
    … …  
    23522342        endif
    23532343
    2354 
    23552344        ! Temporary inclusions for heating diagnostics.
    23562345        call writediagfi(ngrid,"zdtsw","SW heating","T s-1",3,zdtsw)
    … …  
    23592348        call writediagfi(ngrid,"zdtdyn","Dyn. heating","T s-1",3,zdtdyn)
    23602349
    2361         !Pluto specific
    2362       !  call writediagfi(ngrid,"zdtc","tendancy T cond N2","K",3,zdtc)
    2363       !   call writediagfi(ngrid,"zdtch4cloud","tendancy T ch4cloud","K",3,zdtch4cloud)
    2364       !   call writediagfi(ngrid,"zdtcocloud","tendancy T cocloud","K",3,zdtcocloud)
    2365       !   call writediagfi(ngrid,"zq1temp_ch4"," "," ",2,zq1temp_ch4)
    2366       !   call writediagfi(ngrid,"qsat_ch4"," "," ",2,qsat_ch4)
    2367       !   call writediagfi(ngrid,"qsat_ch4_l1"," "," ",2,qsat_ch4_l1)
    2368       !   call writediagfi(ngrid,"senshf1","senshf1"," ",2,sensiblehf1)
    2369       !   call writediagfi(ngrid,"senshf2","senshf2"," ",2,sensiblehf2)
    2370 
    2371 
    23722350        ! For Debugging.
    23732351        !call writediagfi(ngrid,'rnat','Terrain type',' ',2,real(rnat))
    23742352        !call writediagfi(ngrid,'pphi','Geopotential',' ',3,pphi)
    23752353
    2376 
    2377       ! Output aerosols.!AF: TODO: write haze aerosols
    2378       ! if (igcm_n2_ice.ne.0.and.iaero_haze.ne.0) &
    2379       !    call writediagfi(ngrid,'N2ice_reff','N2ice_reff','m',3,reffrad(1,1,iaero_haze))
    2380       ! if (igcm_n2_ice.ne.0.and.iaero_haze.ne.0) &
    2381       !    call writediagfi(ngrid,'N2ice_reffcol','N2ice_reffcol','um kg m^-2',2,reffcol(1,iaero_haze))
    2382       ! if (igcm_h2o_ice.ne.0.and.iaero_h2o.ne.0) &  !AF24: removed
    2383 
    23842354      ! Output tracers.
    23852355      if (tracer) then
     2356         ! call writediagfi(ngrid,"zdtc","tendancy T cond N2","K",3,zdtc)
    23862357
    23872358         do iq=1,nq
    2388           call writediagfi(ngrid,noms(iq),noms(iq),'kg/kg',3,zq(1,1,iq))
    2389          !  call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
     2359            call writediagfi(ngrid,noms(iq),noms(iq),'kg/kg',3,zq(1,1,iq))
     2360            ! call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
    23902361         !                    'kg m^-2',2,qsurf_hist(1,iq) )
    23912362           call writediagfi(ngrid,trim(noms(iq))//'_col',trim(noms(iq))//'_col',    &
    23922363                           'kg m^-2',2,qcol(1,iq) )
    2393          call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
     2364            call writediagfi(ngrid,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
    23942365                         'kg m^-2',2,qsurf(1,iq) )
    2395 
    2396             ! if(watercond.or.CLFvarying)then  !AF24: removed
    23972366
    23982367            if(generic_condensation)then
    … …  
    24012370               call writediagfi(ngrid,"RH_generic","GCS relative humidity"," ",3,RH_generic)
    24022371            endif
    2403 
    2404             ! if(generic_rain)then  !AF24: removed
    2405             ! if((hydrology).and.(.not.ok_slab_ocean))then  !AF24: removed
    2406 
    24072372           ! call writediagfi(ngrid,"tau_col","Total aerosol optical depth","[]",2,tau_col)
    2408 
    24092373         enddo ! end of 'nq' loop
    24102374
    24112375         !Pluto specific
    24122376         call writediagfi(ngrid,'n2_iceflux','n2_iceflux',"kg m^-2 s^-1",2,flusurf(1,igcm_n2) )
    2413          if (aerohaze)then
     2377         if (haze_radproffix)then
    24142378            call writediagfi(ngrid,'haze_reff','haze_reff','m',3,reffrad(1,1,1))
    24152379         end if
    24162380         if (methane) then
     2381            ! call writediagfi(ngrid,"rice_ch4","ch4 ice mass mean radius","m",3,rice_ch4)
     2382            ! call writediagfi(ngrid,"zq1temp_ch4"," "," ",2,zq1temp_ch4)
     2383            ! call writediagfi(ngrid,"qsat_ch4"," "," ",2,qsat_ch4)
     2384            ! call writediagfi(ngrid,"qsat_ch4_l1"," "," ",2,qsat_ch4_l1)
     2385
    24172386            call writediagfi(ngrid,'ch4_iceflux','ch4_iceflux',&
    24182387                              "kg m^-2 s^-1",2,flusurf(1,igcm_ch4_ice) )
    … …  
    24232392
    24242393            ! Tendancies
    2425             !call writediagfi(ngrid,"zdqch4cloud","ch4 cloud","T s-1",&
    2426             !               3,zdqch4cloud(1,1,igcm_ch4_gas))
    2427             ! call writediagfi(ngrid,"zdqcn2_ch4","zdq condn2 ch4","",&
    2428             !                3,zdqc(:,:,igcm_ch4_gas))
    2429             ! call writediagfi(ngrid,"zdqdif_ch4","zdqdif ch4","",&
    2430             !                3,zdqdif(:,:,igcm_ch4_gas))
    2431             ! call writediagfi(ngrid,"zdqsdif_ch4_ice","zdqsdif ch4","",&
    2432             !                2,zdqsdif(:,igcm_ch4_ice))
    2433             ! call writediagfi(ngrid,"zdqadj_ch4","zdqadj ch4","",&
    2434             !                3,zdqadj(:,:,igcm_ch4_gas))
     2394            call writediagfi(ngrid,"zdqcn2_ch4","zdq condn2 ch4","",&
     2395                           3,zdqc(:,:,igcm_ch4_gas))
     2396            call writediagfi(ngrid,"zdqdif_ch4","zdqdif ch4","",&
     2397                           3,zdqdif(:,:,igcm_ch4_gas))
     2398            call writediagfi(ngrid,"zdqsdif_ch4_ice","zdqsdif ch4","",&
     2399                           2,zdqsdif(:,igcm_ch4_ice))
     2400            call writediagfi(ngrid,"zdqadj_ch4","zdqadj ch4","",&
     2401                           3,zdqadj(:,:,igcm_ch4_gas))
    24352402             if (sedimentation) then
    24362403               call writediagfi(ngrid,"zdqsed_ch4","zdqsed ch4","",&
    … …  
    24492416
    24502417         if (carbox) then
     2418            ! call writediagfi(ngrid,"zdtcocloud","tendancy T cocloud","K",3,zdtcocloud)
    24512419            call writediagfi(ngrid,'co_iceflux','co_iceflux',&
    24522420                               "kg m^-2 s^-1",2,flusurf(1,igcm_co_ice) )
    … …  
    25032471   comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid)
    25042472   if (.not.calldifv) comm_LATENT_HF(:)=0.0
    2505    ! if ((tracer).and.(water)) then  !AF24: removed
    25062473
    25072474   if ((tracer).and.(generic_condensation)) then
    2508    ! .and.(.not. water)
    25092475
    25102476      ! If you have set generic_condensation (and not water) and you have set several GCS
    … …  
    25572523   sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid) !!! ????
    25582524   comm_HR_DYN(1:ngrid,1:nlayer) = zdtdyn(1:ngrid,1:nlayer)
    2559 
    2560 !      if (turb_resolved) then
    2561 !        open(17,file='lsf.txt',form='formatted',status='old')
    2562 !        rewind(17)
    2563 !        DO l=1,nlayer
    2564 !          read(17,*) lsf_dt(l),lsf_dq(l)
    2565 !        ENDDO
    2566 !        close(17)
    2567 !        do ig=1,ngrid
    2568 !          if ((tracer).and.(water)) then
    2569 !           pdq(ig,:,igcm_h2o_gas) = pdq(ig,:,igcm_h2o_gas) + lsf_dq(:)
    2570 !          endif
    2571 !          pdt(ig,:) = pdt(ig,:) + lsf_dt(:)
    2572 !          comm_HR_DYN(ig,:) = lsf_dt(:)
    2573 !        enddo
    2574 !      endif
    25752525#endif
    25762526
    … …  
    25872537      CALL send_xios_field("u",zu)
    25882538      CALL send_xios_field("v",zv)
    2589       CALL send_xios_field("omega",omega)
    2590 
    2591 !       IF (calltherm) THEN  !AF24: removed
    2592       ! IF (water) THEN  !AF24: removed
     2539      call send_xios_field("w",pw)
    25932540
    25942541      CALL send_xios_field("ISR",fluxtop_dn)

  • trunk/LMDZ.PLUTO/libf/phypluto/tabfi_mod.F90

    r3184 r3453  
    192192! boundary layer and turbulence
    193193      z0 = tab_cntrl(tab0+19)
    194       lmixmin = tab_cntrl(tab0+20)
    195       emin_turb = tab_cntrl(tab0+21)
     194! for paleoclimate
     195      tpal = tab_cntrl(tab0+20)
     196      adjust = tab_cntrl(tab0+21) ! for Triton albedo adjustment
     197      ! lmixmin = tab_cntrl(tab0+20)
     198      ! emin_turb = tab_cntrl(tab0+21)
     199
     200!!! AF24: parameters below are not used?
    196201! optical properties of polar caps and ground emissivity
    197202      emisice(1) = tab_cntrl(tab0+24)
Note: See TracChangeset for help on using the changeset viewer.