c CALL ioconf_calendar('360d') annee_ref=1988 day_ref= 1 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) ! CALL ioget_calendar(un_an) ! zstophy = un_an/12. ! zout = un_an/12. ! zstophy = 1036800. ! zout = 1036800. zstophy = dtime zout = dtime print*,'annee_ref day_ref zjulian ',annee_ref,day_ref,zjulian ! print*,'un_an itau_phy dtime nhori nid_day ',un_an,itau_phy, ! s dtime,nhori print*,'zstophy,zout ',zstophy,zout c !- itau_phy : time step at which the history tape starts !- zjulian : The Julian date at which the itau was equal to 0 !- dtime : Time step in seconds. Time step of the counter itau !- used in histwrt for instance CALL histbeg("histforcing", iim,rlon_1d, jjm+1,rlat_1d, . 1,iim,1,jjm+1,itau_phy,zjulian,dtime, . nhori, nid_day) CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", . llm, presl_1d/100., nvert) c c Flux nets radiatifs de reference a la tropopause !- zstophy: Frequency of this operation (in seconds) !- zout: Frequency at which the variable should be !- written (in seconds) ! CALL histdef(nid_day, "aire", "Temp sol", "K", ! . iim,jjm+1,nhori, 1,1,1, -99, 32, ! . "once", zstophy,zout) c CALL histdef(nid_day, "albs", "Surf. albedo", "-", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c c Flux nets radiatifs de reference c toa CALL histdef(nid_day, "toa_sw_ref", "SW flux at TOA ", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "toa_sw_ref0", "SW flux at TOA", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "toa_lw_ref","LW flux at toa", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "toa_lw_ref0", "LW flux at toa", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c tps CALL histdef(nid_day, "tps_sw_ref", "SW flux at TPS", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "tps_sw_ref0", "SW flux at TPS", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "tps_lw_ref","LW flux at tps", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "tps_lw_ref0", "LW flux at tps", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c srf CALL histdef(nid_day, "srf_sw_ref", "SW flux at SRF ", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "srf_sw_ref0", "SW flux at SRF", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "srf_lw_ref","LW flux at srf", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "srf_lw_ref0", "LW flux at srf", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) ! ! Flux initiaux après perturbation, sans ajustement strato ! c toa CALL histdef(nid_day, "toa_sw_ini","LW flux at toa", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "toa_sw_ini0","LW flux at toa", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "toa_lw_ini", "LW flux at toa", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "toa_lw_ini0", "LW flux at toa", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c tps c CALL histdef(nid_day, "tps_sw_ini", "SW flux at TPS", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "tps_sw_ini0", "SW flux at TPS", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "tps_lw_ini", "LW flux at TPS", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "tps_lw_ini0", "LW flux at TPS", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) ! ! Forcages initiaux après perturbation, sans ajustement strato ! c toa CALL histdef(nid_day, "d_toa_sw_ini", "I SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_toa_sw_ini0", "I SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_toa_lw_ini", "I SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_toa_lw_ini0", "I SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c tps CALL histdef(nid_day,"d_tps_sw_ini","I SW rad forc", . "w/m2",iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day,"d_tps_sw_ini0","I SW rad forc", . "w/m2",iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day,"d_tps_lw_ini","I LW rad forc","w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day,"d_tps_lw_ini0","I LW rad forc","w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) ! ! Forcages après ajustement strato ! c toa CALL histdef(nid_day, "d_toa_sw_adj", "Adj SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_toa_sw_adj0", "I SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_toa_lw_adj", "Adj LW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_toa_lw_adj0", "I LW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c tps CALL histdef(nid_day,"d_tps_sw_adj","Adj SW forcing","w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day,"d_tps_sw_adj0","Adj SW forcing","w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_tps_lw_adj", "Adj LW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_tps_lw_adj0", "I LW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c srf CALL histdef(nid_day, "d_srf_sw_adj", "Adj SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_srf_sw_adj0", "I SW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_srf_lw_adj", "Adj LW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day, "d_srf_lw_adj0", "I LW forcing", "w/m2", . iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) ! ! Bilan radiatif de la stratosphere ! CALL histdef(nid_day,"bilq_ref","reference radiatiave inbalance", . "w/m2",iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day,"bilq_ini","initial radiatiave inbalance", . "w/m2",iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) CALL histdef(nid_day,"bilq_adj","adjusted radiatiave inbalance", . "w/m2",iim,jjm+1,nhori, 1,1,1, -99, 32, . "inst(X)", zstophy,zout) c c Champs 3D: c c CALL histdef(nid_day, "temp", "Air temperature", "K", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "temp_adj", "Adj. Air temp.", "K", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "dHrad_dT", "Rad. heating rate deriv", "K/s/K", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "heat_ref", "Reference heating rate", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "heat0_ref", "Reference heating rate CS", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "cool_ref", "Reference cooling rate", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "cool0_ref", "Reference cooling rate CS", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_heat_ini", "Heating rate ini diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_heat0_ini", "Heating rate CS ini diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_cool_ini", "Cooling rate ini diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_cool0_ini", "Cooling rate CS ini diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_heat_adj", "Heating rate adj diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_heat0_adj", "Heating rate CS adj diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_cool_adj", "Cooling rate adj diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) c CALL histdef(nid_day, $ "d_cool0_adj", "Cooling rate CS adj diff", "K/s", . iim,jj_nb,nhori, klev,1,klev,nvert, 32, . "inst(X)", zstophy,zout) CALL histend(nid_day)