Ignore:
Timestamp:
2020-03-24T22:25:05+01:00 (8 months ago)
Author:
dcarneir
Message:

Removing commented code in subroutine sit_asm_inc

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_sit/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r11952 r12595  
    10241024      INTEGER  ::   it 
    10251025      REAL(wp) ::   zincwgt   ! IAU weight for current time step 
    1026 ! #if defined key_lim2 
    1027 !       REAL(wp), DIMENSION(jpi,jpj) ::   zofrld, zohicif, zseaicendg, zhicifinc  ! LIM 
    1028 !       REAL(wp) ::   zhicifmin = 0.5_wp      ! ice minimum depth in metres 
    1029 !       !!THICKNESS INCS NOT SET UP FOR LIM 
    1030 ! #endif 
    10311026      !!---------------------------------------------------------------------- 
    10321027 
     
    10511046            ENDIF 
    10521047 
    1053             ! Sea-ice : LIM-3 case (to add) 
    1054  
    1055 ! #if defined key_lim2 
    1056 !             ! Sea-ice : LIM-2 case (to add if needed) 
    1057 !             zofrld (:,:) = frld(:,:) 
    1058 !             zohicif(:,:) = hicif(:,:) 
    1059 !             ! 
    1060 !             frld  = MIN( MAX( frld (:,:) - seaice_bkginc(:,:) * zincwgt, 0.0_wp), 1.0_wp) 
    1061 !             pfrld = MIN( MAX( pfrld(:,:) - seaice_bkginc(:,:) * zincwgt, 0.0_wp), 1.0_wp) 
    1062 !             fr_i(:,:) = 1.0_wp - frld(:,:)        ! adjust ice fraction 
    1063 !             ! 
    1064 !             zseaicendg(:,:) = zofrld(:,:) - frld(:,:)   ! find out actual sea ice nudge applied 
    1065 !             ! 
    1066 !             ! Nudge sea ice depth to bring it up to a required minimum depth 
    1067 !             WHERE( zseaicendg(:,:) > 0.0_wp .AND. hicif(:,:) < zhicifmin )  
    1068 !                zhicifinc(:,:) = (zhicifmin - hicif(:,:)) * zincwgt     
    1069 !             ELSEWHERE 
    1070 !                zhicifinc(:,:) = 0.0_wp 
    1071 !             END WHERE 
    1072 !             ! 
    1073 !             ! nudge ice depth 
    1074 !             hicif (:,:) = hicif (:,:) + zhicifinc(:,:) 
    1075 !             phicif(:,:) = phicif(:,:) + zhicifinc(:,:)        
    1076 !             ! 
    1077 !             ! seaice salinity balancing (to add) 
    1078 ! #endif 
    1079  
    10801048#if defined key_cice && defined key_asminc 
    10811049            ! Sea-ice thickness : CICE case. Pass ice thickness increment tendency into CICE 
     
    11061074            neuler = 0                    ! Force Euler forward step 
    11071075 
    1108             ! Sea-ice : LIM-3 case (to add) 
    1109  
    1110 ! #if defined key_lim2 
    1111 !             ! Sea-ice : LIM-2 case (add if needed) 
    1112 !             zofrld(:,:)=frld(:,:) 
    1113 !             zohicif(:,:)=hicif(:,:) 
    1114 !             !  
    1115 !             ! Initialize the now fields the background + increment 
    1116 !             frld (:,:) = MIN( MAX( frld(:,:) - seaice_bkginc(:,:), 0.0_wp), 1.0_wp) 
    1117 !             pfrld(:,:) = frld(:,:)  
    1118 !             fr_i (:,:) = 1.0_wp - frld(:,:)                ! adjust ice fraction 
    1119 !             zseaicendg(:,:) = zofrld(:,:) - frld(:,:)      ! find out actual sea ice nudge applied 
    1120 !             ! 
    1121 !             ! Nudge sea ice depth to bring it up to a required minimum depth 
    1122 !             WHERE( zseaicendg(:,:) > 0.0_wp .AND. hicif(:,:) < zhicifmin )  
    1123 !                zhicifinc(:,:) = (zhicifmin - hicif(:,:)) * zincwgt     
    1124 !             ELSEWHERE 
    1125 !                zhicifinc(:,:) = 0.0_wp 
    1126 !             END WHERE 
    1127 !             ! 
    1128 !             ! nudge ice depth 
    1129 !             hicif (:,:) = hicif (:,:) + zhicifinc(:,:) 
    1130 !             phicif(:,:) = phicif(:,:)        
    1131 !             ! 
    1132 !             ! seaice salinity balancing (to add) 
    1133 ! #endif 
    1134   
    11351076#if defined key_cice && defined key_asminc 
    11361077            ! Sea-ice thickness : CICE case. Pass ice thickness increment tendency into CICE 
     
    11491090          
    11501091         ENDIF 
    1151  
    1152 !#if defined defined key_lim2 || defined key_cice 
    1153 ! 
    1154 !            IF (ln_seaicebal ) THEN        
    1155 !             !! balancing salinity increments 
    1156 !             !! simple case from limflx.F90 (doesn't include a mass flux) 
    1157 !             !! assumption is that as ice concentration is reduced or increased 
    1158 !             !! the snow and ice depths remain constant 
    1159 !             !! note that snow is being created where ice concentration is being increased 
    1160 !             !! - could be more sophisticated and 
    1161 !             !! not do this (but would need to alter h_snow) 
    1162 ! 
    1163 !             usave(:,:,:)=sb(:,:,:)   ! use array as a temporary store 
    1164 ! 
    1165 !             DO jj = 1, jpj 
    1166 !               DO ji = 1, jpi  
    1167 !           ! calculate change in ice and snow mass per unit area 
    1168 !           ! positive values imply adding salt to the ocean (results from ice formation) 
    1169 !           ! fwf : ice formation and melting 
    1170 ! 
    1171 !                 zfons = ( -nfresh_da(ji,jj)*soce + nfsalt_da(ji,jj) )*rdt 
    1172 ! 
    1173 !           ! change salinity down to mixed layer depth 
    1174 !                 mld=hmld_kara(ji,jj) 
    1175 ! 
    1176 !           ! prevent small mld 
    1177 !           ! less than 10m can cause salinity instability  
    1178 !                 IF (mld < 10) mld=10 
    1179 ! 
    1180 !           ! set to bottom of a level  
    1181 !                 DO jk = jpk-1, 2, -1 
    1182 !                   IF ((mld > gdepw(ji,jj,jk)) .and. (mld < gdepw(ji,jj,jk+1))) THEN  
    1183 !                     mld=gdepw(ji,jj,jk+1) 
    1184 !                     jkmax=jk 
    1185 !                   ENDIF 
    1186 !                 ENDDO 
    1187 ! 
    1188 !            ! avoid applying salinity balancing in shallow water or on land 
    1189 !            !  
    1190 ! 
    1191 !            ! dsal_ocn (psu kg m^-2) / (kg m^-3 * m) 
    1192 ! 
    1193 !                 dsal_ocn=0.0_wp 
    1194 !                 sal_thresh=5.0_wp        ! minimum salinity threshold for salinity balancing 
    1195 ! 
    1196 !                 if (tmask(ji,jj,1) > 0 .AND. tmask(ji,jj,jkmax) > 0 ) & 
    1197 !                              dsal_ocn = zfons / (rhop(ji,jj,1) * mld) 
    1198 ! 
    1199 !           ! put increments in for levels in the mixed layer 
    1200 !           ! but prevent salinity below a threshold value  
    1201 ! 
    1202 !                   DO jk = 1, jkmax               
    1203 ! 
    1204 !                     IF (dsal_ocn > 0.0_wp .or. sb(ji,jj,jk)+dsal_ocn > sal_thresh) THEN  
    1205 !                           sb(ji,jj,jk) = sb(ji,jj,jk) + dsal_ocn 
    1206 !                           sn(ji,jj,jk) = sn(ji,jj,jk) + dsal_ocn 
    1207 !                     ENDIF 
    1208 ! 
    1209 !                   ENDDO 
    1210 ! 
    1211 !      !            !  salt exchanges at the ice/ocean interface 
    1212 !      !            zpmess         = zfons / rdt_ice    ! rdt_ice is ice timestep 
    1213 !      ! 
    1214 !      !! Adjust fsalt. A +ve fsalt means adding salt to ocean 
    1215 !      !!           fsalt(ji,jj) =  fsalt(ji,jj) + zpmess     ! adjust fsalt   
    1216 !      !!                
    1217 !      !!           emps(ji,jj) = emps(ji,jj) + zpmess        ! or adjust emps (see icestp1d)  
    1218 !      !!                                                     ! E-P (kg m-2 s-2) 
    1219 !      !            emp(ji,jj) = emp(ji,jj) + zpmess          ! E-P (kg m-2 s-2) 
    1220 !               ENDDO !ji 
    1221 !             ENDDO !jj! 
    1222 ! 
    1223 !            ENDIF !ln_seaicebal 
    1224 ! 
    1225 !#endif 
    12261092 
    12271093      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.