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 4907 for branches/2014/dev_CNRS_2014 – NEMO

Ignore:
Timestamp:
2014-11-27T22:02:56+01:00 (9 years ago)
Author:
cetlod
Message:

2014/dev_CNRS_2014 : minor bugfix

Location:
branches/2014/dev_CNRS_2014/NEMOGCM
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/SHARED/field_def.xml

    r4906 r4907  
    3030         <field id="heatc"        long_name="Heat content vertically integrated"        unit="J/m2"                      /> 
    3131         <field id="saltc"        long_name="Salt content vertically integrated"        unit="PSU*kg/m2"                /> 
    32          <!-- EOS --!> 
     32         <!-- EOS --> 
    3333         <field id="alpha"        long_name="thermal expansion"                         unit="1/degC" grid_ref="grid_T_3D"/> 
    3434         <field id="beta"         long_name="haline contraction"                        unit="1/psu"  grid_ref="grid_T_3D"/> 
    3535         <field id="bn2"          long_name="squared Brunt-Vaisala frequency"           unit="1/s"    grid_ref="grid_T_3D"/> 
    3636         <field id="rhop"         long_name="potential density (sigma0)"                unit="kg/m3"  grid_ref="grid_T_3D"/> 
    37          <!-- Energy - horizontal divergence --!> 
     37         <!-- Energy - horizontal divergence --> 
    3838         <field id="eken"         long_name="kinetic energy"                            unit="m2/s2"  grid_ref="grid_T_3D"/> 
    3939         <field id="hdiv"         long_name="horizontal divergence"                     unit="s-1"    grid_ref="grid_T_3D"/> 
  • branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limdiahsb.F90

    r4902 r4907  
    7575 
    7676      ! 1/area 
    77       z1_area = 1.d0 / MAX( glob_sum( area(:,:) * tms(:,:) ), epsi06 ) 
    78  
    79       zinda = MAX( 0.d0 , SIGN( 1.d0 , glob_sum( area(:,:) * tms(:,:) ) - epsi06 ) ) 
     77      z1_area = 1.0 / MAX( glob_sum( area(:,:) * tms(:,:) ), epsi06 ) 
     78 
     79      zinda = MAX( 0.0 , SIGN( 1.0 , glob_sum( area(:,:) * tms(:,:) ) - epsi06 ) ) 
    8080      ! ----------------------- ! 
    8181      ! 1 -  Content variations ! 
     
    153153      ! 3 - Diagnostics writing ! 
    154154      ! ----------------------- ! 
    155       zindb = MAX( 0.d0 , SIGN( 1.d0 , zbg_ivo - epsi06 ) ) 
     155      zindb = MAX( 0.0 , SIGN( 1.0 , zbg_ivo - epsi06 ) ) 
    156156      ! 
    157157      CALL iom_put( 'ibgvoltot' , zbg_ivo * rhoic * r1_rau0 * 1.e-9        )   ! ice volume (km3 equivalent liquid)          
  • branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

    r4902 r4907  
    222222            !     Second step in limthd_dh      :  heat remaining if total melt (zq_rema)  
    223223            !     Third  step in limsbc         :  heat from ice-ocean mass exchange (zf_mass) + solar 
    224             hfx_out(ji,jj) = hfx_out(ji,jj)                                                                                                        &  
     224            hfx_out(ji,jj) = hfx_out(ji,jj)                                                                                       &  
    225225               ! Non solar heat flux received by the ocean 
    226                &    +        pfrld(ji,jj) * qns(ji,jj)                                                                                             & 
     226               &    +        pfrld(ji,jj) * qns(ji,jj)                                                                            & 
    227227               ! latent heat of precip (note that precip is included in qns but not in qns_ice) 
    228                &    +      ( pfrld(ji,jj)**betas - pfrld(ji,jj) ) * sprecip(ji,jj) * ( cpic * ( MIN( tatm_ice(ji,jj), rt0_snow ) - rtt ) - lfus )  & 
    229                &    +      ( 1._wp - pfrld(ji,jj) ) * ( tprecip(ji,jj) - sprecip(ji,jj) ) * rcp * ( tatm_ice(ji,jj) - rtt )                        & 
     228               &    +      ( pfrld(ji,jj)**betas - pfrld(ji,jj) ) * sprecip(ji,jj)       & 
     229               &         * ( cpic * ( MIN( tatm_ice(ji,jj), rt0_snow ) - rtt ) - lfus )  & 
     230               &    +      ( 1._wp - pfrld(ji,jj) ) * ( tprecip(ji,jj) - sprecip(ji,jj) ) * rcp * ( tatm_ice(ji,jj) - rtt )       & 
    230231               ! heat flux taken from the ocean where there is open water ice formation 
    231                &    -      qlead(ji,jj) * r1_rdtice                                                                                                & 
     232               &    -      qlead(ji,jj) * r1_rdtice                                                                               & 
    232233               ! heat flux taken from the ocean during bottom growth/melt (fhld should be 0 while bott growth) 
    233                &    -      at_i(ji,jj) * fhtur(ji,jj)                                                                                              & 
     234               &    -      at_i(ji,jj) * fhtur(ji,jj)                                                                             & 
    234235               &    -      at_i(ji,jj) *  fhld(ji,jj) 
    235236 
  • branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90

    r4902 r4907  
    133133                  !Energy of melting q(S,T) [J.m-3] 
    134134                  zindb = 1._wp - MAX(  0._wp , SIGN( 1._wp , -v_i(ji,jj,jl) + epsi10 )  )   !0 if no ice and 1 if yes 
    135                   e_i(ji,jj,jk,jl) = zindb * e_i(ji,jj,jk,jl) / ( area(ji,jj) * MAX( v_i(ji,jj,jl) ,  epsi10 ) ) * REAL( nlay_i, wp ) 
     135                  e_i(ji,jj,jk,jl) = zindb * e_i(ji,jj,jk,jl) & 
     136                      &   / ( area(ji,jj) * MAX( v_i(ji,jj,jl) ,  epsi10 ) ) * REAL( nlay_i, wp ) 
    136137                  e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) * unit_fac 
    137138               END DO 
Note: See TracChangeset for help on using the changeset viewer.