- Timestamp:
- 2015-06-19T17:52:53+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5021_nn_etau_revision/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90
r5239 r5446 245 245 tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:) * tmask(:,:,:) ! salinity 246 246 ! 247 CALL eos ( tsn, rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop 247 ! 248 CALL eos ( tsn, rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop 249 CALL eos_rab( tsn, rab_n ) ! now local thermal/haline expension ratio at T-points 250 CALL bn2 ( tsn, rab_n, rn2 ) ! before Brunt-Vaisala frequency need for zdfmxl 251 252 rn2b(:,:,:) = rn2(:,:,:) ! need for zdfmxl 248 253 CALL zdf_mxl( kt ) ! In any case, we need mxl 249 254 ! … … 259 264 fr_i(:,:) = sf_dyn(jf_ice)%fnow(:,:,1) * tmask(:,:,1) ! Sea-ice fraction 260 265 qsr (:,:) = sf_dyn(jf_qsr)%fnow(:,:,1) * tmask(:,:,1) ! solar radiation 261 IF 266 IF( ln_dynrnf ) & 262 267 rnf (:,:) = sf_dyn(jf_rnf)%fnow(:,:,1) * tmask(:,:,1) ! river runoffs 263 268 … … 383 388 384 389 ! 385 IF 390 IF( ln_dynrnf ) THEN 386 391 jf_rnf = jfld + 1 ; jfld = jf_rnf 387 392 slf_d(jf_rnf) = sn_rnf … … 535 540 !!--------------------------------------------------------------------- 536 541 #if defined key_ldfslp && ! defined key_c1d 542 CALL eos ( pts, rhd, rhop, gdept_0(:,:,:) ) 537 543 CALL eos_rab( pts, rab_n ) ! now local thermal/haline expension ratio at T-points 538 544 CALL bn2 ( pts, rab_n, rn2 ) ! now Brunt-Vaisala 539 IF( ln_zps ) & ! Partial steps: before Horizontal DErivative 540 & CALL zps_hde( kt, jpts, pts, gtsu, gtsv, & ! Partial steps: before horizontal gradient 541 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & ! 542 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level 543 ! only gtsu, gtsv, rhd, gru , grv are used 544 545 546 ! ! of t, s, rd at the bottom ocean level 545 546 ! Partial steps: before Horizontal DErivative 547 IF( ln_zps .AND. .NOT. ln_isfcav) & 548 & CALL zps_hde ( kt, jpts, pts, gtsu, gtsv, & ! Partial steps: before horizontal gradient 549 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 550 IF( ln_zps .AND. ln_isfcav) & 551 & CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv, & ! Partial steps for top cell (ISF) 552 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 553 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the first ocean level 554 555 rn2b(:,:,:) = rn2(:,:,:) ! need for zdfmxl 547 556 CALL zdf_mxl( kt ) ! mixed layer depth 548 557 CALL ldf_slp( kt, rhd, rn2 ) ! slopes
Note: See TracChangeset
for help on using the changeset viewer.