Changeset 2865 for trunk/LMDZ.GENERIC
- Timestamp:
- Jan 16, 2023, 4:47:05 PM (2 years ago)
- Location:
- trunk/LMDZ.GENERIC/libf/phystd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/comm_wrf.F90
r2025 r2865 20 20 REAL,SAVE,ALLOCATABLE :: comm_CLOUDFRAC(:,:) 21 21 REAL,SAVE,ALLOCATABLE :: comm_TOTCLOUDFRAC(:) 22 REAL,SAVE,ALLOCATABLE :: comm_RAIN(:,:) 23 REAL,SAVE,ALLOCATABLE :: comm_SNOW(:,:) 22 REAL,SAVE,ALLOCATABLE :: comm_RH(:,:) 23 REAL,SAVE,ALLOCATABLE :: comm_DQICE(:,:) 24 REAL,SAVE,ALLOCATABLE :: comm_DQVAP(:,:) 24 25 REAL,SAVE,ALLOCATABLE :: comm_ALBEQ(:) 25 26 REAL,SAVE,ALLOCATABLE :: comm_FLUXTOP_DN(:) … … 29 30 REAL,SAVE,ALLOCATABLE :: comm_FLUXSURF_LW(:) 30 31 REAL,SAVE,ALLOCATABLE :: comm_FLXGRD(:) 31 REAL,SAVE,ALLOCATABLE :: comm_LSCEZ(:,:) 32 REAL,SAVE,ALLOCATABLE :: comm_DTRAIN(:,:) 33 REAL,SAVE,ALLOCATABLE :: comm_DTLSC(:,:) 32 34 REAL,SAVE,ALLOCATABLE :: comm_H2OICE_REFF(:,:) 33 35 REAL,SAVE,ALLOCATABLE :: comm_LATENT_HF(:) … … 56 58 allocate(comm_CLOUDFRAC(ngrid,nlayer)) 57 59 allocate(comm_TOTCLOUDFRAC(ngrid)) 58 allocate(comm_RAIN(ngrid,nlayer)) 59 allocate(comm_SNOW(ngrid,nlayer)) 60 allocate(comm_RH(ngrid,nlayer)) 61 allocate(comm_DQICE(ngrid,nlayer)) 62 allocate(comm_DQVAP(ngrid,nlayer)) 60 63 allocate(comm_ALBEQ(ngrid)) 61 64 allocate(comm_FLUXTOP_DN(ngrid)) … … 65 68 allocate(comm_FLUXSURF_LW(ngrid)) 66 69 allocate(comm_FLXGRD(ngrid)) 67 allocate(comm_LSCEZ(ngrid,nlayer)) 70 allocate(comm_DTRAIN(ngrid,nlayer)) 71 allocate(comm_DTLSC(ngrid,nlayer)) 68 72 allocate(comm_H2OICE_REFF(ngrid,nlayer)) 69 73 allocate(comm_LATENT_HF(ngrid)) … … 90 94 deallocate(comm_CLOUDFRAC) 91 95 deallocate(comm_TOTCLOUDFRAC) 92 deallocate(comm_RAIN) 93 deallocate(comm_SNOW) 96 deallocate(comm_RH) 97 deallocate(comm_DQICE) 98 deallocate(comm_DQVAP) 94 99 deallocate(comm_ALBEQ) 95 100 deallocate(comm_FLUXTOP_DN) … … 99 104 deallocate(comm_FLUXSURF_LW) 100 105 deallocate(comm_FLXGRD) 101 deallocate(comm_LSCEZ) 106 deallocate(comm_DTRAIN) 107 deallocate(comm_DTLSC) 102 108 deallocate(comm_H2OICE_REFF) 103 109 deallocate(comm_LATENT_HF) -
trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90
r2839 r2865 78 78 use mod_phys_lmdz_omp_data, ONLY: is_omp_master 79 79 #else 80 use comm_wrf, only : comm_HR_SW, comm_HR_LW, &81 comm_CLOUDFRAC,comm_TOTCLOUDFRAC,&82 comm_SURFRAIN,comm_REEVAP,comm_HR_DYN,&83 comm_RAIN,comm_SNOW,comm_ALBEQ,&84 comm_FLUXTOP_DN,comm_FLUXABS_SW,&85 comm_FLUXTOP_LW,comm_FLUXSURF_SW,&86 comm_FLUXSURF_LW,comm_FLXGRD,&87 comm_LSCEZ,comm_H2OICE_REFF80 use comm_wrf, only : comm_HR_SW, comm_HR_LW, & 81 comm_CLOUDFRAC,comm_TOTCLOUDFRAC, comm_RH, & 82 comm_SURFRAIN,comm_REEVAP,comm_HR_DYN, & 83 comm_DQICE,comm_DQVAP,comm_ALBEQ, & 84 comm_FLUXTOP_DN,comm_FLUXABS_SW, & 85 comm_FLUXTOP_LW,comm_FLUXSURF_SW, & 86 comm_FLUXSURF_LW,comm_FLXGRD, & 87 comm_DTRAIN,comm_DTLSC,comm_H2OICE_REFF 88 88 #endif 89 89 … … 2572 2572 2573 2573 #else 2574 comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer) 2575 comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer) 2576 if ((tracer).and.(water)) then 2577 comm_CLOUDFRAC(1:ngrid,1:nlayer)=cloudfrac(1:ngrid,1:nlayer) 2578 comm_TOTCLOUDFRAC(1:ngrid)=totcloudfrac(1:ngrid) 2579 comm_RAIN(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_vap) 2580 comm_SURFRAIN(1:ngrid)=zdqsrain(1:ngrid) 2581 comm_SNOW(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_ice) 2582 comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid) 2583 comm_H2OICE_REFF(1:ngrid,1:nlayer)=reffrad(1:ngrid,1:nlayer,iaero_h2o) 2584 comm_REEVAP(1:ngrid)=reevap_precip(1:ngrid) 2585 endif 2586 comm_FLUXTOP_DN(1:ngrid)=fluxtop_dn(1:ngrid) 2587 comm_FLUXABS_SW(1:ngrid)=fluxabs_sw(1:ngrid) 2588 comm_FLUXTOP_LW(1:ngrid)=fluxtop_lw(1:ngrid) 2589 comm_FLUXSURF_SW(1:ngrid)=fluxsurf_sw(1:ngrid) 2590 comm_FLUXSURF_LW(1:ngrid)=fluxsurf_lw(1:ngrid) 2591 comm_FLXGRD(1:ngrid)=fluxgrd(1:ngrid) 2592 comm_LSCEZ(1:ngrid,1:nlayer)=lscaledEz(1:ngrid,1:nlayer) 2593 sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid) !!! ???? 2574 comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer) 2575 comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer) 2576 if ((tracer).and.(water)) then 2577 comm_CLOUDFRAC(1:ngrid,1:nlayer)=cloudfrac(1:ngrid,1:nlayer) 2578 comm_TOTCLOUDFRAC(1:ngrid)=totcloudfrac(1:ngrid) 2579 comm_SURFRAIN(1:ngrid)=zdqsrain(1:ngrid) 2580 comm_DQVAP(1:ngrid,1:nlayer)=pdq(1:ngrid,1:nlayer,igcm_h2o_vap) 2581 comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid) 2582 comm_H2OICE_REFF(1:ngrid,1:nlayer)=reffrad(1:ngrid,1:nlayer,iaero_h2o) 2583 comm_REEVAP(1:ngrid)=reevap_precip(1:ngrid) 2584 comm_DTRAIN(1:ngrid,1:nlayer)=zdtrain(1:ngrid,1:nlayer) 2585 comm_DTLSC(1:ngrid,1:nlayer)=dtlscale(1:ngrid,1:nlayer) 2586 comm_RH(1:ngrid,1:nlayer)=RH(1:ngrid,1:nlayer) 2587 endif 2588 !comm_DQICE(1:ngrid,1:nlayer)=zdqdyn(1:ngrid,1:nlayer) 2589 comm_FLUXTOP_DN(1:ngrid)=fluxtop_dn(1:ngrid) 2590 comm_FLUXABS_SW(1:ngrid)=fluxabs_sw(1:ngrid) 2591 comm_FLUXTOP_LW(1:ngrid)=fluxtop_lw(1:ngrid) 2592 comm_FLUXSURF_SW(1:ngrid)=fluxsurf_sw(1:ngrid) 2593 comm_FLUXSURF_LW(1:ngrid)=fluxsurf_lw(1:ngrid) 2594 comm_FLXGRD(1:ngrid)=fluxgrd(1:ngrid) 2595 sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid) !!! ???? 2596 comm_HR_DYN(1:ngrid,1:nlayer) = zdtdyn(1:ngrid,1:nlayer) 2594 2597 2595 2598 if (turb_resolved) then
Note: See TracChangeset
for help on using the changeset viewer.