Changeset 6589 for branches/UKMO/dev_r5518_sshinc_with_VVL
- Timestamp:
- 2016-05-20T14:29:30+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_sshinc_with_VVL/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r6526 r6589 96 96 IF( neuler == 0 .AND. kt == nit000 ) z2dt = rdt 97 97 98 99 #if defined key_asminc 100 ! ! Include the IAU weighted SSH increment 101 IF( lk_asminc .AND. ln_sshinc .AND. ln_asmiau ) THEN 102 CALL ssh_asm_inc( kt ) 103 #if defined key_vvl 104 ! Don't directly adjust ssh but change hdivn at all levels instead 105 ! In trasbc also add in the heat and salt content associated with these changes at each level 106 DO jk = 1, jpkm1 107 hdivn(:,:,jk) = hdivn(:,:,jk) - ( ssh_iau(:,:) / fse3t_n(:,:,jk) ) * ( e3t_0(:,:,jk) / ht_0(:,:) ) * tmask(:,:,jk) 108 END DO 109 ENDIF 110 #endif 111 #endif 112 113 98 114 ! !------------------------------! 99 115 ! ! After Sea Surface Height ! … … 124 140 #endif 125 141 126 #if defined key_asminc127 ! ! Include the IAU weighted SSH increment128 IF( lk_asminc .AND. ln_sshinc .AND. ln_asmiau ) THEN129 CALL ssh_asm_inc( kt )130 ssha(:,:) = ssha(:,:) + z2dt * ssh_iau(:,:)131 #if defined key_vvl132 ! Also add heat and salt fluxes associated with this volume change to keep heat and salt content unchanged.133 qns(:,:) = qns(:,:) + ssh_iau(:,:) * rau0_rcp * tsn(:,:,1,jp_tem)134 sfx(:,:) = sfx(:,:) + ssh_iau(:,:) * rau0 * tsn(:,:,1,jp_sal)135 ENDIF136 #endif137 #endif138 142 139 143 ! !------------------------------!
Note: See TracChangeset
for help on using the changeset viewer.