- Timestamp:
- 2018-04-13T19:03:31+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/NST_SRC/agrif_lim3_update.F90
r9454 r9482 1 1 #define TWO_WAY 2 !!#undef TWO_WAY 2 3 3 4 MODULE agrif_lim3_update … … 32 33 PRIVATE 33 34 34 PUBLIC agrif_update_lim3 ! called by agrif_user.F90 35 PUBLIC agrif_update_lim3 ! called by agrif_user.F90 and icestp.F90 35 36 36 37 !!---------------------------------------------------------------------- … … 51 52 !!---------------------------------------------------------------------- 52 53 ! 53 !! clem: I think the update should take place each time the ocean sees the surface forcings 54 !! (but maybe I am wrong and we should update every rhot time steps) 55 IF( ( MOD( (kt-nit000)/nn_fsbc + 1, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) /=0 ) .AND. (kt /= 0) ) RETURN ! do not update if nb of child time steps differ from time refinement 56 ! i.e. update only at the parent time step 57 IF( nn_ice == 0 ) RETURN ! do not update if child domain does not have ice 54 IF( Agrif_Root() .OR. nn_ice == 0 ) RETURN ! do not update if inside Parent Grid or if child domain does not have ice 55 ! 56 IF( ( MOD( (kt-nit000)/nn_fsbc + 1, Agrif_irhot() * Agrif_Parent(nn_fsbc)/nn_fsbc ) /=0 ) .AND. (kt /= 0) ) RETURN ! update only at the parent ice time step 58 57 ! 59 58 Agrif_SpecialValueFineGrid = -9999. … … 148 147 END DO 149 148 ! 150 ! integrated values151 vt_i(i1:i2,j1:j2) = SUM( v_i(i1:i2,j1:j2,:) , dim=3 )152 vt_s(i1:i2,j1:j2) = SUM( v_s(i1:i2,j1:j2,:) , dim=3 )153 at_i(i1:i2,j1:j2) = SUM( a_i(i1:i2,j1:j2,:) , dim=3 )154 et_s(i1:i2,j1:j2) = SUM( SUM( e_s(i1:i2,j1:j2,:,:), dim=4 ), dim=3 )155 et_i(i1:i2,j1:j2) = SUM( SUM( e_i(i1:i2,j1:j2,:,:), dim=4 ), dim=3 )156 157 at_ip(i1:i2,j1:j2) = SUM( a_ip(i1:i2,j1:j2,:), dim=3 ) ! melt ponds158 vt_ip(i1:i2,j1:j2) = SUM( v_ip(i1:i2,j1:j2,:), dim=3 )159 !160 ato_i(i1:i2,j1:j2) = 1._wp - at_i(i1:i2,j1:j2) ! open water fraction161 162 149 DO jl = 1, jpl 163 150 WHERE( tmask(i1:i2,j1:j2,1) == 0._wp ) t_su(i1:i2,j1:j2,jl) = rt0 ! to avoid a division by 0 in sbcblk.F90
Note: See TracChangeset
for help on using the changeset viewer.