Ignore:
Timestamp:
2020-12-10T18:32:40+01:00 (12 months ago)
Author:
kingr
Message:

Removed mass/salt conservation from SSH IAU code.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_sshiau/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r10302 r14154  
    291291      IF( iom_use('rnf_x_sss') )   CALL iom_put( "rnf_x_sss", rnf*tsn(:,:,1,jp_sal) )   ! runoff term on sss 
    292292 
    293 #if defined key_asminc 
    294 ! WARNING: THIS MAY WELL NOT BE REQUIRED - WE DON'T WANT TO CHANGE T&S BUT THIS MAY COMPENSATE ANOTHER TERM... 
    295 ! Rate of change in e3t for each level is ssh_iau*e3t_0/ht_0 
    296 ! Contribution to tsa should be rate of change in level / per m of ocean? (hence the division by fse3t_n) 
    297       IF( ln_sshinc ) THEN         ! input of heat and salt due to assimilation 
    298          DO jj = 2, jpj  
    299             DO ji = fs_2, fs_jpim1 
    300                zdep = ssh_iau(ji,jj) / ( ht_0(ji,jj) + 1.0 - ssmask(ji, jj) ) 
    301                DO jk = 1, jpkm1 
    302                   tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem)   & 
    303                                         &            + tsn(ji,jj,jk,jp_tem) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) 
    304                   tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal)   & 
    305                                         &            + tsn(ji,jj,jk,jp_sal) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) 
    306                END DO 
    307             END DO   
    308          END DO   
    309       ENDIF 
    310 #endif 
    311   
    312293      IF( l_trdtra )   THEN                      ! send trends for further diagnostics 
    313294         ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.