7 hprogram, hcoupling, ptimec, ptstep, kyear, kmonth, kday, ptime, &
8 ki, ksv, ksw, ptsun, pzenith, pzenith2, pazim, pzref, puref, pu, pv, &
9 pqa, pta, prhoa, psv, pco2, hsv,prain, psnow, plw, pdir_sw, psca_sw, &
10 psw_bands, pps, ppa, psftq, psfth, psfts, psfco2, psfu, psfv, &
11 ptrad, pdir_alb, psca_alb, pemis, ptsurf, pz0, pz0h, pqsurf, &
12 ppew_a_coef, ppew_b_coef, ppet_a_coef, ppeq_a_coef, ppet_b_coef, &
53 USE modi_init_water_sbl
56 USE modi_canopy_grid_update
58 USE modi_coupling_seaflux_n
60 USE yomhook
,ONLY : lhook, dr_hook
61 USE parkind1
,ONLY : jprb
69 TYPE(dst_t),
INTENT(INOUT) :: dst
70 TYPE(slt_t),
INTENT(INOUT) :: slt
72 CHARACTER(LEN=6),
INTENT(IN) :: hprogram
73 CHARACTER(LEN=1),
INTENT(IN) :: hcoupling
76 REAL,
INTENT(IN) :: ptimec
77 INTEGER,
INTENT(IN) :: kyear
78 INTEGER,
INTENT(IN) :: kmonth
79 INTEGER,
INTENT(IN) :: kday
80 REAL,
INTENT(IN) :: ptime
81 INTEGER,
INTENT(IN) :: ki
82 INTEGER,
INTENT(IN) :: ksv
83 INTEGER,
INTENT(IN) :: ksw
84 REAL,
DIMENSION(KI),
INTENT(IN) :: ptsun
85 REAL,
INTENT(IN) :: ptstep
86 REAL,
DIMENSION(KI),
INTENT(IN) :: pzref
87 REAL,
DIMENSION(KI),
INTENT(IN) :: puref
89 REAL,
DIMENSION(KI),
INTENT(IN) :: pta
90 REAL,
DIMENSION(KI),
INTENT(IN) :: pqa
91 REAL,
DIMENSION(KI),
INTENT(IN) :: prhoa
92 REAL,
DIMENSION(KI,KSV),
INTENT(IN) :: psv
95 CHARACTER(LEN=6),
DIMENSION(KSV),
INTENT(IN):: hsv
96 REAL,
DIMENSION(KI),
INTENT(IN) :: pu
97 REAL,
DIMENSION(KI),
INTENT(IN) :: pv
98 REAL,
DIMENSION(KI,KSW),
INTENT(IN) :: pdir_sw
100 REAL,
DIMENSION(KI,KSW),
INTENT(IN) :: psca_sw
102 REAL,
DIMENSION(KSW),
INTENT(IN) :: psw_bands
103 REAL,
DIMENSION(KI),
INTENT(IN) :: pzenith
104 REAL,
DIMENSION(KI),
INTENT(IN) :: pzenith2
105 REAL,
DIMENSION(KI),
INTENT(IN) :: pazim
106 REAL,
DIMENSION(KI),
INTENT(IN) :: plw
108 REAL,
DIMENSION(KI),
INTENT(IN) :: pps
109 REAL,
DIMENSION(KI),
INTENT(IN) :: ppa
110 REAL,
DIMENSION(KI),
INTENT(IN) :: pco2
111 REAL,
DIMENSION(KI),
INTENT(IN) :: psnow
112 REAL,
DIMENSION(KI),
INTENT(IN) :: prain
115 REAL,
DIMENSION(KI),
INTENT(OUT) :: psfth
116 REAL,
DIMENSION(KI),
INTENT(OUT) :: psftq
117 REAL,
DIMENSION(KI),
INTENT(OUT) :: psfu
118 REAL,
DIMENSION(KI),
INTENT(OUT) :: psfv
119 REAL,
DIMENSION(KI),
INTENT(OUT) :: psfco2
120 REAL,
DIMENSION(KI,KSV),
INTENT(OUT):: psfts
122 REAL,
DIMENSION(KI),
INTENT(OUT) :: ptrad
123 REAL,
DIMENSION(KI,KSW),
INTENT(OUT):: pdir_alb
124 REAL,
DIMENSION(KI,KSW),
INTENT(OUT):: psca_alb
125 REAL,
DIMENSION(KI),
INTENT(OUT) :: pemis
127 REAL,
DIMENSION(KI),
INTENT(OUT) :: ptsurf
128 REAL,
DIMENSION(KI),
INTENT(OUT) :: pz0
129 REAL,
DIMENSION(KI),
INTENT(OUT) :: pz0h
130 REAL,
DIMENSION(KI),
INTENT(OUT) :: pqsurf
132 REAL,
DIMENSION(KI),
INTENT(IN) :: ppew_a_coef
133 REAL,
DIMENSION(KI),
INTENT(IN) :: ppew_b_coef
134 REAL,
DIMENSION(KI),
INTENT(IN) :: ppet_a_coef
135 REAL,
DIMENSION(KI),
INTENT(IN) :: ppeq_a_coef
136 REAL,
DIMENSION(KI),
INTENT(IN) :: ppet_b_coef
137 REAL,
DIMENSION(KI),
INTENT(IN) :: ppeq_b_coef
138 CHARACTER(LEN=2),
INTENT(IN) :: htest
144 REAL,
DIMENSION(KI) :: zwind
145 REAL,
DIMENSION(KI) :: zexna
146 REAL,
DIMENSION(KI) :: zta
147 REAL,
DIMENSION(KI) :: zpa
148 REAL,
DIMENSION(KI) :: zzref
149 REAL,
DIMENSION(KI) :: zuref
150 REAL,
DIMENSION(KI) :: zu
151 REAL,
DIMENSION(KI) :: zv
152 REAL,
DIMENSION(KI) :: zqa
153 REAL,
DIMENSION(KI) :: zpeq_a_coef
154 REAL,
DIMENSION(KI) :: zpeq_b_coef
159 REAL,
DIMENSION(KI) :: zsflux_u
160 REAL,
DIMENSION(KI) :: zsflux_t
161 REAL,
DIMENSION(KI) :: zsflux_q
162 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zforc_u
163 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zdforc_udu
165 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zforc_e
166 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zdforc_ede
168 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zforc_t
169 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zdforc_tdt
171 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zforc_q
172 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zdforc_qdq
174 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zlmo
175 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zlm
176 REAL,
DIMENSION(KI,SM%SSB%NLVL) :: zleps
177 REAL,
DIMENSION(KI) :: zh
178 REAL,
DIMENSION(KI) :: zustar
180 REAL,
DIMENSION(KI) :: zpet_a_coef
181 REAL,
DIMENSION(KI) :: zpet_b_coef
182 REAL,
DIMENSION(KI) :: zpew_a_coef
183 REAL,
DIMENSION(KI) :: zpew_b_coef
185 REAL,
DIMENSION(KI) :: zalfau
186 REAL,
DIMENSION(KI) :: zbetau
187 REAL,
DIMENSION(KI) :: zalfath
188 REAL,
DIMENSION(KI) :: zbetath
189 REAL,
DIMENSION(KI) :: zalfaq
190 REAL,
DIMENSION(KI) :: zbetaq
192 CHARACTER(LEN=1) :: gcoupling
193 REAL(KIND=JPRB) :: zhook_handle
200 IF (lhook) CALL dr_hook(
'COUPLING_SEAFLUX_SBL_N',0,zhook_handle)
210 sm%SSB%XZ,sm%SSB%XZF,sm%SSB%XDZ,sm%SSB%XDZF)
217 IF(any(sm%SSB%XT(:,:) == xundef))
THEN
218 CALL
init_water_sbl(sm%SSB%NLVL, ppa, pps, pta, pqa, prhoa, pu, pv, prain, psnow, &
219 psfth, psftq, pzref, puref, sm%S%XSST, sm%S%XZ0, sm%SSB%XZ, &
220 sm%SSB%XT, sm%SSB%XQ, sm%SSB%XU, sm%SSB%XTKE, sm%SSB%XP)
227 CALL
init_forc( zforc_u, zdforc_udu, zforc_e, zdforc_ede, &
228 zforc_t, zdforc_tdt, zforc_q, zdforc_qdq )
234 zlmo = spread(sm%SSB%XLMO,2,sm%SSB%NLVL)
239 zwind = sqrt(pu**2+pv**2)
240 CALL
canopy_evol(ki,sm%SSB%NLVL,ptstep,1,sm%SSB%XZ,zwind,pta,pqa,ppa,prhoa, &
241 zsflux_u,zsflux_t,zsflux_q, &
242 zforc_u,zdforc_udu,zforc_e,zdforc_ede, &
243 zforc_t,zdforc_tdt,zforc_q,zdforc_qdq, &
244 sm%SSB%XZ,sm%SSB%XZF,sm%SSB%XDZ,sm%SSB%XDZF, &
245 sm%SSB%XU,sm%SSB%XTKE,sm%SSB%XT, &
246 sm%SSB%XQ,zlmo,zlm,zleps,sm%SSB%XP,zustar, &
247 zalfau,zbetau,zalfath,zbetath,zalfaq,zbetaq )
256 sm%SSB%XQ(:,1), pu, pv, sm%SSB%XZ(:,1), &
258 prhoa, zalfau, zbetau, zalfath, &
259 zbetath, zalfaq, zbetaq, &
260 zpa, zta, zqa, zu, zv, &
261 zuref, zzref, zexna, &
262 zpew_a_coef, zpew_b_coef, &
263 zpet_a_coef, zpet_b_coef, &
264 zpeq_a_coef, zpeq_b_coef )
273 gcoupling = hcoupling
276 pps, ppa, pta, pqa, pu, pv, &
278 ppew_a_coef, ppew_b_coef, &
279 ppet_a_coef, ppet_b_coef, &
280 ppeq_a_coef, ppeq_b_coef, &
281 zpa, zta, zqa, zu, zv, &
283 zpew_a_coef, zpew_b_coef, &
284 zpet_a_coef, zpet_b_coef, &
285 zpeq_a_coef, zpeq_b_coef )
295 hprogram, gcoupling, ptimec, &
296 ptstep, kyear, kmonth, kday, ptime, &
298 ptsun, pzenith, pzenith2, pazim, &
299 zzref, zuref, zu, zv, zqa, zta, prhoa, psv, pco2, hsv, &
300 prain, psnow, plw, pdir_sw, psca_sw, psw_bands, pps, zpa, &
301 psftq, psfth, psfts, psfco2, psfu, psfv, &
302 ptrad, pdir_alb, psca_alb, pemis, ptsurf, pz0, pz0h, pqsurf, &
303 zpew_a_coef, zpew_b_coef, &
304 zpet_a_coef, zpeq_a_coef, zpet_b_coef, zpeq_b_coef, &
312 IF (.NOT. sm%S%LSBL .AND. lhook) CALL dr_hook(
'COUPLING_SEAFLUX_SBL_N',1,zhook_handle)
313 IF (.NOT. sm%S%LSBL)
RETURN
320 CALL
init_forc( zforc_u, zdforc_udu, zforc_e, zdforc_ede, &
321 zforc_t, zdforc_tdt, zforc_q, zdforc_qdq )
323 zsflux_u = - sqrt(psfu(:)**2+psfv(:)**2) / prhoa(:)
324 zsflux_t(:) = psfth(:) / xcpd * zexna(:) / prhoa(:)
325 zsflux_q(:) = psftq(:)
332 zwind = sqrt(pu**2+pv**2)
333 CALL
canopy_evol(ki,sm%SSB%NLVL,ptstep,2,sm%SSB%XZ,zwind,pta,pqa,ppa,prhoa, &
334 zsflux_u,zsflux_t,zsflux_q, &
335 zforc_u,zdforc_udu,zforc_e,zdforc_ede, &
336 zforc_t,zdforc_tdt,zforc_q,zdforc_qdq, &
337 sm%SSB%XZ,sm%SSB%XZF,sm%SSB%XDZ,sm%SSB%XDZF, &
338 sm%SSB%XU,sm%SSB%XTKE,sm%SSB%XT, &
339 sm%SSB%XQ,zlmo,zlm,zleps,sm%SSB%XP,zustar, &
340 zalfau,zbetau,zalfath,zbetath,zalfaq,zbetaq )
342 sm%SSB%XLMO(:) = zlmo(:,sm%SSB%NLVL)
350 IF (sm%S%LSBL .AND. sm%DGS%N2M>=1) CALL
init_2m_10m( sm%SSB%XP(:,2), sm%SSB%XT(:,2), &
351 sm%SSB%XQ(:,2), sm%SSB%XU, sm%SSB%XZ, &
352 pu, pv, zwind, prhoa, &
353 sm%DGS%XT2M, sm%DGS%XQ2M, sm%DGS%XHU2M, &
354 sm%DGS%XZON10M, sm%DGS%XMER10M,&
355 sm%DGS%XWIND10M, sm%DGS%XWIND10M_MAX, sm%DGS%XT2M_MIN, &
356 sm%DGS%XT2M_MAX, sm%DGS%XHU2M_MIN, sm%DGS%XHU2M_MAX )
358 IF (lhook) CALL dr_hook(
'COUPLING_SEAFLUX_SBL_N',1,zhook_handle)
subroutine init_coupling_canopy(PP, PPA, PT, PQ, PU, PV, PZ, PXU, PRHOA, PALFAU, PBETAU, PALFATH, PBETATH, PALFAQ, PBETAQ, PPPA, PTTA, PQQA, PUU, PVV, PUUREF, PZZREF, PEXNA, PPEW_AA_COEF, PPEW_BB_COEF, PPET_AA_COEF, PPET_BB_COEF, PPEQ_AA_COEF, PPEQ_BB_COEF)
subroutine coupling_seaflux_n(SM, DST, SLT, HPROGRAM, HCOUPLING, PTIMEC, PTSTEP, KYEAR, KMONTH, KDAY, PTIME, KI, KSV, KSW, PTSUN, PZENITH, PZENITH2, PAZIM, PZREF, PUREF, PU, PV, PQA, PTA, PRHOA, PSV, PCO2, HSV, PRAIN, PSNOW, PLW, PDIR_SW, PSCA_SW, PSW_BANDS, PPS, PPA, PSFTQ, PSFTH, PSFTS, PSFCO2, PSFU, PSFV, PTRAD, PDIR_ALB, PSCA_ALB, PEMIS, PTSURF, PZ0, PZ0H, PQSURF, PPEW_A_COEF, PPEW_B_COEF, PPET_A_COEF, PPEQ_A_COEF, PPET_B_COEF, PPEQ_B_COEF, HTEST)
subroutine canopy_grid_update(KI, KLVL, PH, PZFORC, PZ, PZF, PDZ, PDZF)
subroutine init_coupling(HCOUPLING, PPS, PPA, PTA, PQA, PU, PV, PUREF, PZREF, PPEW_A_COEF, PPEW_B_COEF, PPET_A_COEF, PPET_B_COEF, PPEQ_A_COEF, PPEQ_B_COEF, PPPA, PTTA, PQQA, PUU, PVV, PUUREF, PZZREF, PPEW_AA_COEF, PPEW_BB_COEF, PPET_AA_COEF, PPET_BB_COEF, PPEQ_AA_COEF, PPEQ_BB_COEF)
subroutine init_2m_10m(PP, PT, PQ, PXU, PXZ, PU, PV, PWIND, PRHOA, PT2M, PQ2M, PHU2M, PZON10M, PMER10M, PWIND10M, PWIND10M_MAX, PT2M_MIN, PT2M_MAX, PHU2M_MIN, PHU2M_MAX)
subroutine coupling_seaflux_sbl_n(SM, DST, SLT, HPROGRAM, HCOUPLING, PTIMEC, PTSTEP, KYEAR, KMONTH, KDAY, PTIME, KI, KSV, KSW, PTSUN, PZENITH, PZENITH2, PAZIM, PZREF, PUREF, PU, PV, PQA, PTA, PRHOA, PSV, PCO2, HSV, PRAIN, PSNOW, PLW, PDIR_SW, PSCA_SW, PSW_BANDS, PPS, PPA, PSFTQ, PSFTH, PSFTS, PSFCO2, PSFU, PSFV, PTRAD, PDIR_ALB, PSCA_ALB, PEMIS, PTSURF, PZ0, PZ0H, PQSURF, PPEW_A_COEF, PPEW_B_COEF, PPET_A_COEF, PPEQ_A_COEF, PPET_B_COEF, PPEQ_B_COEF, HTEST)
subroutine init_forc(PFORC_U, PDFORC_UDU, PFORC_E, PDFORC_EDE, PFORC_T, PDFORC_TDT, PFORC_Q, PDFORC_QDQ)
subroutine init_water_sbl(KLVL, PPA, PPS, PTA, PQA, PRHOA, PU, PV, PRAIN, PSNOW, PSFTH, PSFTQ, PZREF, PUREF, PTS, PZ0, PZ, PT, PQ, PWIND, PTKE, PP)
subroutine canopy_evol(KI, KLVL, PTSTEP, KIMPL, PZZ, PWIND, PTA, PQA, PPA, PRHOA, PSFLUX_U, PSFLUX_T, PSFLUX_Q, PFORC_U, PDFORC_UDU, PFORC_E, PDFORC_EDE, PFORC_T, PDFORC_TDT, PFORC_Q, PDFORC_QDQ, PZ, PZF, PDZ, PDZF, PU, PTKE, PT, PQ, PLMO, PLM, PLEPS, PP, PUSTAR, PALFAU, PBETAU, PALFATH, PBETATH, PALFAQ, PBETAQ, ONEUTRAL)