Changeset 10365 for NEMO/branches/2018
- Timestamp:
- 2018-12-01T12:33:31+01:00 (5 years ago)
- Location:
- NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/ICE/icedyn_rhg_evp.F90
r10345 r10365 143 143 REAL(wp), DIMENSION(jpi,jpj) :: zds ! shear 144 144 REAL(wp), DIMENSION(jpi,jpj) :: zs1, zs2, zs12 ! stress tensor components 145 REAL(wp), DIMENSION(jpi,jpj) :: zu_ice, zv_ice, zresr ! check convergence145 !!$ REAL(wp), DIMENSION(jpi,jpj) :: zu_ice, zv_ice, zresr ! check convergence 146 146 REAL(wp), DIMENSION(jpi,jpj) :: zssh_lead_m ! array used for the calculation of ice surface slope: 147 147 ! ! ocean surface (ssh_m) if ice is not embedded … … 324 324 DO jter = 1 , nn_nevp ! loop over jter ! 325 325 ! !----------------------! 326 IF(ln_ctl) THEN ! Convergence test327 DO jj = 1, jpjm1328 zu_ice(:,jj) = u_ice(:,jj) ! velocity at previous time step329 zv_ice(:,jj) = v_ice(:,jj)330 END DO331 ENDIF326 !!$ IF(ln_ctl) THEN ! Convergence test 327 !!$ DO jj = 1, jpjm1 328 !!$ zu_ice(:,jj) = u_ice(:,jj) ! velocity at previous time step 329 !!$ zv_ice(:,jj) = v_ice(:,jj) 330 !!$ END DO 331 !!$ ENDIF 332 332 333 333 ! --- divergence, tension & shear (Appendix B of Hunke & Dukowicz, 2002) --- ! … … 636 636 ! 637 637 ENDIF 638 639 IF(ln_ctl) THEN ! Convergence test640 DO jj = 2 , jpjm1641 zresr(:,jj) = MAX( ABS( u_ice(:,jj) - zu_ice(:,jj) ), ABS( v_ice(:,jj) - zv_ice(:,jj) ) )642 END DO643 zresm = MAXVAL( zresr( 1:jpi, 2:jpjm1 ) )644 IF( lk_mpp ) CALL mpp_max( 'icedyn_rhg_evp', zresm ) ! max over the global domain645 ENDIF638 639 !!$ IF(ln_ctl) THEN ! Convergence test 640 !!$ DO jj = 2 , jpjm1 641 !!$ zresr(:,jj) = MAX( ABS( u_ice(:,jj) - zu_ice(:,jj) ), ABS( v_ice(:,jj) - zv_ice(:,jj) ) ) 642 !!$ END DO 643 !!$ zresm = MAXVAL( zresr( 1:jpi, 2:jpjm1 ) ) 644 !!$ IF( lk_mpp ) CALL mpp_max( 'icedyn_rhg_evp', zresm ) ! max over the global domain 645 !!$ ENDIF 646 646 ! 647 647 ! ! ==================== ! -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/SBC/sbcfwb.F90
r10314 r10365 154 154 WHERE( erp < 0._wp ) ztmsk_pos = 0._wp 155 155 ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 156 !157 zsurf_neg = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_neg(:,:) ) ! Area filled by <0 and >0 erp158 zsurf_pos = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_pos(:,:) )159 156 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 160 157 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 161 158 ! 162 159 IF( z_fwf < 0._wp ) THEN ! spread out over >0 erp area to increase evaporation 163 zsurf_tospread = zsurf_pos 164 ztmsk_tospread(:,:) = ztmsk_pos(:,:) 160 zsurf_pos = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_pos(:,:) ) 161 zsurf_tospread = zsurf_pos 162 ztmsk_tospread(:,:) = ztmsk_pos(:,:) 165 163 ELSE ! spread out over <0 erp area to increase precipitation 166 zsurf_tospread = zsurf_neg 167 ztmsk_tospread(:,:) = ztmsk_neg(:,:) 164 zsurf_neg = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_neg(:,:) ) ! Area filled by <0 and >0 erp 165 zsurf_tospread = zsurf_neg 166 ztmsk_tospread(:,:) = ztmsk_neg(:,:) 168 167 ENDIF 169 168 ! -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zflx.F90
r10314 r10365 172 172 END DO 173 173 174 t_oce_co2_flx = glob_sum( 'p4zflx', oce_co2(:,:) ) ! Total Flux of Carbon 174 IF( iom_use("tcflx") .OR. iom_use("tcflxcum") .OR. kt == nitrst & 175 & .OR. (ln_check_mass .AND. kt == nitend) ) & 176 t_oce_co2_flx = glob_sum( 'p4zflx', oce_co2(:,:) ) ! Total Flux of Carbon 175 177 t_oce_co2_flx_cum = t_oce_co2_flx_cum + t_oce_co2_flx ! Cumulative Total Flux of Carbon 176 178 ! t_atm_co2_flx = glob_sum( 'p4zflx', satmco2(:,:) * e1e2t(:,:) ) ! Total atmospheric pCO2 … … 205 207 CALL iom_put( "Dpo2" , zw2d ) 206 208 ENDIF 207 IF( iom_use( "tcflx" ) )CALL iom_put( "tcflx" , t_oce_co2_flx * rfact2r ) ! molC/s208 CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum )! molC209 CALL iom_put( "tcflx" , t_oce_co2_flx * rfact2r ) ! molC/s 210 CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum ) ! molC 209 211 ! 210 212 DEALLOCATE( zw2d ) -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/TRP/trcrad.F90
r10314 r10365 153 153 ztrcorn = glob_sum( 'trcrad', MIN( 0., ptrn(:,:,:,jn) ) * cvol(:,:,:) ) 154 154 ! 155 ztrmasb = glob_sum( 'trcrad', MAX( 0., ptrb(:,:,:,jn) ) * cvol(:,:,:) )156 ztrmasn = glob_sum( 'trcrad', MAX( 0., ptrn(:,:,:,jn) ) * cvol(:,:,:) )157 !158 155 IF( ztrcorb /= 0 ) THEN 156 ztrmasb = glob_sum( 'trcrad', MAX( 0., ptrb(:,:,:,jn) ) * cvol(:,:,:) ) 159 157 zcoef = 1. + ztrcorb / ztrmasb 160 158 DO jk = 1, jpkm1 … … 165 163 ! 166 164 IF( ztrcorn /= 0 ) THEN 165 ztrmasn = glob_sum( 'trcrad', MAX( 0., ptrn(:,:,:,jn) ) * cvol(:,:,:) ) 167 166 zcoef = 1. + ztrcorn / ztrmasn 168 167 DO jk = 1, jpkm1
Note: See TracChangeset
for help on using the changeset viewer.