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 824 – NEMO

Changeset 824


Ignore:
Timestamp:
2008-03-06T12:37:37+01:00 (16 years ago)
Author:
rblod
Message:

Delete LIM_SRC and add corresponding modifications in Ocean part

Location:
branches/dev_002_LIM/NEMO
Files:
1 deleted
17 edited

Legend:

Unmodified
Added
Removed

  • branches/dev_002_LIM/NEMO/NST_SRC/agrif_user.F90

    r782 r824  
    6262      USE sol_oce 
    6363      USE in_out_manager 
    64 #if defined key_ice_lim 
     64#if defined key_ice_lim || defined key_ice_lim_old 
    6565      USE ice_oce 
    6666#endif 

  • branches/dev_002_LIM/NEMO/OPA_SRC/DIA/diafwb.F90

    r719 r824  
    2020   USE blk_oce         ! ??? 
    2121   USE flxblk          ! atmospheric surface quantity 
     22   USE flxblk_2          ! atmospheric surface quantity 
    2223   USE lib_mpp         ! distributed memory computing library 
    2324 

  • branches/dev_002_LIM/NEMO/OPA_SRC/DIA/diawri.F90

    r719 r824  
    245245            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    246246#endif 
    247 #if ! defined key_dynspg_rl && defined key_ice_lim 
     247#if ! defined key_dynspg_rl && ( defined key_ice_lim || defined key_lim2 ) 
    248248         ! sowaflup = sowaflep + sorunoff + sowafldp + a term associated to 
    249249         !    internal damping to Levitus that can be diagnosed from others 
    … …  
    291291 
    292292 
    293 #if ( defined key_coupled && ! defined key_ice_lim )  
     293#if defined key_coupled &&  ! defined key_ice_lim && ! defined key_lim2   
    294294         CALL histdef( nid_T, "sohefldp", "Surface Heat Flux: Damping"         , "W/m2"   ,   &  ! qrp 
    295295            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    … …  
    317317#endif 
    318318 
    319 #if defined key_ice_lim && defined key_coupled 
     319#if ( defined key_ice_lim  || defined key_lim2 ) && defined key_coupled 
    320320         CALL histdef( nid_T,"soicetem" , "Ice Surface Temperature"            , "K"      ,   &  ! tn_ice 
    321321            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    … …  
    423423      CALL histwrite( nid_T, "sossheig", it, sshn          , ndim_hT, ndex_hT )   ! sea surface height 
    424424#endif 
    425 #if ! defined key_dynspg_rl && defined key_ice_lim 
     425#if ! defined key_dynspg_rl && ( defined key_ice_lim || defined key_lim2 ) 
    426426      CALL histwrite( nid_T, "iowaflup", it, fsalt(:,:)    , ndim_hT, ndex_hT )   ! ice=>ocean water flux 
    427427      CALL histwrite( nid_T, "sowaflep", it, fmass(:,:)    , ndim_hT, ndex_hT )   ! atmos=>ocean water flux 
    … …  
    448448      CALL histwrite( nid_T, "sosbhfup", it, qsb           , ndim_hT, ndex_hT )   ! sensible heat flux 
    449449#endif 
    450 #if ( defined key_coupled && ! defined key_ice_lim )  
     450#if  defined key_coupled && ! defined key_ice_lim && ! defined key_lim2  
    451451      CALL histwrite( nid_T, "sohefldp", it, qrp           , ndim_hT, ndex_hT )   ! heat flux damping 
    452452      CALL histwrite( nid_T, "sowafldp", it, erp           , ndim_hT, ndex_hT )   ! freshwater flux damping 
    … …  
    466466      CALL histwrite( nid_T, "sohtc300", it, htc3          , ndim_hT, ndex_hT )   ! first 300m heaat content 
    467467#endif 
    468 #if defined key_ice_lim &&  defined key_coupled  
     468#if ( defined key_ice_lim  ||  defined key_lim2 ) &&  defined key_coupled  
    469469      CALL histwrite( nid_T, "soicetem", it, tn_ice        , ndim_hT, ndex_hT )   ! surf. ice temperature 
    470470      CALL histwrite( nid_T, "soicealb", it, alb_ice       , ndim_hT, ndex_hT )   ! ice albedo 

  • branches/dev_002_LIM/NEMO/OPA_SRC/DIA/diawri_dimg.h90

    r719 r824  
    187187       !        fsel(:,:,15) = fsel(:,:,15) + fbt(:,:) 
    188188       fsel(:,:,16) = fsel(:,:,16) + emps(:,:) 
    189 #if defined key_ice_lim 
     189#if defined key_ice_lim || defined key_ice_lim_old 
    190190       fsel(:,:,17) = fsel(:,:,17) + fsalt(:,:) 
    191191#endif 
    … …  
    277277          !         fsel(:,:,15) =  fbt(:,:) 
    278278          fsel(:,:,16) =  emps(:,:) * tmask(:,:,1) 
    279 #if defined key_ice_lim 
     279#if defined key_ice_lim || defined key_ice_lim_old 
    280280          fsel(:,:,17) =  fsalt(:,:) * tmask(:,:,1) 
    281281#endif 

  • branches/dev_002_LIM/NEMO/OPA_SRC/SBC/bulk.F90

    r719 r824  
    1616   USE ocfzpt          ! ocean freezing point 
    1717   USE flxblk          ! bulk formulae 
     18   USE flxblk_2        ! bulk formulae 
    1819   USE blk_oce         ! bulk variable  
    1920   USE flx_oce 
    … …  
    5051      !! * Local declarations     
    5152      REAL(wp), DIMENSION(jpi,jpj) ::   zsst  
    52 # if ! defined key_ice_lim 
     53# if ( ! defined key_ice_lim && !defined key_lim2 ) 
    5354      INTEGER  ::   ji, jj         ! dummy loop indices   
    5455      REAL(wp) ::   ztgel, zicopa 
    … …  
    6970     ENDIF 
    7071 
    71 # if ! defined key_ice_lim 
     72# if ( ! defined key_ice_lim && ! defined key_lim2 ) 
    7273      ! opa model ice freeze()       
    7374      DO jj = 1, jpj 
    … …  
    8788 
    8889         zsst(:,:) = gsst(:,:) / REAL( nfbulk ) * tmask(:,:,1) 
    89          CALL flx_blk( zsst )     
     90         IF ( .NOT. lk_lim2 )  CALL flx_blk( zsst )     
     91 
     92#if defined key_lim2 
     93         CALL flx_blk_2( zsst )     
     94#endif 
    9095      
    9196         gsst(:,:) = 0.     
    9297 
    93 # if ! defined key_ice_lim 
     98# if ( ! defined key_ice_lim && ! defined key_lim2 ) 
    9499         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
    95100            CALL prt_ctl_info(' Forcings ') 

  • branches/dev_002_LIM/NEMO/OPA_SRC/SBC/flx_oce.F90

    r719 r824  
    3030      p_emp             !: evaporation minus precipitation             
    3131 
    32 #elif defined key_ice_lim || defined key_flx_bulk_monthly || defined key_flx_bulk_daily || defined key_flx_core 
     32#elif defined key_ice_lim || defined key_lim2 || defined key_flx_bulk_monthly || defined key_flx_bulk_daily || defined key_flx_core 
    3333   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)    ::   &  !: 
    3434      qsr_ice  ,      &  !: solar flux over ice 

  • branches/dev_002_LIM/NEMO/OPA_SRC/SBC/flxblk.F90

    r789 r824  
    44   !! Ocean forcing:  bulk thermohaline forcing of the ocean (or ice) 
    55   !!===================================================================== 
    6 #if defined key_flx_bulk_monthly || defined key_flx_bulk_daily 
     6#if  defined key_flx_bulk_monthly || defined key_flx_bulk_daily  
    77   !!---------------------------------------------------------------------- 
    88   !!   'key_flx_bulk_monthly'   or                            MONTHLY bulk 

  • branches/dev_002_LIM/NEMO/OPA_SRC/SBC/flxfwb.F90

    r719 r824  
    2323   USE blk_oce 
    2424   USE flxblk          ! bulk formulea 
     25   USE flxblk_2          ! bulk formulea 
    2526   USE daymod          ! calendar 
    2627 

  • branches/dev_002_LIM/NEMO/OPA_SRC/SBC/ocesbc.F90

    r719 r824  
    6464CONTAINS 
    6565 
    66 #if defined key_ice_lim 
    67    !!---------------------------------------------------------------------- 
    68    !!   'key_ice_lim' :                                   LIM sea-ice model 
     66#if defined key_ice_lim || defined key_lim2 
     67   !!---------------------------------------------------------------------- 
     68   !!   'key_ice_lim'     :                              LIM2 sea-ice model 
     69   !!   'key_lim2'        :                              LIM3 sea-ice model 
    6970   !!---------------------------------------------------------------------- 
    7071# if defined key_coupled 
    … …  
    816817      !!---------------------------------------------------------------------- 
    817818 
    818 #if defined key_ice_lim 
     819#if defined key_ice_lim || defined key_lim2 
    819820      ! sea ice indicator (1 or 0) 
    820821      DO jj = 1, jpj 

  • branches/dev_002_LIM/NEMO/OPA_SRC/SBC/taumod.F90

    r719 r824  
    5959   !!---------------------------------------------------------------------- 
    6060   !!   'key_oasis3' or 'key_oasis4' and           Coupled Ocean/Atmosphere 
    61    !!   'key_ice_lim'                                   LIM sea-ice 
     61   !!   'key_ice_lim'or 'key_lim2'                              LIM sea-ice 
    6262   !!---------------------------------------------------------------------- 
    6363   ! New way: 3D referential link to the earth (avoid north pole pb) 

  • branches/dev_002_LIM/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r789 r824  
    147147               zind(ji,jj,jk) =  MAX ( upsrnfh(ji,jj) * upsrnfz(jk),     &  ! changing advection scheme near runoff 
    148148                  &                    upsadv(ji,jj)                     &  ! in the vicinity of some straits 
    149 #if defined key_ice_lim 
     149#if defined key_ice_lim || defined key_lim2 
    150150                  &                  , tmask(ji,jj,jk)                   &  ! half upstream tracer fluxes 
    151151                  &                  * MAX( 0., SIGN( 1., fzptn(ji,jj)   &  ! if tn < ("freezing"+0.1 ) 

  • branches/dev_002_LIM/NEMO/OPA_SRC/cpl.F90

    r719 r824  
    569569      sieoc(:,:) = sieoc(:,:) + freeze(:,:) 
    570570 
    571 #if defined key_ice_lim 
     571#if defined key_ice_lim || defined key_lim2 
    572572      alboc(:,:) = alboc(:,:) + freeze(:,:) * alb_ice(:,:) 
    573573      ticoc(:,:) = ticoc(:,:) + freeze(:,:) * tn_ice(:,:)  

  • branches/dev_002_LIM/NEMO/OPA_SRC/cpl_oce.F90

    r719 r824  
    218218      qsrc           !: solar radiation (w m-2) 
    219219 
    220 #  if defined key_ice_lim 
     220#  if defined key_ice_lim || defined key_lim3 
    221221   REAL(wp), DIMENSION(jpi,jpj) ::   &  !: 
    222222      watm        ,    &  !: 

  • branches/dev_002_LIM/NEMO/OPA_SRC/ice_oce.F90

    r719 r824  
    1111   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    1212   !!---------------------------------------------------------------------- 
    13 #if defined key_ice_lim 
     13#if defined key_lim2 || defined key_ice_lim 
    1414   !!---------------------------------------------------------------------- 
    15    !!   'key_ice_lim'   :                                     LIM ice model 
     15   !!   'key_lim2 or key_ice_lim '   :             LIM 2.0 or 3.0 ice model 
    1616   !!---------------------------------------------------------------------- 
    1717   !! * Modules used 
    … …  
    2323  
    2424   !! Shared module variables 
    25    LOGICAL, PUBLIC, PARAMETER ::   lk_ice_lim = .TRUE.    !: LIM ice model 
     25# if defined  key_lim2 
     26   LOGICAL, PUBLIC, PARAMETER ::   lk_lim2        = .TRUE.    !: LIM2 ice model 
     27   LOGICAL, PUBLIC, PARAMETER ::   lk_ice_lim     = .FALSE.   !: LIM3 ice model 
     28# else 
     29   LOGICAL, PUBLIC, PARAMETER ::   lk_lim2        = .FALSE.   !: LIM2 ice model 
     30   LOGICAL, PUBLIC, PARAMETER ::   lk_ice_lim     = .TRUE.    !: LIM3 ice model 
     31# endif 
    2632 
    2733   !!---------------------------------------------------------------------- 
    … …  
    5763   !!   Default option                                 NO LIM sea-ice model 
    5864   !!---------------------------------------------------------------------- 
    59    LOGICAL, PUBLIC, PARAMETER ::   lk_ice_lim = .FALSE.  !: No LIM ice model 
     65   LOGICAL, PUBLIC, PARAMETER ::   lk_lim2        = .FALSE.  !: No LIM 2.0 ice model 
     66   LOGICAL, PUBLIC, PARAMETER ::   lk_ice_lim     = .FALSE.  !: No LIM 3.0 ice model 
    6067#endif 
    6168 

  • branches/dev_002_LIM/NEMO/OPA_SRC/opa.F90

    r719 r824  
    5959 
    6060   USE phycst          ! physical constant                  (par_cst routine) 
     61#if defined key_ice_lim 
    6162   USE iceini          ! initialization of sea-ice         (ice_init routine) 
     63#endif 
     64#if defined key_lim2 
     65   USE iceini_2        ! initialization of sea-ice         (ice_init_2 routine) 
     66#endif 
    6267   USE cpl             ! coupled ocean/atmos.              (cpl_init routine) 
    6368   USE ocfzpt          ! ocean freezing point              (oc_fz_pt routine) 
    … …  
    283288 
    284289#if defined key_ice_lim 
    285       CALL ice_init                         ! Sea ice model 
     290      CALL ice_init                         ! Sea ice model LIM3 
     291#endif 
     292 
     293#if defined key_lim2 
     294      CALL ice_init_2                       ! Sea ice model LIM2 
    286295#endif 
    287296 

  • branches/dev_002_LIM/NEMO/OPA_SRC/restart.F90

    r783 r824  
    144144      CALL iom_rstput( kt, nitrst, numrow, 'hdivn'  , hdivn   ) 
    145145 
    146 #if defined key_ice_lim         
     146#if defined key_ice_lim || defined key_lim2 
    147147      CALL iom_rstput( kt, nitrst, numrow, 'nfice'  , REAL( nfice, wp) )   !  ice computation frequency 
    148148      CALL iom_rstput( kt, nitrst, numrow, 'sst_io' , sst_io  ) 
    … …  
    210210      !!---------------------------------------------------------------------- 
    211211      REAL(wp) ::   zcoef, zkt, zrdt, zrdttra1, zndastp, znfice, znfbulk 
    212 #if defined key_ice_lim 
     212#if defined key_ice_lim || defined key_lim2 
    213213      INTEGER  ::   ji, jj 
    214214#endif 
    … …  
    306306      !!sm: TO BE MOVED IN NEW SURFACE MODULE... 
    307307 
    308 #if defined key_ice_lim 
     308#if defined key_ice_lim || defined key_lim2 
    309309      ! Louvain La Neuve Sea Ice Model 
    310310      IF( iom_varid( numror, 'nfice', ldstop = .FALSE. ) > 0 ) then  

  • branches/dev_002_LIM/NEMO/OPA_SRC/step.F90

    r789 r824  
    9696   USE zpshde          ! partial step: hor. derivative     (zps_hde routine) 
    9797   USE ice_oce         ! sea-ice variable 
     98#if defined key_ice_lim 
    9899   USE icestp          ! sea-ice time-stepping             (ice_stp routine) 
    99  
     100#endif 
     101#if defined key_lim2 
     102   USE icestp_2        ! sea-ice time-stepping             (ice_stp_2 routine) 
     103#endif 
    100104   USE diawri          ! Standard run outputs             (dia_wri routine) 
    101105   USE trdicp          ! Ocean momentum/tracers trends    (trd_wri routine) 
    … …  
    206210                         CALL flx( kstp )             ! heat and freshwater fluxes 
    207211 
    208       IF( lk_ice_lim )   CALL ice_stp( kstp )         ! sea-ice model (Update stress & fluxes) 
     212#if defined key_ice_lim 
     213      CALL ice_stp( kstp )           ! sea-ice model (Update stress & fluxes) 
     214#endif 
     215#if defined key_lim2 
     216      CALL ice_stp_2( kstp )         ! sea-ice model (Update stress & fluxes) 
     217#endif 
    209218 
    210219                         CALL oce_sbc( kstp )         ! ocean surface boudaries 
Note: See TracChangeset for help on using the changeset viewer.