- Timestamp:
- 2015-02-23T17:07:54+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5094_UKMO_ISFCLEAN/NEMOGCM/NEMO/OPA_SRC/step.F90
r5098 r5104 146 146 IF( lk_ldfslp ) THEN ! slope of lateral mixing 147 147 CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density 148 IF( ln_zps ) CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient 148 IF( ln_zps .AND. .NOT. ln_isfcav) & 149 & CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient 149 150 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 150 IF( ln_zps .AND. ln_isfcav)&151 IF( ln_zps .AND. ln_isfcav) & 151 152 & CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps for top cell (ISF) 152 153 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & … … 180 181 ! is necessary to compute momentum advection for the rhs of barotropic loop: 181 182 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 182 IF( ln_zps ) CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient 183 IF( ln_zps .AND. .NOT. ln_isfcav) & 184 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient 183 185 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 184 IF( ln_zps .AND. ln_isfcav)&186 IF( ln_zps .AND. ln_isfcav) & 185 187 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps for top cell (ISF) 186 188 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & … … 259 261 CALL tra_nxt( kstp ) ! tracer fields at next time step 260 262 CALL eos ( tsa, rhd, rhop, fsdept_n(:,:,:) ) ! Time-filtered in situ density for hpg computation 261 IF( ln_zps ) CALL zps_hde ( kstp, jpts, tsa, gtsu, gtsv, & ! Partial steps: before horizontal gradient 263 IF( ln_zps .AND. .NOT. ln_isfcav) & 264 & CALL zps_hde ( kstp, jpts, tsa, gtsu, gtsv, & ! Partial steps: before horizontal gradient 262 265 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 263 IF( ln_zps .AND. ln_isfcav)&266 IF( ln_zps .AND. ln_isfcav) & 264 267 & CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv, & ! Partial steps for top cell (ISF) 265 268 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & … … 268 271 IF ( .NOT. lk_dynspg_ts ) THEN ! eos already called in time-split case 269 272 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 270 IF( ln_zps ) CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient 271 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 272 IF( ln_zps .AND. ln_isfcav) & 273 IF( ln_zps .AND. .NOT. ln_isfcav) & 274 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient 275 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 276 IF( ln_zps .AND. ln_isfcav) & 273 277 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps for top cell (ISF) 274 278 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , &
Note: See TracChangeset
for help on using the changeset viewer.