Changeset 6754
- Timestamp:
- 2016-06-29T11:15:45+02:00 (9 years ago)
- Location:
- branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r6626 r6754 74 74 INTEGER, INTENT(in) :: kt ! time step 75 75 ! 76 INTEGER :: j k ! dummy loop indice76 INTEGER :: ji, jj, jk ! dummy loop indices 77 77 REAL(wp) :: z2dt, z1_rau0 ! local scalars 78 78 !!---------------------------------------------------------------------- … … 94 94 z2dt = 2._wp * rdt ! set time step size (Euler/Leapfrog) 95 95 IF( neuler == 0 .AND. kt == nit000 ) z2dt = rdt 96 97 #if defined key_asminc 98 ! ! Include the IAU weighted SSH increment 99 IF( lk_asminc .AND. ln_sshinc .AND. ln_asmiau ) THEN 100 CALL ssh_asm_inc( kt ) 101 #if defined key_vvl 102 ! Don't directly adjust ssh but change hdivn at all levels instead 103 ! In trasbc also add in the heat and salt content associated with these changes at each level 104 DO jk = 1, jpkm1 105 hdivn(:,:,jk) = hdivn(:,:,jk) - ( ssh_iau(:,:) / ( ht_0(:,:) + 1.0 - ssmask(:,:) ) ) * ( e3t_0(:,:,jk) / fse3t_n(:,:,jk) ) * tmask(:,:,jk) 106 END DO 107 CALL lbc_lnk( hdivn, 'T', 1. ) ! Not sure that's necessary 108 ENDIF 109 #endif 110 #endif 96 111 97 112 ! !------------------------------! … … 123 138 #endif 124 139 125 #if defined key_asminc126 ! ! Include the IAU weighted SSH increment127 IF( lk_asminc .AND. ln_sshinc .AND. ln_asmiau ) THEN128 CALL ssh_asm_inc( kt )129 ssha(:,:) = ssha(:,:) + z2dt * ssh_iau(:,:)130 ENDIF131 #endif132 133 140 ! !------------------------------! 134 141 ! ! outputs ! -
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r6626 r6754 33 33 USE timing ! Timing 34 34 USE eosbn2 35 #if defined key_asminc 36 USE asminc ! Assimilation increment 37 #endif 35 38 36 39 IMPLICIT NONE … … 278 281 END DO 279 282 ENDIF 283 284 #if defined key_asminc 285 ! WARNING: THIS MAY WELL NOT BE REQUIRED - WE DON'T WANT TO CHANGE T&S BUT THIS MAY COMPENSATE ANOTHER TERM... 286 ! Rate of change in e3t for each level is ssh_iau*e3t_0/ht_0 287 ! Contribution to tsa should be rate of change in level / per m of ocean? (hence the division by fse3t_n) 288 IF( ln_sshinc ) THEN ! input of heat and salt due to assimilation 289 DO jj = 2, jpj 290 DO ji = fs_2, fs_jpim1 291 zdep = ssh_iau(ji,jj) / ( ht_0(ji,jj) + 1.0 - ssmask(ji, jj) ) 292 DO jk = 1, jpkm1 293 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 294 & + tsn(ji,jj,jk,jp_tem) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) 295 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) & 296 & + tsn(ji,jj,jk,jp_sal) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) 297 END DO 298 END DO 299 END DO 300 ENDIF 301 #endif 280 302 281 303 IF( l_trdtra ) THEN ! send trends for further diagnostics
Note: See TracChangeset
for help on using the changeset viewer.