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 49 for trunk/NEMO/OPA_SRC/eosbn2.F90 – NEMO

Ignore:
Timestamp:
2004-04-22T11:13:02+02:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE023 : Addition of new diagnostics controled with logical key l_ctl

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/eosbn2.F90

    r15 r49  
    222222      END SELECT 
    223223 
     224      IF(l_ctl)   WRITE(numout,*) ' eos  : ', SUM( prd(:,:,:) ) 
     225 
    224226   END SUBROUTINE eos_insitu 
    225227 
     
    371373                  zs = psal(ji,jj,jk) 
    372374                  !   ... density and potential volumic mass 
    373                   prd(ji,jj,jk) = ( 0.0285 - ralpha * zt ) * tmask(ji,jj,jk) 
     375                  prd  (ji,jj,jk) = ( 0.0285 - ralpha * zt )        * tmask(ji,jj,jk) 
    374376                  prhop(ji,jj,jk) = ( rau0 * prd(ji,jj,jk) + rau0 ) * tmask(ji,jj,jk) 
    375377               END DO 
     
    390392                  zs = psal(ji,jj,jk) 
    391393                  !   ... density and potential volumic mass 
    392                   prd(ji,jj,jk) = (  rbeta  * zs - ralpha * zt - 1. ) * tmask(ji,jj,jk) 
    393                   prhop(ji,jj,jk) = ( rau0 * prd(ji,jj,jk) + rau0 ) * tmask(ji,jj,jk) 
     394                  prd  (ji,jj,jk) = ( rbeta  * zs - ralpha * zt - 1. ) * tmask(ji,jj,jk) 
     395                  prhop(ji,jj,jk) = ( rau0 * prd(ji,jj,jk) + rau0 )    * tmask(ji,jj,jk) 
    394396               END DO 
    395397            END DO 
     
    405407 
    406408      END SELECT 
     409 
     410      IF(l_ctl)   WRITE(numout,*) ' eos-p: ', SUM( prd(:,:,:) ), ' pot : ', SUM( prhop(:,:,:) ) 
    407411 
    408412   END SUBROUTINE eos_insitu_pot 
     
    580584      END SELECT 
    581585 
     586      IF(l_ctl)   WRITE(numout,*) ' eos2d: ', SUM( prd(1:jpim1,1:jpjm1) ) 
     587 
    582588   END SUBROUTINE eos_insitu_2d 
    583589 
     
    753759      END SELECT 
    754760 
     761      IF(l_ctl)   WRITE(numout,*) ' bn2  : ', SUM( pn2 (:,:,:) ) 
     762#if defined key_zdfddm 
     763      IF(l_ctl)   WRITE(numout,*) ' rrau : ', SUM( rrau(:,:,:) ) 
     764#endif 
     765 
    755766   END SUBROUTINE eos_bn2 
    756767 
Note: See TracChangeset for help on using the changeset viewer.