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 3981 for branches/2013/dev_r3406_CNRS_LIM3 – NEMO

Ignore:
Timestamp:
2013-07-17T14:36:59+02:00 (11 years ago)
Author:
clem
Message:

bug correction

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90

    r3963 r3981  
    320320      REAL(wp)   ::   ztmp 
    321321      INTEGER    ::   ji, jj   ! dummy loop indices 
    322       INTEGER    ::   zjpi, zjpj ! local variables: size of ptab 
    323       !!----------------------------------------------------------------------- 
    324       zjpi = SIZE(ptab,1) 
    325       zjpj = SIZE(ptab,2) 
    326       ! 
    327       ztmp = 0.e0 
    328       ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    329       DO jj = 1, zjpj 
    330          DO ji =1, zjpi 
    331          ztmp =  ptab(ji,jj) * tmask_i(ji,jj) 
    332          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     322      !!----------------------------------------------------------------------- 
     323      ! 
     324      ztmp = 0.e0 
     325      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     326      DO jj = 1, jpj 
     327         DO ji = 1, jpi 
     328            ztmp =  ptab(ji,jj) * tmask_i(ji,jj) 
     329            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    333330         END DO 
    334331      END DO 
    … …  
    351348      REAL(wp)   ::   ztmp 
    352349      INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    353       INTEGER    ::   zjpi, zjpj, zjpk ! local variables: size of ptab 
    354       !!----------------------------------------------------------------------- 
    355       ! 
    356       zjpi = SIZE(ptab,1) 
    357       zjpj = SIZE(ptab,2) 
     350      INTEGER    ::   zjpk ! local variables: size of ptab 
     351      !!----------------------------------------------------------------------- 
     352      ! 
    358353      zjpk = SIZE(ptab,3) 
    359354      ! 
    … …  
    361356      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    362357      DO jk = 1, zjpk 
    363          DO jj = 1, zjpj 
    364             DO ji =1, zjpi 
    365             ztmp =  ptab(ji,jj,jk) * tmask_i(ji,jj) 
    366             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     358         DO jj = 1, jpj 
     359            DO ji = 1, jpi 
     360               ztmp =  ptab(ji,jj,jk) * tmask_i(ji,jj) 
     361               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    367362            END DO 
    368363         END DO     
    … …  
    386381      REAL(wp)   ::   ztmp 
    387382      INTEGER    ::   ji, jj   ! dummy loop indices 
    388       INTEGER    ::   zjpi, zjpj ! local variables: size of ptab 
    389       !!----------------------------------------------------------------------- 
    390       ! 
    391       zjpi = SIZE(ptab1,1) 
    392       zjpj = SIZE(ptab1,2) 
    393       ! 
    394       ztmp = 0.e0 
    395       ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    396       DO jj = 1, zjpj 
    397          DO ji =1, zjpi 
    398          ztmp =  ptab1(ji,jj) * tmask_i(ji,jj) 
    399          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    400          ztmp =  ptab2(ji,jj) * tmask_i(ji,jj) 
    401          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     383      !!----------------------------------------------------------------------- 
     384      ! 
     385      ztmp = 0.e0 
     386      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     387      DO jj = 1, jpj 
     388         DO ji = 1, jpi 
     389            ztmp =  ptab1(ji,jj) * tmask_i(ji,jj) 
     390            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     391            ztmp =  ptab2(ji,jj) * tmask_i(ji,jj) 
     392            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    402393         END DO 
    403394      END DO 
    … …  
    420411      REAL(wp)   ::   ztmp 
    421412      INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    422       INTEGER    ::   zjpi, zjpj, zjpk ! local variables: size of ptab 
    423       !!----------------------------------------------------------------------- 
    424       ! 
    425       zjpi = SIZE(ptab1,1) 
    426       zjpj = SIZE(ptab1,2) 
     413      INTEGER    ::   zjpk ! local variables: size of ptab 
     414      !!----------------------------------------------------------------------- 
     415      ! 
    427416      zjpk = SIZE(ptab1,3) 
    428417      ! 
    … …  
    430419      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    431420      DO jk = 1, zjpk 
    432          DO jj = 1, zjpj 
    433             DO ji =1, zjpi 
    434             ztmp =  ptab1(ji,jj,jk) * tmask_i(ji,jj) 
    435             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    436             ztmp =  ptab2(ji,jj,jk) * tmask_i(ji,jj) 
    437             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     421         DO jj = 1, jpj 
     422            DO ji = 1, jpi 
     423               ztmp =  ptab1(ji,jj,jk) * tmask_i(ji,jj) 
     424               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     425               ztmp =  ptab2(ji,jj,jk) * tmask_i(ji,jj) 
     426               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    438427            END DO 
    439428         END DO     
    … …  
    458447      REAL(wp)   ::   ztmp 
    459448      INTEGER    ::   ji, jj   ! dummy loop indices 
    460       INTEGER    ::   zjpi, zjpj ! local variables: size of ptab 
    461       !!----------------------------------------------------------------------- 
    462       zjpi = SIZE(ptab,1) 
    463       zjpj = SIZE(ptab,2) 
    464       ! 
    465       ztmp = 0.e0 
    466       ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    467       DO jj = 1, zjpj 
    468          DO ji =1, zjpi 
    469          ztmp =  ptab(ji,jj) * tmask_i(ji,jj) 
    470          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     449      !!----------------------------------------------------------------------- 
     450      ! 
     451      ztmp = 0.e0 
     452      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     453      DO jj = 1, jpj 
     454         DO ji = 1, jpi 
     455            ztmp =  ptab(ji,jj) * tmask_i(ji,jj) 
     456            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    471457         END DO 
    472458      END DO 
    … …  
    489475      REAL(wp)   ::   ztmp 
    490476      INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    491       INTEGER    ::   zjpi, zjpj, zjpk ! local variables: size of ptab 
    492       !!----------------------------------------------------------------------- 
    493       ! 
    494       zjpi = SIZE(ptab,1) 
    495       zjpj = SIZE(ptab,2) 
     477      INTEGER    ::   zjpk ! local variables: size of ptab 
     478      !!----------------------------------------------------------------------- 
     479      ! 
    496480      zjpk = SIZE(ptab,3) 
    497481      ! 
    … …  
    499483      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    500484      DO jk = 1, zjpk 
    501          DO jj = 1, zjpj 
    502             DO ji =1, zjpi 
    503             ztmp =  ptab(ji,jj,jk) * tmask_i(ji,jj) 
    504             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     485         DO jj = 1, jpj 
     486            DO ji = 1, jpi 
     487               ztmp =  ptab(ji,jj,jk) * tmask_i(ji,jj) 
     488               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    505489            END DO 
    506490         END DO     
    … …  
    524508      REAL(wp)   ::   ztmp 
    525509      INTEGER    ::   ji, jj   ! dummy loop indices 
    526       INTEGER    ::   zjpi, zjpj ! local variables: size of ptab 
    527       !!----------------------------------------------------------------------- 
    528       ! 
    529       zjpi = SIZE(ptab1,1) 
    530       zjpj = SIZE(ptab1,2) 
    531       ! 
    532       ztmp = 0.e0 
    533       ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    534       DO jj = 1, zjpj 
    535          DO ji =1, zjpi 
    536          ztmp =  ptab1(ji,jj) * tmask_i(ji,jj) 
    537          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    538          ztmp =  ptab2(ji,jj) * tmask_i(ji,jj) 
    539          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     510      !!----------------------------------------------------------------------- 
     511      ! 
     512      ! 
     513      ztmp = 0.e0 
     514      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     515      DO jj = 1, jpj 
     516         DO ji = 1, jpi 
     517            ztmp =  ptab1(ji,jj) * tmask_i(ji,jj) 
     518            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     519            ztmp =  ptab2(ji,jj) * tmask_i(ji,jj) 
     520            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    540521         END DO 
    541522      END DO 
    … …  
    558539      REAL(wp)   ::   ztmp 
    559540      INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    560       INTEGER    ::   zjpi, zjpj, zjpk ! local variables: size of ptab 
    561       !!----------------------------------------------------------------------- 
    562       ! 
    563       zjpi = SIZE(ptab1,1) 
    564       zjpj = SIZE(ptab1,2) 
     541      INTEGER    ::   zjpk ! local variables: size of ptab 
     542      !!----------------------------------------------------------------------- 
     543      ! 
    565544      zjpk = SIZE(ptab1,3) 
    566545      ! 
    … …  
    568547      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    569548      DO jk = 1, zjpk 
    570          DO jj = 1, zjpj 
    571             DO ji =1, zjpi 
    572             ztmp =  ptab1(ji,jj,jk) * tmask_i(ji,jj) 
    573             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    574             ztmp =  ptab2(ji,jj,jk) * tmask_i(ji,jj) 
    575             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     549         DO jj = 1, jpj 
     550            DO ji = 1, jpi 
     551               ztmp =  ptab1(ji,jj,jk) * tmask_i(ji,jj) 
     552               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     553               ztmp =  ptab2(ji,jj,jk) * tmask_i(ji,jj) 
     554               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    576555            END DO 
    577556         END DO     
    … …  
    596575      REAL(wp)   ::   ztmp 
    597576      INTEGER    ::   ji, jj   ! dummy loop indices 
    598       INTEGER    ::   zjpi, zjpj ! local variables: size of ptab 
    599       !!----------------------------------------------------------------------- 
    600       zjpi = SIZE(ptab,1) 
    601       zjpj = SIZE(ptab,2) 
    602       ! 
    603       ztmp = 0.e0 
    604       ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    605       DO jj = 1, zjpj 
    606          DO ji =1, zjpi 
    607          ztmp =  ptab(ji,jj) * tmask_i(ji,jj) 
    608          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     577      !!----------------------------------------------------------------------- 
     578      ! 
     579      ztmp = 0.e0 
     580      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     581      DO jj = 1, jpj 
     582         DO ji = 1, jpi 
     583            ztmp =  ptab(ji,jj) * tmask_i(ji,jj) 
     584            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    609585         END DO 
    610586      END DO 
    … …  
    627603      REAL(wp)   ::   ztmp 
    628604      INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    629       INTEGER    ::   zjpi, zjpj, zjpk ! local variables: size of ptab 
    630       !!----------------------------------------------------------------------- 
    631       ! 
    632       zjpi = SIZE(ptab,1) 
    633       zjpj = SIZE(ptab,2) 
     605      INTEGER    ::   zjpk ! local variables: size of ptab 
     606      !!----------------------------------------------------------------------- 
     607      ! 
    634608      zjpk = SIZE(ptab,3) 
    635609      ! 
    … …  
    637611      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    638612      DO jk = 1, zjpk 
    639          DO jj = 1, zjpj 
    640             DO ji =1, zjpi 
    641             ztmp =  ptab(ji,jj,jk) * tmask_i(ji,jj) 
    642             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     613         DO jj = 1, jpj 
     614            DO ji = 1, jpi 
     615               ztmp =  ptab(ji,jj,jk) * tmask_i(ji,jj) 
     616               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    643617            END DO 
    644618         END DO     
    … …  
    662636      REAL(wp)   ::   ztmp 
    663637      INTEGER    ::   ji, jj   ! dummy loop indices 
    664       INTEGER    ::   zjpi, zjpj ! local variables: size of ptab 
    665       !!----------------------------------------------------------------------- 
    666       ! 
    667       zjpi = SIZE(ptab1,1) 
    668       zjpj = SIZE(ptab1,2) 
    669       ! 
    670       ztmp = 0.e0 
    671       ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    672       DO jj = 1, zjpj 
    673          DO ji =1, zjpi 
    674          ztmp =  ptab1(ji,jj) * tmask_i(ji,jj) 
    675          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    676          ztmp =  ptab2(ji,jj) * tmask_i(ji,jj) 
    677          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     638      !!----------------------------------------------------------------------- 
     639      ! 
     640      ! 
     641      ztmp = 0.e0 
     642      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     643      DO jj = 1, jpj 
     644         DO ji = 1, jpi 
     645            ztmp =  ptab1(ji,jj) * tmask_i(ji,jj) 
     646            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     647            ztmp =  ptab2(ji,jj) * tmask_i(ji,jj) 
     648            CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    678649         END DO 
    679650      END DO 
    … …  
    696667      REAL(wp)   ::   ztmp 
    697668      INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    698       INTEGER    ::   zjpi, zjpj, zjpk ! local variables: size of ptab 
    699       !!----------------------------------------------------------------------- 
    700       ! 
    701       zjpi = SIZE(ptab1,1) 
    702       zjpj = SIZE(ptab1,2) 
     669      INTEGER    ::   zjpk ! local variables: size of ptab 
     670      !!----------------------------------------------------------------------- 
     671      ! 
    703672      zjpk = SIZE(ptab1,3) 
    704673      ! 
    … …  
    706675      ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    707676      DO jk = 1, zjpk 
    708          DO jj = 1, zjpj 
    709             DO ji =1, zjpi 
    710             ztmp =  ptab1(ji,jj,jk) * tmask_i(ji,jj) 
    711             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    712             ztmp =  ptab2(ji,jj,jk) * tmask_i(ji,jj) 
    713             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     677         DO jj = 1, jpj 
     678            DO ji = 1, jpi 
     679               ztmp =  ptab1(ji,jj,jk) * tmask_i(ji,jj) 
     680               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     681               ztmp =  ptab2(ji,jj,jk) * tmask_i(ji,jj) 
     682               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    714683            END DO 
    715684         END DO     
Note: See TracChangeset for help on using the changeset viewer.