- Timestamp:
- 2021-06-17T14:14:52+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/ticket2680_C1D_PAPA/src/OFF/dtadyn.F90
r14310 r15010 21 21 !!---------------------------------------------------------------------- 22 22 USE oce ! ocean dynamics and tracers variables 23 USE c1d ! 1D configuration: l k_c1d23 USE c1d ! 1D configuration: ln_c1d 24 24 USE dom_oce ! ocean domain: variables 25 25 #if defined key_qco … … 135 135 CALL fld_read( kt, 1, sf_dyn ) != read data at kt time step ==! 136 136 ! 137 IF( l_ldfslp .AND. .NOT.l k_c1d ) CALL dta_dyn_slp( kt, Kbb, Kmm ) ! Computation of slopes137 IF( l_ldfslp .AND. .NOT.ln_c1d ) CALL dta_dyn_slp( kt, Kbb, Kmm ) ! Computation of slopes 138 138 ! 139 139 ts(:,:,:,jp_tem,Kmm) = sf_dyn(jf_tem)%fnow(:,:,:) * tmask(:,:,:) ! temperature … … 187 187 avs(:,:,:) = avt(:,:,:) 188 188 ! 189 IF( ln_trabbl .AND. .NOT.l k_c1d ) THEN ! diffusive Bottom boundary layer param189 IF( ln_trabbl .AND. .NOT.ln_c1d ) THEN ! diffusive Bottom boundary layer param 190 190 ahu_bbl(:,:) = sf_dyn(jf_ubl)%fnow(:,:,1) * umask(:,:,1) ! bbl diffusive coef 191 191 ahv_bbl(:,:) = sf_dyn(jf_vbl)%fnow(:,:,1) * vmask(:,:,1) … … 321 321 END DO 322 322 ! 323 IF( l_ldfslp .AND. .NOT.l k_c1d ) THEN ! slopes323 IF( l_ldfslp .AND. .NOT.ln_c1d ) THEN ! slopes 324 324 IF( sf_dyn(jf_tem)%ln_tint ) THEN ! time interpolation 325 325 ALLOCATE( uslpdta (jpi,jpj,jpk,2), vslpdta (jpi,jpj,jpk,2), & … … 628 628 & / REAL( sf_dyn(jf_tem)%nrec(2,sf_dyn(jf_tem)%naa) - sf_dyn(jf_tem)%nrec(2,sf_dyn(jf_tem)%nbb), wp ) 629 629 ztintb = 1. - ztinta 630 IF( l_ldfslp .AND. .NOT.l k_c1d ) THEN ! Computes slopes (here avt is used as workspace)630 IF( l_ldfslp .AND. .NOT.ln_c1d ) THEN ! Computes slopes (here avt is used as workspace) 631 631 uslp (:,:,:) = ztintb * uslpdta (:,:,:,1) + ztinta * uslpdta (:,:,:,2) 632 632 vslp (:,:,:) = ztintb * vslpdta (:,:,:,1) + ztinta * vslpdta (:,:,:,2) … … 640 640 CALL compute_slopes( kt, zts, zuslp, zvslp, zwslpi, zwslpj, Kbb, Kmm ) 641 641 ! 642 IF( l_ldfslp .AND. .NOT.l k_c1d ) THEN ! Computes slopes (here avt is used as workspace)642 IF( l_ldfslp .AND. .NOT.ln_c1d ) THEN ! Computes slopes (here avt is used as workspace) 643 643 uslp (:,:,:) = zuslp (:,:,:) 644 644 vslp (:,:,:) = zvslp (:,:,:) … … 665 665 !!--------------------------------------------------------------------- 666 666 ! 667 IF( l_ldfslp .AND. .NOT.l k_c1d ) THEN ! Computes slopes (here avt is used as workspace)667 IF( l_ldfslp .AND. .NOT.ln_c1d ) THEN ! Computes slopes (here avt is used as workspace) 668 668 CALL eos ( pts, rhd, rhop, gdept_0(:,:,:) ) 669 669 CALL eos_rab( pts, rab_n, Kmm ) ! now local thermal/haline expension ratio at T-points
Note: See TracChangeset
for help on using the changeset viewer.