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 5895 for branches/UKMO/dev_r5518_flux_adjust/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx_adj.F90 – NEMO

Ignore:
Timestamp:
2015-11-18T17:05:56+01:00 (8 years ago)
Author:
kuniko
Message:

Made treatment of river mouth grids consistent in nn_flxadjfw 1&2.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_flux_adjust/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx_adj.F90

    r5889 r5895  
    105105                  DO ji = 1, jpi 
    106106                     zerp = ( 1. - 2.*rnfmsk(ji,jj) )  &              ! No damping in vicinity of river mouths 
    107                         &        * sf_erp(1)%fnow(ji,jj,1)     &       
    108                         &        * MAX( sss_m(ji,jj), 1.e-20 ) 
     107                        &        * sf_erp(1)%fnow(ji,jj,1) 
    109108                     sfx(ji,jj) = sfx(ji,jj) + zerp                 ! salt flux 
    110109                  END DO 
     
    117116               DO jj = 1, jpj 
    118117                  DO ji = 1, jpi                             
    119                      zerp = sf_erp(1)%fnow(ji,jj,1) 
     118                     zerp = ( 1. - 2.*rnfmsk(ji,jj) )  &              ! No damping in vicinity of river mouths 
     119                        &        * sf_erp(1)%fnow(ji,jj,1) 
    120120                     IF( ln_sssr_bnd )   zerp = SIGN( 1., zerp ) * MIN( zerp_bnd, ABS(zerp) ) 
    121121                     emp(ji,jj) = emp (ji,jj) + zerp 
     
    191191         ! 
    192192         ! fill sf_qrp with sn_qrp and control print 
    193          CALL fld_fill( sf_qrp, (/ sn_qrp /), cn_dir, 'sbc_flx_adj', 'Heat flux adjustment', 'namsbc_flx_adj' ) 
     193         CALL fld_fill( sf_qrp, (/ sn_qrp /), cn_dir, 'sbc_qrp', 'Heat flux adjustment', 'namsbc_flx_adj' ) 
    194194         IF( sf_qrp(1)%ln_tint )   ALLOCATE( sf_qrp(1)%fdta(jpi,jpj,1,2), STAT=ierror ) 
    195195         IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_flx_adj: unable to allocate sf_qrp data array' ) 
Note: See TracChangeset for help on using the changeset viewer.