Changeset 4034
 Timestamp:
 20130925T14:02:06+02:00 (10 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r3963 r4034 63 63 REAL(wp), PARAMETER :: krdgmin = 1.1_wp ! min ridge thickness multiplier 64 64 REAL(wp), PARAMETER :: kraft = 2.0_wp ! rafting multipliyer 65 REAL(wp), PARAMETER :: kamax = 1.0 65 66 66 67 REAL(wp) :: Cp ! … … 227 228 ! to give asum = 1.0 after ridging. 228 229 229 divu_adv(ji,jj) = ( amax  asum(ji,jj) ) / rdt_ice ! asum found in ridgeprep230 divu_adv(ji,jj) = ( kamax  asum(ji,jj) ) / rdt_ice ! asum found in ridgeprep 230 231 231 232 IF( divu_adv(ji,jj) < 0._wp ) closing_net(ji,jj) = MAX( closing_net(ji,jj), divu_adv(ji,jj) ) … … 312 313 DO jj = 1, jpj 313 314 DO ji = 1, jpi 314 IF (ABS(asum(ji,jj)  amax ) .LT. epsi11) THEN315 IF (ABS(asum(ji,jj)  kamax ) .LT. epsi11) THEN 315 316 closing_net(ji,jj) = 0._wp 316 317 opning (ji,jj) = 0._wp 317 318 ELSE 318 319 iterate_ridging = 1 319 divu_adv (ji,jj) = ( amax  asum(ji,jj) ) / rdt_ice320 divu_adv (ji,jj) = ( kamax  asum(ji,jj) ) / rdt_ice 320 321 closing_net(ji,jj) = MAX( 0._wp, divu_adv(ji,jj) ) 321 322 opning (ji,jj) = MAX( 0._wp, divu_adv(ji,jj) ) … … 356 357 DO ji = 1, jpi 357 358 358 IF (ABS(asum(ji,jj)  amax) .GT. epsi11) asum_error = .true.359 IF (ABS(asum(ji,jj)  kamax) .GT. epsi11) asum_error = .true. 359 360 360 361 dardg1dt(ji,jj) = dardg1dt(ji,jj) * dti … … 375 376 DO jj = 1, jpj 376 377 DO ji = 1, jpi 377 IF (ABS(asum(ji,jj)  amax) .GT. epsi11) THEN ! there is a bug378 IF (ABS(asum(ji,jj)  kamax) .GT. epsi11) THEN ! there is a bug 378 379 WRITE(numout,*) ' ' 379 380 WRITE(numout,*) ' ALERTE : Ridging error: total area = ', asum(ji,jj) … … 482 483 IF ( ABS( zchk_smv ) > 1.e4 ) WRITE(numout,*) 'violation saline [psu*m3/day] (limitd_me) = ',(zchk_smv * 86400.) 483 484 IF ( zchk_vmin < 0. ) WRITE(numout,*) 'violation v_i<0 [mm] (limitd_me) = ',(zchk_vmin * 1.e3) 484 IF ( zchk_amax > amax+epsi10 ) WRITE(numout,*) 'violation a_i>amax (limitd_me) = ',zchk_amax485 IF ( zchk_amax > kamax+epsi10 ) WRITE(numout,*) 'violation a_i>amax (limitd_me) = ',zchk_amax 485 486 IF ( zchk_amin < 0. ) WRITE(numout,*) 'violation a_i<0 (limitd_me) = ',zchk_amin 486 487 ENDIF … … 1153 1154 afrft(ji,jj) = arft1(ji,jj) / aicen_init(ji,jj,jl1) !rafting 1154 1155 1155 IF (afrac(ji,jj) > amax + epsi11) THEN !riging1156 IF (afrac(ji,jj) > kamax + epsi11) THEN !riging 1156 1157 large_afrac = .true. 1157 ELSEIF (afrac(ji,jj) > amax) THEN ! roundoff error1158 afrac(ji,jj) = amax1158 ELSEIF (afrac(ji,jj) > kamax) THEN ! roundoff error 1159 afrac(ji,jj) = kamax 1159 1160 ENDIF 1160 IF (afrft(ji,jj) > amax + epsi11) THEN !rafting1161 IF (afrft(ji,jj) > kamax + epsi11) THEN !rafting 1161 1162 large_afrft = .true. 1162 ELSEIF (afrft(ji,jj) > amax) THEN ! roundoff error1163 afrft(ji,jj) = amax1163 ELSEIF (afrft(ji,jj) > kamax) THEN ! roundoff error 1164 afrft(ji,jj) = kamax 1164 1165 ENDIF 1165 1166 … … 1309 1310 ji = indxi(ij) 1310 1311 jj = indxj(ij) 1311 IF( afrac(ji,jj) > amax + epsi11 ) THEN1312 IF( afrac(ji,jj) > kamax + epsi11 ) THEN 1312 1313 WRITE(numout,*) '' 1313 1314 WRITE(numout,*) ' ardg > a_i' … … 1321 1322 ji = indxi(ij) 1322 1323 jj = indxj(ij) 1323 IF( afrft(ji,jj) > amax + epsi11 ) THEN1324 IF( afrft(ji,jj) > kamax + epsi11 ) THEN 1324 1325 WRITE(numout,*) '' 1325 1326 WRITE(numout,*) ' arft > a_i'
Note: See TracChangeset
for help on using the changeset viewer.