- Timestamp:
- 2012-11-27T15:42:24+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90
r3625 r3680 61 61 INTEGER , SAVE :: jf_mld ! index of mixed layer deptht 62 62 INTEGER , SAVE :: jf_emp ! index of water flux 63 INTEGER , SAVE :: jf_emps ! index of water flux - concentr/dilution 63 64 INTEGER , SAVE :: jf_qsr ! index of solar radiation 64 65 INTEGER , SAVE :: jf_wnd ! index of wind speed … … 242 243 ENDIF 243 244 ! 244 tsn(:,:,:,jp_tem) = sf_dyn(jf_tem)%fnow(:,:,:) * tmask(:,:,:) ! temperature245 tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:) * tmask(:,:,:) ! salinity245 tsn(:,:,:,jp_tem) = sf_dyn(jf_tem)%fnow(:,:,:) * tmask(:,:,:) ! temperature 246 tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:) * tmask(:,:,:) ! salinity 246 247 ! 247 248 CALL eos ( tsn, rhd, rhop ) ! In any case, we need rhop 248 249 CALL zdf_mxl( kt ) ! In any case, we need mxl 249 250 ! 250 avt(:,:,:) = sf_dyn(jf_avt)%fnow(:,:,:) * tmask(:,:,:) ! vertical diffusive coefficient251 un (:,:,:) = sf_dyn(jf_uwd)%fnow(:,:,:) * umask(:,:,:) ! u-velocity252 vn (:,:,:) = sf_dyn(jf_vwd)%fnow(:,:,:) * vmask(:,:,:) ! v-velocity251 avt(:,:,:) = sf_dyn(jf_avt)%fnow(:,:,:) * tmask(:,:,:) ! vertical diffusive coefficient 252 un (:,:,:) = sf_dyn(jf_uwd)%fnow(:,:,:) * umask(:,:,:) ! u-velocity 253 vn (:,:,:) = sf_dyn(jf_vwd)%fnow(:,:,:) * vmask(:,:,:) ! v-velocity 253 254 IF( .NOT.ln_dynwzv ) & ! w-velocity read in file 254 255 wn (:,:,:) = sf_dyn(jf_wwd)%fnow(:,:,:) * tmask(:,:,:) … … 305 306 CALL prt_ctl(tab2d_1=hmld , clinfo1=' hmld - : ', mask1=tmask, ovlap=1 ) 306 307 CALL prt_ctl(tab2d_1=sfx , clinfo1=' sfx - : ', mask1=tmask, ovlap=1 ) 308 CALL prt_ctl(tab2d_1=emp , clinfo1=' emp - : ', mask1=tmask, ovlap=1 ) 307 309 CALL prt_ctl(tab2d_1=wndm , clinfo1=' wspd - : ', mask1=tmask, ovlap=1 ) 308 310 CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask, ovlap=1 ) … … 349 351 sn_sal = FLD_N( 'dyna_grid_T' , 120 , 'vosaline' , .true. , .true. , 'yearly' , '' , '' ) 350 352 sn_mld = FLD_N( 'dyna_grid_T' , 120 , 'somixght' , .true. , .true. , 'yearly' , '' , '' ) 351 sn_emp = FLD_N( 'dyna_grid_T' , 120 , 'sowaflcd' , .true. , .true. , 'yearly' , '' , '' ) 353 sn_emp = FLD_N( 'dyna_grid_T' , 120 , 'sowaflup' , .true. , .true. , 'yearly' , '' , '' ) 354 sn_emps = FLD_N( 'dyna_grid_T' , 120 , 'sowaflcd' , .true. , .true. , 'yearly' , '' , '' ) 352 355 !! sn_emp = FLD_N( 'dyna_grid_T' , 120 , 'sowaflup' , .true. , .true. , 'yearly' , '' , '' ) ! v3.5+ 353 356 sn_sfx = FLD_N( 'dyna_grid_T' , 120 , 'sosfldow' , .true. , .true. , 'yearly' , '' , '' ) ! v3.5+ … … 392 395 ENDIF 393 396 394 jf_tem = 1 ; jf_sal = 2 ; jf_mld = 3 ; jf_emp = 4 ; jf_ ice = 5 ; jf_qsr = 6395 jf_wnd = 7 ; jf_uwd = 8 ; jf_vwd = 9 ; jf_wwd = 10 ; jf_avt = 11 ; jfld = 11396 ! 397 slf_d(jf_tem) = sn_tem ; slf_d(jf_sal) = sn_sal ; slf_d(jf_mld) = sn_mld398 slf_d(jf_emp) = sn_emp ; slf_d(jf_ ice) = sn_ice ; slf_d(jf_qsr) = sn_qsr399 slf_d(jf_ wnd) = sn_wnd ; slf_d(jf_uwd) = sn_uwd ; slf_d(jf_vwd) = sn_vwd400 slf_d(jf_ wwd) = sn_wwd ; slf_d(jf_avt) = sn_avt397 jf_tem = 1 ; jf_sal = 2 ; jf_mld = 3 ; jf_emp = 4 ; jf_emps = 5 ; jf_ice = 6 ; jf_qsr = 7 398 jf_wnd = 8 ; jf_uwd = 9 ; jf_vwd = 10 ; jf_wwd = 11 ; jf_avt = 12 ; jfld = 12 399 ! 400 slf_d(jf_tem) = sn_tem ; slf_d(jf_sal) = sn_sal ; slf_d(jf_mld) = sn_mld 401 slf_d(jf_emp) = sn_emp ; slf_d(jf_emps) = sn_emps ; slf_d(jf_ice) = sn_ice 402 slf_d(jf_qsr) = sn_qsr ; slf_d(jf_wnd) = sn_wnd ; slf_d(jf_avt) = sn_avt 403 slf_d(jf_uwd) = sn_uwd ; slf_d(jf_vwd) = sn_vwd ; slf_d(jf_wwd) = sn_wwd 401 404 ! 402 405 IF( .NOT.ln_degrad ) THEN ! no degrad option 403 406 IF( lk_traldf_eiv .AND. ln_dynbbl ) THEN ! eiv & bbl 404 jf_ubl = 1 2 ; jf_vbl = 13 ; jf_eiw = 14 ; jfld = 14407 jf_ubl = 13 ; jf_vbl = 14 ; jf_eiw = 15 ; jfld = 15 405 408 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl ; slf_d(jf_eiw) = sn_eiw 406 409 ENDIF 407 410 IF( .NOT.lk_traldf_eiv .AND. ln_dynbbl ) THEN ! no eiv & bbl 408 jf_ubl = 1 2 ; jf_vbl = 13 ; jfld = 13411 jf_ubl = 13 ; jf_vbl = 14 ; jfld = 14 409 412 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl 410 413 ENDIF 411 414 IF( lk_traldf_eiv .AND. .NOT.ln_dynbbl ) THEN ! eiv & no bbl 412 jf_eiw = 1 2 ; jfld = 12; slf_d(jf_eiw) = sn_eiw415 jf_eiw = 13 ; jfld = 13 ; slf_d(jf_eiw) = sn_eiw 413 416 ENDIF 414 417 ELSE 415 jf_ahu = 1 2 ; jf_ahv = 13 ; jf_ahw = 14 ; jfld = 14418 jf_ahu = 13 ; jf_ahv = 14 ; jf_ahw = 15 ; jfld = 15 416 419 slf_d(jf_ahu) = sn_ahu ; slf_d(jf_ahv) = sn_ahv ; slf_d(jf_ahw) = sn_ahw 417 420 IF( lk_traldf_eiv .AND. ln_dynbbl ) THEN ! eiv & bbl 418 jf_ubl = 1 5 ; jf_vbl = 16421 jf_ubl = 16 ; jf_vbl = 17 419 422 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl 420 jf_eiu = 1 7 ; jf_eiv = 18 ; jf_eiw = 19 ; jfld = 19423 jf_eiu = 18 ; jf_eiv = 19 ; jf_eiw = 20 ; jfld = 20 421 424 slf_d(jf_eiu) = sn_eiu ; slf_d(jf_eiv) = sn_eiv ; slf_d(jf_eiw) = sn_eiw 422 425 ENDIF 423 426 IF( .NOT.lk_traldf_eiv .AND. ln_dynbbl ) THEN ! no eiv & bbl 424 jf_ubl = 1 5 ; jf_vbl = 16 ; jfld = 16427 jf_ubl = 16 ; jf_vbl = 17 ; jfld = 17 425 428 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl 426 429 ENDIF 427 430 IF( lk_traldf_eiv .AND. .NOT.ln_dynbbl ) THEN ! eiv & no bbl 428 jf_eiu = 1 5 ; jf_eiv = 16 ; jf_eiw = 17 ; jfld = 17431 jf_eiu = 16 ; jf_eiv = 17 ; jf_eiw = 18 ; jfld = 18 429 432 slf_d(jf_eiu) = sn_eiu ; slf_d(jf_eiv) = sn_eiv ; slf_d(jf_eiw) = sn_eiw 430 433 ENDIF … … 440 443 ! Open file for each variable to get his number of dimension 441 444 DO ifpr = 1, jfld 442 CALL iom_open( slf_d(ifpr)%clname, inum )445 CALL iom_open( TRIM( cn_dir )//TRIM( slf_d(ifpr)%clname ), inum ) 443 446 idv = iom_varid( inum , slf_d(ifpr)%clvar ) ! id of the variable sdjf%clvar 444 447 idimv = iom_file ( inum )%ndims(idv) ! number of dimension for variable sdjf%clvar
Note: See TracChangeset
for help on using the changeset viewer.