- Timestamp:
- 2012-04-30T13:24:25+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90
r3372 r3374 106 106 DO ji = 1, jpi 107 107 src_calving_hflx(ji,jj) = narea 108 src_calving (ji,jj) = nicbpack * (njmpp+jj-1) + nimpp+ji-1 !!gm use here mjg(jj) and mig(ji)... 109 !!gm src_calving (ji,jj) = nicbpack * mjg(jj) + mig(ji) 108 src_calving (ji,jj) = nicbpack * mjg(jj) + mig(ji) 110 109 END DO 111 110 END DO … … 124 123 i1 = i3 - i2*nicbpack 125 124 i3 = INT( src_calving_hflx(ji,jj) ) 126 IF( i1 == nimpp+ji-1.AND. i3 == narea ) THEN125 IF( i1 == mig(ji) .AND. i3 == narea ) THEN 127 126 IF( nicbdi < 0 ) THEN ; nicbdi = ji 128 127 ELSE ; nicbei = ji … … 140 139 i1 = i3 - i2*nicbpack 141 140 i3 = INT( src_calving_hflx(ji,jj) ) 142 IF( i2 == njmpp+jj-1.AND. i3 == narea ) THEN141 IF( i2 == mjg(jj) .AND. i3 == narea ) THEN 143 142 IF( nicbdj < 0 ) THEN ; nicbdj = jj 144 143 ELSE ; nicbej = jj … … 300 299 rn_test_box(3) < gphit(ji,jj) .AND. gphit(ji,jj) < rn_test_box(4) ) THEN 301 300 localberg%mass_scaling = rn_mass_scaling(iberg) 302 localpt%xi = REAL( nimpp+ji-1, wp )303 localpt%yj = REAL( njmpp+jj-1, wp )301 localpt%xi = REAL( mig(ji), wp ) 302 localpt%yj = REAL( mjg(jj), wp ) 304 303 localpt%lon = icb_utl_bilin(glamt, localpt%xi, localpt%yj, 'T', 0, 0 ) 305 304 localpt%lat = icb_utl_bilin(gphit, localpt%xi, localpt%yj, 'T', 0, 0 )
Note: See TracChangeset
for help on using the changeset viewer.