7 hprogram,ki,pt2m_o,htest)
41 USE yomhook
, ONLY : lhook, dr_hook
42 USE parkind1
, ONLY : jprb
52 TYPE(teb_t),
INTENT(INOUT) :: t
55 CHARACTER(LEN=6),
INTENT(IN) :: hprogram
56 INTEGER,
INTENT(IN) :: ki
57 REAL,
DIMENSION(KI),
INTENT(IN) :: pt2m_o
58 CHARACTER(LEN=2),
INTENT(IN) :: htest
64 REAL,
DIMENSION (KI) :: ztrd3
65 REAL,
DIMENSION (KI) :: zt2inc
66 REAL,
DIMENSION (KI) :: ztcls
68 REAL(KIND=JPRB) :: zhook_handle
70 IF (lhook) CALL dr_hook(
'ASSIM_TEB_N',0,zhook_handle)
73 CALL
abor1_sfx(
'ASSIM_TEB_n: FATAL ERROR DURING ARGUMENT TRANSFER')
76 WRITE(*,*)
'UPDATING TOWN FOR SCHEME: ',trim(u%CTOWN)
78 IF ( top%NROAD_LAYER < 3 ) CALL
abor1_sfx(
'ASSIM_TEB_n: Only imlemented with 3 or more layers')
80 ztrd3(:) = t%CUR%XT_ROAD(:,3)
81 ztcls(:) = xat2m_teb(:)
86 WHERE ( pt2m_o(:) /= 999. )
87 zt2inc(:) = pt2m_o(:) - ztcls(:)
90 IF ( nprintlev > 0 .AND. ki>0)
WRITE(*,*)
'Mean T2m increments over TOWN ',sum(zt2inc)/ki
95 WHERE (ztrd3(:)/=xundef)
96 ztrd3(:) = ztrd3(:) + zt2inc(:)/(2.0*xpi)
100 IF (ki>0)
WRITE(*,*)
'Mean T_ROAD3 increments over TOWN ',sum(zt2inc)/ki
103 t%CUR%XT_ROAD(:,3) = ztrd3
105 IF (lhook) CALL dr_hook(
'ASSIM_TEB_N',1,zhook_handle)
subroutine assim_teb_n(U, T, TOP, HPROGRAM, KI, PT2M_O, HTEST)
subroutine abor1_sfx(YTEXT)