- Timestamp:
- 2015-07-10T13:28:53+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4765_CNRS_agrif/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90
r4624 r5581 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 537 CALL eos( pts, rhd, rhop, gdept_0(:,:,:) ) ! Time-filtered in situ density 538 CALL bn2( pts, rn2 ) ! before Brunt-Vaisala frequency 539 IF( ln_zps ) & 540 & CALL zps_hde( kt, jpts, pts, gtsu, gtsv, rhd, gru, grv ) ! Partial steps: before Horizontal DErivative 541 ! ! of t, s, rd at the bottom ocean level 542 CALL eos ( pts, rhd, rhop, gdept_0(:,:,:) ) 543 CALL eos_rab( pts, rab_n ) ! now local thermal/haline expension ratio at T-points 544 CALL bn2 ( pts, rab_n, rn2 ) ! now Brunt-Vaisala 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 542 556 CALL zdf_mxl( kt ) ! mixed layer depth 543 557 CALL ldf_slp( kt, rhd, rn2 ) ! slopes
Note: See TracChangeset
for help on using the changeset viewer.