New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3374 for branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90 – NEMO

Ignore:
Timestamp:
2012-04-30T13:24:25+02:00 (12 years ago)
Author:
sga
Message:

NEMO branch dev_r3337_NOCS10_ICB: adopt mig, mjg, mi1 and mj1 for improved obscurity

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90

    r3372 r3374  
    106106         DO ji = 1, jpi 
    107107            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) 
    110109         END DO 
    111110      END DO 
     
    124123         i1 = i3 - i2*nicbpack 
    125124         i3 = INT( src_calving_hflx(ji,jj) ) 
    126          IF( i1 == nimpp+ji-1 .AND. i3 == narea ) THEN 
     125         IF( i1 == mig(ji) .AND. i3 == narea ) THEN 
    127126            IF( nicbdi < 0 ) THEN   ;   nicbdi = ji 
    128127            ELSE                    ;   nicbei = ji 
     
    140139         i1 = i3 - i2*nicbpack 
    141140         i3 = INT( src_calving_hflx(ji,jj) ) 
    142          IF( i2 == njmpp+jj-1 .AND. i3 == narea ) THEN 
     141         IF( i2 == mjg(jj) .AND. i3 == narea ) THEN 
    143142            IF( nicbdj < 0 ) THEN   ;   nicbdj = jj 
    144143            ELSE                    ;   nicbej = jj 
     
    300299                rn_test_box(3) < gphit(ji,jj) .AND. gphit(ji,jj) < rn_test_box(4) ) THEN 
    301300               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 ) 
    304303               localpt%lon = icb_utl_bilin(glamt, localpt%xi, localpt%yj, 'T', 0, 0 ) 
    305304               localpt%lat = icb_utl_bilin(gphit, localpt%xi, localpt%yj, 'T', 0, 0 ) 
Note: See TracChangeset for help on using the changeset viewer.