Changeset 8225


Ignore:
Timestamp:
2017-06-27T15:26:53+02:00 (3 years ago)
Author:
jwhile
Message:

Reintroducing a bug fix as recomended by Tim Graham's review

File:
1 edited

Legend:

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

    r7947 r8225  
    236236               DO jk = ikt, ikb - 1 
    237237               ! compute tfreez for the temperature correction (we add water at freezing temperature) 
    238 !                  zpress = grav*rau0*fsdept(ji,jj,jk)*1.e-04 
    239                   zt_frz = -1.9 !eos_fzp( tsn(ji,jj,jk,jp_sal), zpress ) 
    240238               ! compute trend 
    241239                  tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem)                                          & 
    242                      &           + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)          & 
    243                      &               - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) & 
    244                      &           * r1_hisf_tbl(ji,jj) 
     240                        &           + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) 
    245241                  tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal)                                          & 
    246242                     &           + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) 
     
    249245               ! level partially include in ice shelf boundary layer  
    250246               ! compute tfreez for the temperature correction (we add water at freezing temperature) 
    251 !               zpress = grav*rau0*fsdept(ji,jj,ikb)*1.e-04 
    252                zt_frz = -1.9 !eos_fzp( tsn(ji,jj,ikb,jp_sal), zpress ) 
    253247               ! compute trend 
    254248               tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem)                                           & 
    255                   &              + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)          & 
    256                   &                  - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) &  
    257                   &              * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 
     249                     &              + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 
    258250               tsa(ji,jj,ikb,jp_sal) = tsa(ji,jj,ikb,jp_sal)                                           & 
    259251                  &              + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj)  
Note: See TracChangeset for help on using the changeset viewer.