Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynhpg_tam.F90
- Timestamp:
- 2011-02-15T12:58:59+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynhpg_tam.F90
r1885 r2587 109 109 PUBLIC dyn_hpg_adj ! routine called by step_tam module 110 110 PUBLIC dyn_hpg_adj_tst! routine called by test module 111 #if defined key_tst_tlm 111 112 PUBLIC dyn_hpg_tlm_tst! routine called by test module 113 #endif 112 114 113 115 !!* Namelist nam_dynhpg : Choice of horizontal pressure gradient computation … … 1059 1061 & ) 1060 1062 END SUBROUTINE dyn_hpg_adj_tst 1061 1063 #if defined key_tst_tlm 1062 1064 SUBROUTINE dyn_hpg_tlm_tst( kumadt ) 1063 1065 !!----------------------------------------------------------------------- … … 1089 1091 USE tamtrj ! writing out state trajectory 1090 1092 USE par_tlm, ONLY: & 1093 & tlm_bch, & 1091 1094 & cur_loop, & 1092 1095 & h_ratio … … 1160 1163 & jk 1161 1164 CHARACTER(LEN=14) :: cl_name 1162 CHARACTER (LEN=128) :: file_out, file_wop 1165 CHARACTER (LEN=128) :: file_out, file_wop, file_xdx 1163 1166 CHARACTER (LEN=90) :: FMT 1164 1167 REAL(KIND=wp), DIMENSION(100):: & … … 1228 1231 ! Output filename Xn=F(X0) 1229 1232 !-------------------------------------------------------------------- 1230 file_wop='trj_wop_dynhpg'1231 1232 1233 CALL tlm_namrd 1233 1234 gamma = h_ratio 1234 1235 file_wop='trj_wop_dynhpg' 1236 file_xdx='trj_xdx_dynhpg' 1235 1237 !-------------------------------------------------------------------- 1236 1238 ! Initialize the tangent input with random noise: dx … … 1262 1264 ! Complete Init for Direct 1263 1265 !------------------------------------------------------------------- 1264 CALL istate_p1266 IF ( tlm_bch /= 2 ) CALL istate_p 1265 1267 1266 1268 ! *** initialize the reference trajectory … … 1296 1298 ! Compute the direct model F(X0,t=n) = Xn 1297 1299 !-------------------------------------------------------------------- 1298 CALL dyn_hpg(nit000) 1299 IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop) 1300 IF ( tlm_bch /= 2 ) CALL dyn_hpg(nit000) 1301 IF ( tlm_bch == 0 ) CALL trj_wri_spl(file_wop) 1302 IF ( tlm_bch == 1 ) CALL trj_wri_spl(file_xdx) 1300 1303 !-------------------------------------------------------------------- 1301 1304 ! Compute the Tangent 1302 1305 !-------------------------------------------------------------------- 1303 IF ( cur_loop .NE. 0) THEN 1304 !-------------------------------------------------------------------- 1305 ! Storing data 1306 !-------------------------------------------------------------------- 1307 zua_out (:,:,:) = ua (:,:,:) 1308 zva_out (:,:,:) = va (:,:,:) 1309 1306 IF ( tlm_bch == 2 ) THEN 1310 1307 !-------------------------------------------------------------------- 1311 1308 ! Initialize the tangent variables … … 1324 1321 ! Compute the scalar product: ( L(t0,tn) gamma dx0 ) ) 1325 1322 !-------------------------------------------------------------------- 1326 1327 1323 zsp2_1 = DOT_PRODUCT( ua_tl, ua_tl ) 1328 1324 zsp2_2 = DOT_PRODUCT( va_tl, va_tl ) … … 1334 1330 zua_wop (:,:,:) = ua (:,:,:) 1335 1331 zva_wop (:,:,:) = va (:,:,:) 1332 CALL trj_rd_spl(file_xdx) 1333 zua_out (:,:,:) = ua (:,:,:) 1334 zva_out (:,:,:) = va (:,:,:) 1336 1335 !-------------------------------------------------------------------- 1337 1336 ! Compute the Linearization Error … … 1474 1473 !!====================================================================== 1475 1474 #endif 1475 #endif 1476 1476 END MODULE dynhpg_tam
Note: See TracChangeset
for help on using the changeset viewer.