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 2874 for branches/2011/dev_r2802_UKMO8_cice/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90 – NEMO

Ignore:
Timestamp:
2011-09-28T12:19:59+02:00 (13 years ago)
Author:
charris
Message:

Code for running NEMO with CICE (for fully coupled mode this should be used in combination with dev_r2802_UKMO8_sbccpl). Changes are described briefly below.

physct: Constants modified to be consistent with CICE

nemogcm / prtctl / mppini: Changes to NEMO decomposition (activated using key_nemocice_decomp) to produce 'square' options in CICE. Can run without this key / code but this requires a global gather / scatter in the NEMO-CICE coupling which gets very slow on large processors numbers.

sbc_ice: CICE options and arrays added

sbcmod: CICE option added, including calls for initialising and finalising CICE.

sbcblk_core: Make sure necessary forcing field are available for CICE

sbcice_cice: Main CICE coupling code.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_r2802_UKMO8_cice/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

    r2528 r2874  
    4848#endif 
    4949 
     50#if defined key_cice 
     51   REAL(wp), PUBLIC ::   rau0     = 1026._wp      !: reference volumic mass (density)  (kg/m3) 
     52#else 
    5053   REAL(wp), PUBLIC ::   rau0     = 1035._wp      !: reference volumic mass (density)  (kg/m3) 
     54#endif 
    5155   REAL(wp), PUBLIC ::   rau0r                    !: reference specific volume         (m3/kg) 
    5256   REAL(wp), PUBLIC ::   rcp      =    4.e+3_wp   !: ocean specific heat 
    5357   REAL(wp), PUBLIC ::   ro0cpr                   !: = 1. / ( rau0 * rcp ) 
    5458 
    55 #if defined key_lim3 
     59#if defined key_lim3 || defined key_cice 
    5660   REAL(wp), PUBLIC ::   rcdsn   =   0.31_wp      !: thermal conductivity of snow 
    5761   REAL(wp), PUBLIC ::   rcdic   =   2.034396_wp  !: thermal conductivity of fresh ice 
     
    100104      rsiyea = 365.25 * rday * 2. * rpi / 6.283076 
    101105      rsiday = rday / ( 1. + rday / rsiyea ) 
     106#if defined key_cice 
     107      omega =  7.292116e-05 
     108#else 
    102109      omega  = 2. * rpi / rsiday  
     110#endif 
    103111 
    104112      rau0r  = 1. /   rau0   
Note: See TracChangeset for help on using the changeset viewer.