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 10196 for branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90 – NEMO

Ignore:
Timestamp:
2018-10-16T12:15:14+02:00 (6 years ago)
Author:
jpalmier
Message:

add DMS flux --

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

    r10047 r10196  
    3333      !!---------------------------------------------------------------------- 
    3434      USE bio_medusa_mod 
    35       USE dom_oce,           ONLY: atfp, atfp1, neuler, rdt, tmask 
     35      USE dom_oce,           ONLY: atfp, atfp1, neuler, rdt, tmask, e3t_0 
     36# if defined key_vvl 
     37      USE dom_oce,           ONLY: e3t_n  
     38# endif 
    3639      USE in_out_manager,    ONLY: lwp, numout 
    3740      USE iom,               ONLY: iom_put 
    … …  
    6265      INTEGER :: jn 
    6366 
    64       REAL(wp) :: fq0,fq1,fq2,fq3 
     67      REAL(wp) :: fq0,fq1,fq2,fq3,fq4 
     68 
     69 
     70   !!* Substitution 
     71#  include "domzgr_substitute.h90" 
    6572 
    6673# if defined key_roam                      
    … …  
    149156                  fq2 = fq0 + fq1 
    150157                  fq3 = f_benout_n(ji,jj) 
    151                   if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     158                  if (lwp) write (numout,'(a,2i3,a,4f15.5)')                   & 
    152159                     'AXY N   cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',      & 
    153160                     fq0,fq1,fq2,fq3 
    … …  
    163170                  fq2 = fq0 + fq1 
    164171                  fq3 = f_benout_si(ji,jj) 
    165                   if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     172                  if (lwp) write (numout,'(a,2i3,a,4f15.5)')                   & 
    166173                     'AXY Si  cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',     & 
    167174                     fq0,fq1,fq2,fq3 
    … …  
    178185                  fq3 = fq0 + fq1 
    179186                  fq4 = f_benout_c(ji,jj) + f_benout_ca(ji,jj) 
    180                   if (lwp) write (numout,'a,2i3,a,5f15,5)')                   & 
     187                  if (lwp) write (numout,'(a,2i3,a,5f15.5)')                   & 
    181188                     'AXY C   cons: (i,j)=',ji,jj,', (flx,ben,asf,err,out)=', & 
    182189                     fq0,fq1,fq2,fq3,fq4 
    … …  
    192199                   fq2 = fq0 + fq1 
    193200                   fq3 = 2.0 * f_benout_ca(ji,jj) 
    194                    if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     201                   if (lwp) write (numout,'(a,2i3,a,4f15.5)')                   & 
    195202                      'AXY alk cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',     & 
    196203                      fq0,fq1,fq2,fq3 
Note: See TracChangeset for help on using the changeset viewer.