- Timestamp:
- 2017-07-10T20:24:21+02:00 (7 years ago)
- Location:
- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/ASM
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90
r8306 r8313 146 146 CALL iom_rstput( kt, nitdin_r, inum, 'sshn' , sshn ) 147 147 #if defined key_lim3 148 IF( nn_ice == 2 .OR. nn_ice ==3 ) THEN149 IF( ALLOCATED( frld) ) THEN150 CALL iom_rstput( kt, nitdin_r, inum, 'iceconc', 1._wp - frld(:,:) )148 IF( nn_ice == 3 ) THEN 149 IF( ALLOCATED(at_i) ) THEN 150 CALL iom_rstput( kt, nitdin_r, inum, 'iceconc', at_i(:,:) ) 151 151 ELSE 152 CALL ctl_warn('Ice concentration not written to background as ice variable frldnot allocated on this timestep')152 CALL ctl_warn('Ice concentration not written to background as ice variable at_i not allocated on this timestep') 153 153 ENDIF 154 154 ENDIF -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r8306 r8313 36 36 USE lib_mpp ! MPP library 37 37 #if defined key_lim3 38 USE ice, ONLY: htm_i, at_i, pfrld, phicif38 USE ice, ONLY: htm_i, at_i, at_i_b 39 39 #endif 40 40 USE sbc_oce ! Surface boundary condition variables. … … 837 837 zohicif(:,:) = htm_i(:,:) 838 838 ! 839 at_i (:,:) = 1. - MIN( MAX( 1.-at_i(:,:) - seaice_bkginc(:,:) * zincwgt, 0.0_wp), 1.0_wp)840 pfrld = MIN( MAX( pfrld(:,:)- seaice_bkginc(:,:) * zincwgt, 0.0_wp), 1.0_wp)839 at_i (:,:) = 1. - MIN( MAX( 1.-at_i (:,:) - seaice_bkginc(:,:) * zincwgt, 0.0_wp), 1.0_wp) 840 at_i_b(:,:) = 1. - MIN( MAX( 1.-at_i_b(:,:) - seaice_bkginc(:,:) * zincwgt, 0.0_wp), 1.0_wp) 841 841 fr_i(:,:) = at_i(:,:) ! adjust ice fraction 842 842 ! … … 852 852 ! nudge ice depth 853 853 htm_i (:,:) = htm_i (:,:) + zhicifinc(:,:) 854 phicif(:,:) = phicif(:,:) + zhicifinc(:,:)855 854 ! 856 855 ! seaice salinity balancing (to add) … … 889 888 ! Initialize the now fields the background + increment 890 889 at_i(:,:) = 1. - MIN( MAX( 1.-at_i(:,:) - seaice_bkginc(:,:), 0.0_wp), 1.0_wp) 891 pfrld(:,:) = 1.-at_i(:,:)890 at_i_b(:,:) = at_i(:,:) 892 891 fr_i(:,:) = at_i(:,:) ! adjust ice fraction 893 892 ! … … 903 902 ! nudge ice depth 904 903 htm_i (:,:) = htm_i (:,:) + zhicifinc(:,:) 905 phicif(:,:) = phicif(:,:)906 904 ! 907 905 ! seaice salinity balancing (to add)
Note: See TracChangeset
for help on using the changeset viewer.