Changeset 4031 for LMDZ6/trunk


Ignore:
Timestamp:
Nov 29, 2021, 3:29:39 PM (3 years ago)
Author:
idelkadi
Message:

Corrections in the LMDZ-ECRAN interface:

  • gas concentrations in VMR instead of MMR
  • use of a constant value for each gas instead of a table of values
  • specify in the 14 SW bands of ECRAD those corresponding to the 6 SW bands for diffuse and direct albedo in LMDZ
  • passage of all processing of input fields for Ecrad (conversions, calculations, ...) from radiation.F90 to radlw.F90
Location:
LMDZ6/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/ecrad/radiation_scheme.F90

    r3946 r4031  
    2222     &  PPRESSURE, PTEMPERATURE, &
    2323     &  PPRESSURE_H, PTEMPERATURE_H, PQ, PQSAT, &
    24      &  PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4, PO3_DP, &
     24     &  PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, &
     25     &  PCCL4, PO3, PO2, &
    2526     &  PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW, &
    2627     &  ZRE_LIQUID_UM, ZRE_ICE_UM, &
     
    151152! AI
    152153REAL(KIND=JPRB),   INTENT(IN) :: PQSAT(KLON,KLEV)
    153 REAL(KIND=JPRB),   INTENT(IN) :: PCO2(KLON,KLEV)
    154 REAL(KIND=JPRB),   INTENT(IN) :: PCH4(KLON,KLEV)
    155 REAL(KIND=JPRB),   INTENT(IN) :: PN2O(KLON,KLEV)
    156 REAL(KIND=JPRB),   INTENT(IN) :: PNO2(KLON,KLEV)
    157 REAL(KIND=JPRB),   INTENT(IN) :: PCFC11(KLON,KLEV)
    158 REAL(KIND=JPRB),   INTENT(IN) :: PCFC12(KLON,KLEV)
    159 REAL(KIND=JPRB),   INTENT(IN) :: PHCFC22(KLON,KLEV)
    160 REAL(KIND=JPRB),   INTENT(IN) :: PCCL4(KLON,KLEV)
    161 REAL(KIND=JPRB),   INTENT(IN) :: PO3_DP(KLON,KLEV) ! AI (kg/kg) ATTENTION (Pa*kg/kg)
     154REAL(KIND=JPRB),   INTENT(IN) :: PCO2
     155REAL(KIND=JPRB),   INTENT(IN) :: PCH4
     156REAL(KIND=JPRB),   INTENT(IN) :: PN2O
     157REAL(KIND=JPRB),   INTENT(IN) :: PNO2
     158REAL(KIND=JPRB),   INTENT(IN) :: PCFC11
     159REAL(KIND=JPRB),   INTENT(IN) :: PCFC12
     160REAL(KIND=JPRB),   INTENT(IN) :: PHCFC22
     161REAL(KIND=JPRB),   INTENT(IN) :: PCCL4
     162REAL(KIND=JPRB),   INTENT(IN) :: PO3(KLON,KLEV) ! AI (kg/kg) ATTENTION (Pa*kg/kg)
     163REAL(KIND=JPRB),   INTENT(IN) :: PO2
    162164
    163165! *** Cloud fraction and hydrometeor mass mixing ratios
     
    326328  print*,'PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4 =', &
    327329        PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4
    328   print*,'PO3_DP =',PO3_DP
     330  print*,'PO3 =',PO3
    329331  print*,'PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW =', &
    330332        PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW
     
    345347! AI appel radiation_setup
    346348call SETUP_RADIATION_SCHEME(loutput)
     349!! Les 6 bandes SW pour l'albedo :
     350!! 0.185-0.25, 0.25-0.4, 0.4-0.69 , 0.69-1.19, 1.19-2.38, 2.38-4.00 micro-metre
     351  call rad_config%define_sw_albedo_intervals(6, &
     352             &  (/ 0.25e-6_jprb, 0.44e-6_jprb, 1.19e-6_jprb, &
     353             &     2.38e-6_jprb, 4.00e-6_jprb /),  (/ 1,2,3,4,5,6 /))
    347354
    348355if (lprint_config) then
     
    372379  print*,'rad_config%i_emiss_from_band_lw =', rad_config%i_emiss_from_band_lw
    373380endif
    374 !stop
    375 ! A EFFACER
    376 !print*,'n_g_lw, n_g_sw =', rad_config%n_g_lw, rad_config%n_g_sw
    377 !print*,'use_canopy_full_spectrum_lw = ', rad_config%use_canopy_full_spectrum_lw
    378 !print*,'rad_config%i_band_from_reordered_g_lw =', &
    379 !       rad_config%i_band_from_reordered_g_lw
    380 !print*,'use_canopy_full_spectrum_lw =', rad_config%use_canopy_full_spectrum_lw
    381 !rad_config%use_canopy_full_spectrum_lw = .TRUE.
    382 ! AI ATTENTION
    383 !rad_config%i_band_from_reordered_g_lw = 1
    384 !rad_config%use_spectral_solar_scaling = .true.
    385 !endif
    386 ! AI ATTENTION test
    387 !rad_config%i_gas_model = IGasModelMonochromatic
    388381
    389382! AI ATTENTION
    390383! Allocate memory in radiation objects
    391 CALL single_level%allocate(KLON, NSW, 2, &
     384CALL single_level%allocate(KLON, NSW, 1, &
    392385     &                     use_sw_albedo_direct=.TRUE.)
    393386
     
    420413! temperature at the half-level corresponding to the surface as
    421414! follows:
    422 thermodynamics%temperature_hl(KIDIA:KFDIA,KLEV+1) &
    423      &  = PTEMPERATURE(KIDIA:KFDIA,KLEV) &
    424      &  + 0.5_JPRB * (PTEMPERATURE_H(KIDIA:KFDIA,KLEV+1) &
    425      &               -PTEMPERATURE_H(KIDIA:KFDIA,KLEV))
     415!thermodynamics%temperature_hl(KIDIA:KFDIA,KLEV+1) &
     416!     &  = PTEMPERATURE(KIDIA:KFDIA,KLEV) &
     417!     &  + 0.5_JPRB * (PTEMPERATURE_H(KIDIA:KFDIA,KLEV+1) &
     418!     &               -PTEMPERATURE_H(KIDIA:KFDIA,KLEV))
    426419
    427420! Alternatively we respect the model's atmospheric temperature in the
     
    461454!single_level%lw_emissivity(KIDIA:KFDIA,1)  = 1.0_JPRB
    462455single_level%lw_emissivity(KIDIA:KFDIA,1)  = PEMIS(KIDIA:KFDIA)
    463 single_level%lw_emissivity(KIDIA:KFDIA,2)  = PEMIS_WINDOW(KIDIA:KFDIA)
     456!single_level%lw_emissivity(KIDIA:KFDIA,2)  = PEMIS_WINDOW(KIDIA:KFDIA)
    464457
    465458! Create the relevant seed from date and time get the starting day
     
    499492!     &  PLAND_SEA_MASK, PCCN_LAND, PCCN_SEA, &
    500493!     &  ZRE_LIQUID_UM)
    501 cloud%re_liq(KIDIA:KFDIA,:) = ZRE_LIQUID_UM(KIDIA:KFDIA,:) * 1.0e-6_JPRB
     494cloud%re_liq(KIDIA:KFDIA,:) = ZRE_LIQUID_UM(KIDIA:KFDIA,:)
    502495
    503496!CALL ICE_EFFECTIVE_RADIUS(KIDIA, KFDIA, KLON, KLEV, &
    504497!     &  PPRESSURE, PTEMPERATURE, PCLOUD_FRAC, PQ_ICE, PQ_SNOW, PGEMU, &
    505498!     &  ZRE_ICE_UM)
    506 cloud%re_ice(KIDIA:KFDIA,:) = ZRE_ICE_UM(KIDIA:KFDIA,:) * 1.0e-6_JPRB
     499cloud%re_ice(KIDIA:KFDIA,:) = ZRE_ICE_UM(KIDIA:KFDIA,:)
    507500
    508501! Get the cloud overlap decorrelation length (for cloud boundaries),
     
    539532  ENDDO
    540533ENDDO
     534!AI ATTENTION meme traitement dans le version offline
     535!call cloud%create_inv_cloud_effective_size_eta(ncol, nlev, &
     536!               &  thermodynamics%pressure_hl, &
     537!               &  low_inv_effective_size, &
     538!               &  middle_inv_effective_size, &
     539!               &  high_inv_effective_size, 0.8_jprb, 0.45_jprb)     
    541540
    542541print*,'******** AEROSOLS (allocate + input) **************************************'
     
    608607!  ENDDO
    609608!ENDDO
    610 ZO3 = PO3_DP
    611609
    612610!  Insert gas mixing ratios
    613611print*,'Insert gas mixing ratios'
    614612CALL gas%put(IH2O,    IMassMixingRatio, PQ)
    615 CALL gas%put(ICO2,    IMassMixingRatio, PCO2)
    616 CALL gas%put(ICH4,    IMassMixingRatio, PCH4)
    617 CALL gas%put(IN2O,    IMassMixingRatio, PN2O)
    618 CALL gas%put(ICFC11,  IMassMixingRatio, PCFC11)
    619 CALL gas%put(ICFC12,  IMassMixingRatio, PCFC12)
    620 CALL gas%put(IHCFC22, IMassMixingRatio, PHCFC22)
    621 CALL gas%put(ICCL4,   IMassMixingRatio, PCCL4)
    622 CALL gas%put(IO3,     IMassMixingRatio, ZO3)
    623 CALL gas%put_well_mixed(IO2, IVolumeMixingRatio, 0.20944_JPRB)
     613CALL gas%put(IO3,     IMassMixingRatio, PO3)
     614CALL gas%put_well_mixed(ICO2,    IVolumeMixingRatio, PCO2)
     615CALL gas%put_well_mixed(ICH4,    IVolumeMixingRatio, PCH4)
     616CALL gas%put_well_mixed(IN2O,    IVolumeMixingRatio, PN2O)
     617CALL gas%put_well_mixed(ICFC11,  IVolumeMixingRatio, PCFC11)
     618CALL gas%put_well_mixed(ICFC12,  IVolumeMixingRatio, PCFC12)
     619CALL gas%put_well_mixed(IHCFC22, IVolumeMixingRatio, PHCFC22)
     620CALL gas%put_well_mixed(ICCL4,   IVolumeMixingRatio, PCCL4)
     621CALL gas%put_well_mixed(IO2,     IVolumeMixingRatio, PO2)
    624622! Ensure the units of the gas mixing ratios are what is required by
    625623! the gas absorption model
  • LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

    r3989 r4031  
    309309  REAL(KIND=8) ZGELAM(klon), &              ! longitudes en rad
    310310               ZGEMU(klon)                  ! sin(latitude)
    311   REAL(KIND=8) ZCO2(klon,klev), &           ! CO2 mass mixing ratios on full levels
    312                ZCH4(klon,klev), &           ! CH4 mass mixing ratios on full levels
    313                ZN2O(klon,klev), &           ! N2O mass mixing ratios on full levels
    314                ZNO2(klon,klev), &           ! NO2 mass mixing ratios on full levels
    315                ZCFC11(klon,klev), &         ! CFC11
    316                ZCFC12(klon,klev), &         ! CFC12
    317                ZHCFC22(klon,klev), &        ! HCFC22
    318                ZCCL4(klon,klev)           ! CCL4
    319 !               ZO3_DP(klon,klev), ZO3_DP_i(klon,klev)            ! Ozone
     311  REAL(KIND=8) ZCO2, &           ! CO2 mass mixing ratios on full levels
     312               ZCH4, &           ! CH4 mass mixing ratios on full levels
     313               ZN2O, &           ! N2O mass mixing ratios on full levels
     314               ZNO2, &           ! NO2 mass mixing ratios on full levels
     315               ZCFC11, &         ! CFC11
     316               ZCFC12, &         ! CFC12
     317               ZHCFC22, &        ! HCFC22
     318               ZCCL4, &          ! CCL4
     319               ZO2               ! O2
     320
    320321  REAL(KIND=8) ZQ_RAIN(klon,klev), &        ! Rain cloud mass mixing ratio (kg/kg) ?
    321322               ZQ_SNOW(klon,klev)           ! Snow cloud mass mixing ratio (kg/kg) ?
     
    486487
    487488  !
     489! AI 02.2021
     490#ifdef CPP_ECRAD
     491  ZEMIS = 1.0
     492  ZEMISW = 1.0
     493  ZGELAM = longitude
     494  ZGEMU = sin(latitude)
     495  ZCO2 = RCO2
     496  ZCH4 = RCH4
     497  ZN2O = RN2O
     498  ZNO2 = 0.0
     499  ZCFC11 = RCFC11
     500  ZCFC12 = RCFC12
     501  ZHCFC22 = 0.0
     502  ZO2 = 0.0
     503  ZCCL4 = 0.0
     504  ZQ_RAIN = 0.0
     505  ZQ_SNOW = 0.0
     506  ZAEROSOL_OLD = 0.0
     507  ZAEROSOL = 0.0
     508#endif
     509
    488510  !-------------------------------------------
    489511  nb_gr = KLON / kdlon
     
    627649      ENDDO
    628650    ENDDO
    629 !
    630 ! AI 02.2021
    631 #ifdef CPP_ECRAD
    632   ZEMIS = 1.0
    633   ZEMISW = 1.0
    634   ZGELAM = longitude
    635   ZGEMU = sin(latitude)
    636   ZCO2 = RCO2
    637   ZCH4 = RCH4
    638   ZN2O = RN2O
    639   ZNO2 = 0.0
    640   ZCFC11 = RCFC11
    641   ZCFC12 = RCFC12
    642   ZHCFC22 = 0.0
    643   ZCCL4 = 0.0
    644   ZQ_RAIN = 0.0
    645   ZQ_SNOW = 0.0
    646   ZAEROSOL_OLD = 0.0
    647   ZAEROSOL = 0.0
    648 #endif
    649651!
    650652!===== iflag_rrtm ================================================
     
    12241226            flwc_i(1:klon,k)    =flwc(1:klon,klev+1-k)
    12251227            fiwc_i(1:klon,k)    =fiwc(1:klon,klev+1-k)
    1226             ref_liq_i(1:klon,k) =ref_liq(1:klon,klev+1-k)
    1227             ref_ice_i(1:klon,k) =ref_ice(1:klon,klev+1-k)
     1228            ref_liq_i(1:klon,k) =ref_liq(1:klon,klev+1-k)*1.0e-6
     1229            ref_ice_i(1:klon,k) =ref_ice(1:klon,klev+1-k)*1.0e-6
    12281230!-OB
    12291231            ref_liq_pi_i(1:klon,k) =ref_liq_pi(1:klon,klev+1-k)
     
    12371239!            ENDDO
    12381240         ENDDO
    1239 ! AI 02.2021
     1241
     1242! AI 11.2021
    12401243! Calcul de ZTH_i (temp aux interfaces 1:klev+1)
     1244! IFS currently sets the half-level temperature at the surface to be
     1245! equal to the skin temperature. The radiation scheme takes as input
     1246! only the half-level temperatures and assumes the Planck function to
     1247! vary linearly in optical depth between half levels. In the lowest
     1248! atmospheric layer, where the atmospheric temperature can be much
     1249! cooler than the skin temperature, this can lead to significant
     1250! differences between the effective temperature of this lowest layer
     1251! and the true value in the model.
     1252! We may approximate the temperature profile in the lowest model level
     1253! as piecewise linear between the top of the layer T[k-1/2], the
     1254! centre of the layer T[k] and the base of the layer Tskin.  The mean
     1255! temperature of the layer is then 0.25*T[k-1/2] + 0.5*T[k] +
     1256! 0.25*Tskin, which can be achieved by setting the atmospheric
     1257! temperature at the half-level corresponding to the surface as
     1258! follows:
     1259! AI ATTENTION fais dans interface radlw
     1260!thermodynamics%temperature_hl(KIDIA:KFDIA,KLEV+1) &
     1261!     &  = PTEMPERATURE(KIDIA:KFDIA,KLEV) &
     1262!     &  + 0.5_JPRB * (PTEMPERATURE_H(KIDIA:KFDIA,KLEV+1) &
     1263!     &               -PTEMPERATURE_H(KIDIA:KFDIA,KLEV))
     1264
    12411265         DO K=2,KLEV
    1242             ZTH_i(:,K)=&
    1243               & (t_i(:,K-1)*pplay_i(:,K-1)*(pplay_i(:,K)-paprs_i(:,K))&
    1244               & +t_i(:,K)*pplay_i(:,K)*(paprs_i(:,K)-pplay_i(:,K-1)))&
    1245               & *(1.0/(paprs_i(:,K)*(pplay_i(:,K)-pplay_i(:,K-1))))
     1266          DO i = 1, kdlon
     1267            ZTH_i(i,K)=&
     1268              & (t_i(i,K-1)*pplay_i(i,K-1)*(pplay_i(i,K)-paprs_i(i,K))&
     1269              & +t_i(i,K)*pplay_i(i,K)*(paprs_i(i,K)-pplay_i(i,K-1)))&
     1270              & *(1.0/(paprs_i(i,K)*(pplay_i(i,K)-pplay_i(i,K-1))))
     1271           ENDDO
    12461272         ENDDO
    1247             ZTH_i(:,KLEV+1)=tsol(:)
    1248             ZTH_i(:,1)=t_i(:,1)-pplay_i(:,1)*(t_i(:,1)-ZTH_i(:,2))&
    1249                       & /(pplay_i(:,1)-paprs_i(:,2))
     1273         DO i = 1, kdlon
     1274! Sommet
     1275            ZTH_i(i,1)=t_i(i,1)-pplay_i(i,1)*(t_i(i,1)-ZTH_i(i,2))&
     1276                      & /(pplay_i(i,1)-paprs_i(i,2))
     1277! Vers le sol
     1278            ZTH_i(i,KLEV+1)=t_i(i,KLEV) + 0.5 * &
     1279                            (tsol(i) - ZTH_i(i,KLEV))
     1280         ENDDO
     1281
     1282! AI ATTENTION TESTS
     1283!      PALBD_NEW = 0.0
     1284!      PALBP_NEW = 0.0
     1285!      ZCO2 = RCO2
     1286!      ZCH4 = RCH4
     1287!      ZN2O = RN2O
     1288!      ZNO2 = 0.0
     1289!      ZCFC11 = RCFC11
     1290!      ZCFC12 = RCFC12
     1291!      ZHCFC22 = 0.0
     1292!      ZO2 = 0.0
     1293!      ZCCL4 = 0.0
    12501294
    12511295      print *,'RADLWSW: avant RADIATION_SCHEME '
     1296!      print*,'RCFC11=',RCFC11
     1297!      print*,'RCFC12=',RCFC12
     1298
    12521299      IF (lldebug) THEN
    12531300        CALL writefield_phy('rmu0',rmu0,1)
    12541301        CALL writefield_phy('tsol',tsol,1)
    12551302        CALL writefield_phy('emissiv_out',ZEMIS,1)
    1256         CALL writefield_phy('emissiv_in',ZEMISW,1)
    1257         CALL writefield_phy('pctsrf_ter',pctsrf(:,is_ter),1)
    1258         CALL writefield_phy('pctsrf_oce',pctsrf(:,is_oce),1)
    1259         CALL writefield_phy('ZGELAM',ZGELAM,1)
    1260         CALL writefield_phy('ZGEMU',ZGEMU,1)
    1261         CALL writefield_phy('zmasq',zmasq,1)
     1303!        CALL writefield_phy('emissiv_in',ZEMISW,1)
     1304!        CALL writefield_phy('pctsrf_ter',pctsrf(:,is_ter),1)
     1305!        CALL writefield_phy('pctsrf_oce',pctsrf(:,is_oce),1)
     1306!        CALL writefield_phy('ZGELAM',ZGELAM,1)
     1307!        CALL writefield_phy('ZGEMU',ZGEMU,1)
     1308!        CALL writefield_phy('zmasq',zmasq,1)
    12621309        CALL writefield_phy('paprs_i',paprs_i,klev+1)
    1263         CALL writefield_phy('pplay_i',pplay_i,klev)
    1264         CALL writefield_phy('t_i',t_i,klev)
     1310!        CALL writefield_phy('pplay_i',pplay_i,klev)
     1311!        CALL writefield_phy('t_i',t_i,klev)
    12651312        CALL writefield_phy('ZTH_i',ZTH_i,klev+1)
    12661313        CALL writefield_phy('cldfra_i',cldfra_i,klev)
    1267         CALL writefield_phy('paer_i',PAER_i,klev)
    12681314        CALL writefield_phy('q_i',q_i,klev)
    12691315        CALL writefield_phy('fiwc_i',fiwc_i,klev)
     
    12711317        CALL writefield_phy('palbd_new',PALBD_NEW,NSW)
    12721318        CALL writefield_phy('palbp_new',PALBP_NEW,NSW)
    1273 !        CALL writefield_phy('ZO3_DP',ZO3_DP,klev)
     1319        CALL writefield_phy('POZON',POZON_i(:,:,1),klev)
     1320!        CALL writefield_phy('ZCO2',ZCO2,klev)
     1321!        CALL writefield_phy('ZCH4',ZCH4,klev)
     1322!        CALL writefield_phy('ZN2O',ZN2O,klev)
     1323!        CALL writefield_phy('ZO2',ZO2,klev)
     1324!        CALL writefield_phy('ZNO2',ZNO2,klev)
     1325!        CALL writefield_phy('ZCFC11',ZCFC11,klev)
     1326!        CALL writefield_phy('ZCFC12',ZCFC12,klev)
     1327!        CALL writefield_phy('ZHCFC22',ZHCFC22,klev)
     1328!        CALL writefield_phy('ZCCL4',ZCCL4,klev)
     1329        CALL writefield_phy('ref_liq_i',ref_liq_i,klev)
     1330        CALL writefield_phy('ref_ice_i',ref_ice_i,klev)
    12741331      ENDIF
    12751332 
     
    12781335! ??? naero_tot
    12791336      & day_cur, current_time, &
    1280 !      & solaire, &
    12811337      & PSCT, &
    12821338      & rmu0, tsol, PALBD_NEW,PALBP_NEW, &   
     
    12921348      & paprs_i, ZTH_i, q_i, qsat_i, &
    12931349!       Gas
    1294        & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, ZCCL4, POZON_i(:,:,1), &
     1350       & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, &
     1351       & ZCCL4, POZON_i(:,:,1), ZO2, &
    12951352!       nuages :
    12961353      & cldfra_i, flwc_i, fiwc_i, ZQ_RAIN, ZQ_SNOW, & 
     
    13161373
    13171374      IF (lldebug) THEN
    1318         CALL writefield_phy('zlwft_i',ZLWFT_i,klev+1)
    1319         CALL writefield_phy('zlwft0_ii',ZLWFT0_ii,klev+1)
    1320         CALL writefield_phy('zswft_i',ZSWFT_i,klev+1)
    1321         CALL writefield_phy('zswft0_i',ZSWFT0_ii,klev+1)
    1322         CALL writefield_phy('zfsdwn_i',ZFSDWN_i,klev+1)
    1323         CALL writefield_phy('zflux2_i',ZFLUX_i(:,2,:),klev+1)
    1324         CALL writefield_phy('zfcdwn_i',ZFCDWN_i,klev+1)
    1325         CALL writefield_phy('zfluc2_i',ZFLUC_i(:,2,:),klev+1)
    1326         CALL writefield_phy('psfswdir',PSFSWDIR,6)
    1327         CALL writefield_phy('psfswdif',PSFSWDIF,6)
    1328         CALL writefield_phy('zflux1_i',ZFLUX_i(:,1,:),klev+1)
    1329         CALL writefield_phy('zfluc1_i',ZFLUC_i(:,1,:),klev+1)
    1330         CALL writefield_phy('zfsup_i',ZFSUP_i,klev+1)
    1331         CALL writefield_phy('zfcup_i',ZFCUP_i,klev+1)
     1375        CALL writefield_phy('FLUX_LW',ZLWFT_i,klev+1)
     1376        CALL writefield_phy('FLUX_LW_CLEAR',ZLWFT0_ii,klev+1)
     1377        CALL writefield_phy('FLUX_SW',ZSWFT_i,klev+1)
     1378        CALL writefield_phy('FLUX_SW_CLEAR',ZSWFT0_ii,klev+1)
     1379        CALL writefield_phy('FLUX_DN_SW',ZFSDWN_i,klev+1)
     1380        CALL writefield_phy('FLUX_DN_LW',ZFLUX_i(:,2,:),klev+1)
     1381        CALL writefield_phy('FLUX_DN_SW_CLEAR',ZFCDWN_i,klev+1)
     1382        CALL writefield_phy('FLUX_DN_LW_CLEAR',ZFLUC_i(:,2,:),klev+1)
     1383        CALL writefield_phy('PSFSWDIR',PSFSWDIR,6)
     1384        CALL writefield_phy('PSFSWDIF',PSFSWDIF,6)
     1385        CALL writefield_phy('FLUX_UP_LW',ZFLUX_i(:,1,:),klev+1)
     1386        CALL writefield_phy('FLUX_UP_LW_CLEAR',ZFLUC_i(:,1,:),klev+1)
     1387        CALL writefield_phy('FLUX_UP_SW',ZFSUP_i,klev+1)
     1388        CALL writefield_phy('FLUX_UP_SW_CLEAR',ZFCUP_i,klev+1)
    13321389      ENDIF
    13331390! ---------
Note: See TracChangeset for help on using the changeset viewer.