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 1129 for trunk/NEMO/OPA_SRC/TRD/trdicp.F90 – NEMO

Ignore:
Timestamp:
2008-06-24T08:45:37+02:00 (16 years ago)
Author:
ctlod
Message:

trunk: compilation error when key_trddyn is active related to the dynamics advection flux formulation, see ticket: #211

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/TRD/trdicp.F90

    r719 r1129  
    472472            WRITE (numout,9504) umo(jpicpd_pvo) / tvolu, vmo(jpicpd_pvo) / tvolv 
    473473            WRITE (numout,9505) umo(jpicpd_ldf) / tvolu, vmo(jpicpd_ldf) / tvolv 
    474             WRITE (numout,9506) umo(jpicpd_zad) / tvolu, vmo(jpicpd_zad) / tvolv 
    475             WRITE (numout,9507) umo(jpicpd_zdf) / tvolu, vmo(jpicpd_zdf) / tvolv 
    476             WRITE (numout,9508) umo(jpicpd_spg) / tvolu, vmo(jpicpd_spg) / tvolv 
    477             WRITE (numout,9509) umo(jpicpd_swf) / tvolu, vmo(jpicpd_swf) / tvolv 
    478             WRITE (numout,9510) umo(jpicpd_dat) / tvolu, vmo(jpicpd_dat) / tvolv 
    479             WRITE (numout,9511) umo(jpicpd_bfr) / tvolu, vmo(jpicpd_bfr) / tvolv 
    480             WRITE (numout,9512) 
    481             WRITE (numout,9513)                                                 & 
     474            WRITE (numout,9506) umo(jpicpd_had) / tvolu, vmo(jpicpd_had) / tvolv 
     475            WRITE (numout,9507) umo(jpicpd_zad) / tvolu, vmo(jpicpd_zad) / tvolv 
     476            WRITE (numout,9508) umo(jpicpd_zdf) / tvolu, vmo(jpicpd_zdf) / tvolv 
     477            WRITE (numout,9509) umo(jpicpd_spg) / tvolu, vmo(jpicpd_spg) / tvolv 
     478            WRITE (numout,9510) umo(jpicpd_swf) / tvolu, vmo(jpicpd_swf) / tvolv 
     479            WRITE (numout,9511) umo(jpicpd_dat) / tvolu, vmo(jpicpd_dat) / tvolv 
     480            WRITE (numout,9512) umo(jpicpd_bfr) / tvolu, vmo(jpicpd_bfr) / tvolv 
     481            WRITE (numout,9513) 
     482            WRITE (numout,9514)                                                 & 
    482483            &     (  umo(jpicpd_hpg) + umo(jpicpd_keg) + umo(jpicpd_rvo) + umo(jpicpd_pvo) + umo(jpicpd_ldf)   & 
    483             &      + umo(jpicpd_zad) + umo(jpicpd_zdf) + umo(jpicpd_spg) + umo(jpicpd_dat) + umo(jpicpd_swf)   & 
    484             &      + umo(jpicpd_bfr) ) / tvolu,   & 
     484            &      + umo(jpicpd_had) + umo(jpicpd_zad) + umo(jpicpd_zdf) + umo(jpicpd_spg) + umo(jpicpd_dat)   & 
     485            &      + umo(jpicpd_swf) + umo(jpicpd_bfr) ) / tvolu,   & 
    485486            &     (  vmo(jpicpd_hpg) + vmo(jpicpd_keg) + vmo(jpicpd_rvo) + vmo(jpicpd_pvo) + vmo(jpicpd_ldf)   & 
    486             &      + vmo(jpicpd_zad) + vmo(jpicpd_zdf) + vmo(jpicpd_spg) + vmo(jpicpd_dat) + vmo(jpicpd_swf)   & 
    487             &      + vmo(jpicpd_bfr) ) / tvolv 
     487            &      + vmo(jpicpd_had) + vmo(jpicpd_zad) + vmo(jpicpd_zdf) + vmo(jpicpd_spg) + vmo(jpicpd_dat)   & 
     488            &      + vmo(jpicpd_swf) + vmo(jpicpd_bfr) ) / tvolv 
    488489         ENDIF 
    489490 
    … …  
    494495 9504    FORMAT(' coriolis term              u= ', e20.13, '    v= ', e20.13) 
    495496 9505    FORMAT(' horizontal diffusion       u= ', e20.13, '    v= ', e20.13) 
    496  9506    FORMAT(' vertical advection         u= ', e20.13, '    v= ', e20.13) 
    497  9507    FORMAT(' vertical diffusion         u= ', e20.13, '    v= ', e20.13) 
    498  9508    FORMAT(' surface pressure gradient  u= ', e20.13, '    v= ', e20.13) 
    499  9509    FORMAT(' surface wind forcing       u= ', e20.13, '    v= ', e20.13) 
    500  9510    FORMAT(' dampimg term               u= ', e20.13, '    v= ', e20.13) 
    501  9511    FORMAT(' bottom flux                u= ', e20.13, '    v= ', e20.13) 
    502  9512    FORMAT(' -----------------------------------------------------------------------------') 
    503  9513    FORMAT(' total trend                u= ', e20.13, '    v= ', e20.13) 
     497 9506    FORMAT(' horizontal advection       u= ', e20.13, '    v= ', e20.13) 
     498 9507    FORMAT(' vertical advection         u= ', e20.13, '    v= ', e20.13) 
     499 9508    FORMAT(' vertical diffusion         u= ', e20.13, '    v= ', e20.13) 
     500 9509    FORMAT(' surface pressure gradient  u= ', e20.13, '    v= ', e20.13) 
     501 9510    FORMAT(' surface wind forcing       u= ', e20.13, '    v= ', e20.13) 
     502 9511    FORMAT(' dampimg term               u= ', e20.13, '    v= ', e20.13) 
     503 9512    FORMAT(' bottom flux                u= ', e20.13, '    v= ', e20.13) 
     504 9513    FORMAT(' -----------------------------------------------------------------------------') 
     505 9514    FORMAT(' total trend                u= ', e20.13, '    v= ', e20.13) 
    504506 
    505507         IF(lwp) THEN 
    … …  
    512514            WRITE (numout,9524) hke(jpicpd_pvo) / tvolt 
    513515            WRITE (numout,9525) hke(jpicpd_ldf) / tvolt 
    514             WRITE (numout,9526) hke(jpicpd_zad) / tvolt 
    515             WRITE (numout,9527) hke(jpicpd_zdf) / tvolt 
    516             WRITE (numout,9528) hke(jpicpd_spg) / tvolt 
    517             WRITE (numout,9529) hke(jpicpd_swf) / tvolt 
    518             WRITE (numout,9530) hke(jpicpd_dat) / tvolt 
    519             WRITE (numout,9531) 
    520             WRITE (numout,9532)   & 
     516            WRITE (numout,9526) hke(jpicpd_had) / tvolt 
     517            WRITE (numout,9527) hke(jpicpd_zad) / tvolt 
     518            WRITE (numout,9528) hke(jpicpd_zdf) / tvolt 
     519            WRITE (numout,9529) hke(jpicpd_spg) / tvolt 
     520            WRITE (numout,9530) hke(jpicpd_swf) / tvolt 
     521            WRITE (numout,9531) hke(jpicpd_dat) / tvolt 
     522            WRITE (numout,9532) 
     523            WRITE (numout,9533)   & 
    521524            &     (  hke(jpicpd_hpg) + hke(jpicpd_keg) + hke(jpicpd_rvo) + hke(jpicpd_pvo) + hke(jpicpd_ldf)   & 
    522             &      + hke(jpicpd_zad) + hke(jpicpd_zdf) + hke(jpicpd_spg) + hke(jpicpd_dat) + hke(jpicpd_swf) ) / tvolt 
     525            &      + hke(jpicpd_had) + hke(jpicpd_zad) + hke(jpicpd_zdf) + hke(jpicpd_spg) + hke(jpicpd_dat)   & 
     526            &      + hke(jpicpd_swf) ) / tvolt 
    523527         ENDIF 
    524528 
    … …  
    529533 9524    FORMAT(' coriolis term             u2= ', e20.13) 
    530534 9525    FORMAT(' horizontal diffusion      u2= ', e20.13) 
    531  9526    FORMAT(' vertical advection        u2= ', e20.13) 
    532  9527    FORMAT(' vertical diffusion        u2= ', e20.13) 
    533  9528    FORMAT(' surface pressure gradient u2= ', e20.13) 
    534  9529    FORMAT(' surface wind forcing      u2= ', e20.13) 
    535  9530    FORMAT(' dampimg term              u2= ', e20.13) 
    536  9531    FORMAT(' --------------------------------------------------') 
    537  9532    FORMAT(' total trend               u2= ', e20.13) 
     535 9526    FORMAT(' horizontal advection      u2= ', e20.13) 
     536 9527    FORMAT(' vertical advection        u2= ', e20.13) 
     537 9528    FORMAT(' vertical diffusion        u2= ', e20.13) 
     538 9529    FORMAT(' surface pressure gradient u2= ', e20.13) 
     539 9530    FORMAT(' surface wind forcing      u2= ', e20.13) 
     540 9531    FORMAT(' dampimg term              u2= ', e20.13) 
     541 9532    FORMAT(' --------------------------------------------------') 
     542 9533    FORMAT(' total trend               u2= ', e20.13) 
    538543 
    539544         IF(lwp) THEN 
    … …  
    541546            WRITE (numout,*) 
    542547            WRITE (numout,9540) kt 
    543             WRITE (numout,9541) ( hke(jpicpd_keg) + hke(jpicpd_rvo) + hke(jpicpd_zad) ) / tvolt 
    544             WRITE (numout,9542) ( hke(jpicpd_keg) + hke(jpicpd_zad) ) / tvolt 
     548            WRITE (numout,9541) ( hke(jpicpd_keg) + hke(jpicpd_rvo) + hke(jpicpd_had) + hke(jpicpd_zad) ) / tvolt 
     549            WRITE (numout,9542) ( hke(jpicpd_keg) + hke(jpicpd_had) + hke(jpicpd_zad) ) / tvolt 
    545550            WRITE (numout,9543) ( hke(jpicpd_pvo) ) / tvolt 
    546551            WRITE (numout,9544) ( hke(jpicpd_rvo) ) / tvolt 
    … …  
    555560 9540    FORMAT(' energetic consistency at it= ', i6, ' :', /' =========================================') 
    556561 9541    FORMAT(' 0 = non linear term(true if key_vorenergy or key_combined): ', e20.13) 
    557  9542    FORMAT(' 0 = ke gradient + vertical advection                      : ', e20.13) 
     562 9542    FORMAT(' 0 = ke gradient + horizontal + vertical advection         : ', e20.13) 
    558563 9543    FORMAT(' 0 = coriolis term  (true if key_vorenergy or key_combined): ', e20.13) 
    559564 9544    FORMAT(' 0 = uh.( rot(u) x uh ) (true if enstrophy conser.)        : ', e20.13) 
Note: See TracChangeset for help on using the changeset viewer.