Ignore:
Timestamp:
2019-10-18T16:49:33+02:00 (12 months ago)
Author:
jpalmier
Message:

add Changes to enable Alkalinity geo-ingeneering tests within MEDUSA

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_GeoAlk_MEDUSA/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

    r10149 r11726  
    173173                  fq2 = fq0 + fq1 
    174174                  fq3 = f_benout_n(ji,jj) 
    175                   if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     175                  if (lwp) write (numout,'(a,2i3,a,4f15.5)')                   & 
    176176                     'AXY N   cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',      & 
    177177                     fq0,fq1,fq2,fq3 
     
    187187                  fq2 = fq0 + fq1 
    188188                  fq3 = f_benout_si(ji,jj) 
    189                   if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     189                  if (lwp) write (numout,'(a,2i3,a,4f15.5)')                   & 
    190190                     'AXY Si  cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',     & 
    191191                     fq0,fq1,fq2,fq3 
     
    202202                  fq3 = fq0 + fq1 
    203203                  fq4 = f_benout_c(ji,jj) + f_benout_ca(ji,jj) 
    204                   if (lwp) write (numout,'a,2i3,a,5f15,5)')                   & 
     204                  if (lwp) write (numout,'(a,2i3,a,5f15.5)')                   & 
    205205                     'AXY C   cons: (i,j)=',ji,jj,', (flx,ben,asf,err,out)=', & 
    206206                     fq0,fq1,fq2,fq3,fq4 
     
    216216                   fq2 = fq0 + fq1 
    217217                   fq3 = 2.0 * f_benout_ca(ji,jj) 
    218                    if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     218                   if (lwp) write (numout,'(a,2i3,a,4f15.5)')                   & 
    219219                      'AXY alk cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',     & 
    220220                      fq0,fq1,fq2,fq3 
     
    593593            DEALLOCATE( rivalk2d ) 
    594594         ENDIF 
     595         !! AXY (16/10/19) 
     596         IF( med_diag%GEOENG_ALK%dgsave ) THEN 
     597            CALL iom_put( "GEOENG_ALK"  , geoengalk2d ) 
     598            DEALLOCATE( geoengalk2d ) 
     599         ENDIF 
    595600         IF( med_diag%DETC%dgsave ) THEN 
    596601            CALL iom_put( "DETC"  , fslowc2d ) 
Note: See TracChangeset for help on using the changeset viewer.