diff --git a/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/cos_sza.F90 b/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/cos_sza.F90 index f02aeb1dc0d8b66f526824678b7bcea1428cb5f0..215743cfb0e67a32125a895a48f7114e5099eba5 100644 --- a/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/cos_sza.F90 +++ b/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/cos_sza.F90 @@ -135,7 +135,7 @@ DO JCOL = KSTART,KEND ENDDO ! Computation of solar hour angle from sunposn -ZTIME = TDTRAD_FULL%TIME + 0.5*XDTRAD +ZTIME = TDTRAD_FULL%XTIME + 0.5*XDTRAD ZUT = MOD( 24.0+MOD(ZTIME/3600.,24.0),24.0 ) ZTUT = ZUT - XTSIDER ZSOLANG = (ZTUT-12.0)*15.0*(RPI/180.) ! hour angle in radians diff --git a/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/radiation_scheme.F90 b/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/radiation_scheme.F90 index 4f00fba6f7b381a844ba27db13d6042dc9668d28..e0cdd612165c9401974af01bd5de8a1092087509 100644 --- a/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/radiation_scheme.F90 +++ b/src/LIB/RAD/ecrad-1.0.1_mnh/ifs/radiation_scheme.F90 @@ -79,7 +79,7 @@ USE YOMHOOK , ONLY : LHOOK, DR_HOOK USE MODD_PARAM_ECRAD_n, ONLY : NAERMACC, NDECOLAT, XCLOUD_FRAC_STD, & ! replace YRERAD to get attributes & LAPPROXLWUPDATE, LAPPROXSWUPDATE USE MODD_RADIATIONS_n , ONLY : NSWB_MNH,NSWB_OLD -USE MODE_THERMO ! , ONLY QSATW_2D +USE MODE_THERMO , ONLY : QSAT USE MODD_DYN_n , ONLY : XTSTEP, NSTOP USE MODD_TIME , ONLY : TDTEXP USE MODD_TIME_n , ONLY : TDTMOD,TDTCUR @@ -336,7 +336,7 @@ thermodynamics%temperature_hl(KIDIA:KFDIA,KLEV+1) & ! & PPRESSURE, PTEMPERATURE, thermodynamics%h2o_sat_liq, 2) !MNH -thermodynamics%h2o_sat_liq(:,:) = QSATW_2D(REAL(PPRESSURE), REAL(PTEMPERATURE)) +thermodynamics%h2o_sat_liq(:,:) = QSAT(REAL(PPRESSURE), REAL(PTEMPERATURE)) thermodynamics%h2o_sat_liq(:,:) = thermodynamics%h2o_sat_liq(:,:) & & / (1.+thermodynamics%h2o_sat_liq(:,:)) ! mixing ratio => spec humid !MNH @@ -361,8 +361,8 @@ single_level%lw_emissivity(KIDIA:KFDIA,2) = PEMIS_WINDOW(KIDIA:KFDIA) ! Create the relevant seed from date and time get the starting day ! and number of minutes since start ! IDAY = NDD(NINDAT) ! NINDAT is AAAAMMDD initial date NDD extract DD as INTEGER -IDAY = TDTEXP%TDATE%DAY ! MNH -ITIM = NINT((TDTMOD%time-TDTCUR%time) / 60.0_JPRB) ! YRRIP contains timestep infos ; number of minutes since beginning +IDAY = TDTEXP%DATE%NDAY ! MNH +ITIM = NINT((TDTMOD%xtime-TDTCUR%xtime) / 60.0_JPRB) ! YRRIP contains timestep infos ; number of minutes since beginning DO JLON = KIDIA, KFDIA ! This method gives a unique value for roughly every 1-km square ! on the globe and every minute. ASIN(PGEMU)*60 gives rough