- Timestamp:
- 2021-09-10T11:22:28+02:00 (21 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_FOAM_IAU_SI3_SIC/src/OCE/ASM/asminc.F90
r15241 r15242 933 933 END WHERE 934 934 935 WRITE(numout,*) "zhi_damin = ",zhi_damin936 935 ! add positive concentration increments with zhi_damin to regions where ice 937 936 ! is already present and bound them to 1 … … 1004 1003 rn_hinew = rn_hinew_save 1005 1004 1006 ! compute thickness now1007 WHERE( a_i(:,:,:) > epsi10 )1008 z1_a_i(:,:,:) = 1.0_wp/a_i(:,:,:)1009 ELSEWHERE1010 z1_a_i(:,:,:) = 0.0_wp1011 END WHERE1012 h_i(:,:,:) = v_i(:,:,:) * z1_a_i(:,:,:)1013 1014 1005 incr_newice(:,:,:) = .FALSE. 1015 1006 initial_step = .FALSE. 1016 1007 END IF 1017 1018 ! make volume consistent with concentration and thickness1019 ! note where conc is zero this will lead to zero volume1020 ! and soon zero thickness after the call to ice_var_zapsmall1021 v_i (:,:,:) = h_i (:,:,:) * a_i (:,:,:)1022 1008 1023 1009 ! as ice concentration has changed, maintain equivalent values for prognostic variables
Note: See TracChangeset
for help on using the changeset viewer.