Changeset 10106
- Timestamp:
- 2018-09-07T15:53:19+02:00 (6 years ago)
- Location:
- NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crs.F90
r7795 r10106 120 120 REAL(wp), DIMENSION(:,:), ALLOCATABLE,SAVE :: ff_crs 121 121 INTEGER, DIMENSION(:,:), ALLOCATABLE,SAVE :: mbathy_crs, mbkt_crs, mbku_crs, mbkv_crs 122 INTEGER, DIMENSION(:,:), ALLOCATABLE,SAVE :: mikt_crs 122 123 123 124 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE,SAVE :: gdept_0_crs, gdepu_0_crs, gdepv_0_crs, gdepw_0_crs … … 169 170 REAL(wp), PUBLIC, DIMENSION(:,:) , ALLOCATABLE :: rnf_crs,rnf_b_crs,h_rnf_crs 170 171 INTEGER , PUBLIC, DIMENSION(:,:) , ALLOCATABLE :: nk_rnf_crs 172 REAL(wp), PUBLIC, DIMENSION(:,:) , ALLOCATABLE :: fwfisf_b_crs, fwfisf_crs 171 173 REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: trc_i_crs,trc_o_crs 172 174 REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: sbc_trc_crs, sbc_trc_b_crs … … 284 286 285 287 ALLOCATE( mbathy_crs(jpi_crs,jpj_crs), mbkt_crs(jpi_crs,jpj_crs) , & 286 & mbku_crs(jpi_crs,jpj_crs) , mbkv_crs(jpi_crs,jpj_crs) , STAT=ierr(9)) 288 & mbku_crs(jpi_crs,jpj_crs) , mbkv_crs(jpi_crs,jpj_crs) , & 289 & mikt_crs(jpi_crs,jpj_crs) , STAT=ierr(9)) 287 290 288 291 ALLOCATE( gdept_0_crs(jpi_crs,jpj_crs,jpk), gdepu_0_crs(jpi_crs,jpj_crs,jpk) , & … … 310 313 & rnf_crs (jpi_crs,jpj_crs), rnf_b_crs(jpi_crs ,jpj_crs), nk_rnf_crs(jpi_crs ,jpj_crs), h_rnf_crs(jpi_crs ,jpj_crs), & 311 314 & emp_crs (jpi_crs,jpj_crs), emp_b_crs(jpi_crs,jpj_crs), & 315 & fwfisf_crs(jpi_crs,jpj_crs), fwfisf_b_crs(jpi_crs,jpj_crs), & 312 316 & sbc_trc_crs (jpi_crs,jpj_crs,jpts), sbc_trc_b_crs(jpi_crs,jpj_crs,jpts), & 313 317 & trc_i_crs (jpi_crs,jpj_crs,jpts), trc_o_crs(jpi_crs,jpj_crs,jpts), & -
NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90
r10101 r10106 1736 1736 END SUBROUTINE crs_dom_bat 1737 1737 1738 SUBROUTINE crs_dom_top 1739 !!---------------------------------------------------------------- 1740 !! *** SUBROUTINE crs_dom_top *** 1741 !! 1742 !! ** Purpose : coarsening mikt 1743 !! 1744 !!---------------------------------------------------------------- 1745 INTEGER :: ji, jj, jk 1746 INTEGER :: ijis, ijie, ijjs, ijje 1747 1748 !!---------------------------------------------------------------- 1749 1750 DO jj = nldj_crs,nlej_crs 1751 1752 ijjs = mjs_crs(jj) 1753 ijje = mje_crs(jj) 1754 1755 DO ji = nldi_crs, nlei_crs 1756 1757 ijis = mis_crs(ji) 1758 ijie = mie_crs(ji) 1759 1760 mikt_crs(ji,jj) = MINVAL( mikt(ijis:ijie,ijjs:ijje) ) 1761 ENDDO 1762 ENDDO 1763 1764 mikt_crs = MAX( mikt_crs , 1 ) 1765 1766 END SUBROUTINE crs_dom_top 1767 1738 1768 SUBROUTINE PIKSRT(N,ARR) 1739 1769 INTEGER ,INTENT(IN) :: N -
NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsfld.F90
r8894 r10106 156 156 CALL crs_dom_ope( rnf_b ,'SUM', 'T', tmask, rnf_b_crs , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 157 157 CALL crs_dom_ope( emp_b ,'SUM', 'T', tmask, emp_b_crs , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 158 CALL crs_dom_ope( fwfisf_b ,'SUM', 'T', tmask, fwfisf_b_crs , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 158 159 159 160 ! Temperature … … 303 304 CALL crs_dom_ope( emp ,'SUM', 'T', tmask, emp_crs , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 304 305 CALL crs_dom_ope( fmmflx,'SUM', 'T', tmask, fmmflx_crs, p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 306 CALL crs_dom_ope( fwfisf,'SUM', 'T', tmask, fwfisf_crs , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 305 307 306 308 CALL crs_dom_ope( fr_i ,'SUM', 'T', tmask, fr_i_crs , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) -
NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90
r7806 r10106 460 460 ztc_f = ztc_n + atfp * ztc_d 461 461 ! 462 IF( jk == 1 ) THEN ! first level 463 ze3t_f = ze3t_f - zfact2 * ( emp_b_crs(ji,jj) - emp_crs(ji,jj) + rnf_crs(ji,jj) - rnf_b_crs(ji,jj) ) 462 IF( jk == mikt_crs(ji,jj) ) THEN ! first level 463 ze3t_f = ze3t_f - zfact2 * ( (emp_b_crs(ji,jj) - emp_crs(ji,jj) ) & 464 & - (rnf_b_crs(ji,jj) - rnf_crs(ji,jj) ) & 465 & + (fwfisf_b_crs(ji,jj) - fwfisf_crs(ji,jj)) ) 464 466 ztc_f = ztc_f - zfact1 * ( psbc_tc(ji,jj,jn) - psbc_tc_b(ji,jj,jn) ) 465 467 ENDIF 466 !cbr as it is a subroutine dedicated to crs, TRA options are not necessary 467 !cbr IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) & ! solar penetration (temperature only) 468 !cbr & ztc_f = ztc_f - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) ) 469 !cbr 470 !cbr IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) ) & ! river runoffs 471 !cbr & ztc_f = ztc_f - zfact1 * ( rnf_tsc(ji,jj,jn) - rnf_tsc_b(ji,jj,jn) ) & 472 !cbr & * fse3t_n(ji,jj,jk) / h_rnf(ji,jj) 468 469 !cbr: not adapted for coarsening for the moment 470 ! ice shelf 471 !IF( ll_isf ) THEN 472 ! ! level fully include in the Losch_2008 ice shelf boundary layer 473 ! IF ( jk >= misfkt(ji,jj) .AND. jk < misfkb(ji,jj) ) & 474 ! ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) & 475 ! & * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) 476 ! ! level partially include in Losch_2008 ice shelf boundary layer 477 ! IF ( jk == misfkb(ji,jj) ) & 478 ! ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) & 479 ! & * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) * ralpha(ji,jj) 480 !END IF 473 481 474 482 ze3t_f = 1.e0 / ze3t_f … … 476 484 ptn(ji,jj,jk,jn) = pta(ji,jj,jk,jn) ! ptn <-- pta 477 485 ! 478 IF( ll_tra_hpg ) THEN ! semi-implicit hpg (T & S only)479 ze3t_d = 1.e0 / ( ze3t_n + rbcp * ze3t_d )480 pta(ji,jj,jk,jn) = ze3t_d * ( ztc_n + rbcp * ztc_d ) ! ta <-- Brown & Campana average481 ENDIF482 486 END DO 483 487 END DO -
NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90
r10104 r10106 122 122 USE dom_oce , ONLY : ln_sco => ln_sco 123 123 USE dom_oce , ONLY : neuler => neuler 124 USE dom_oce , ONLY : mikt => mikt !: f-points (m)125 124 126 125 USE crs, ONLY : mi0 => mi0 … … 193 192 USE crs , ONLY : e3t_0 => e3t_0_crs !: depth of t-points (m) 194 193 USE crs , ONLY : mbkt => mbkt_crs !: vertical index of the bottom last T- ocean level 194 USE crs , ONLY : mikt => mikt_crs !: vertical index of the top first T- ocean level 195 195 USE crs , ONLY : tmask_i => tmask_i_crs !: Interior mask at t-points 196 196 USE crs , ONLY : tmask => tmask_crs !: land/ocean mask at t-points
Note: See TracChangeset
for help on using the changeset viewer.