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 11397 for NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icevar.F90 – NEMO

Ignore:
Timestamp:
2019-08-02T18:31:52+02:00 (5 years ago)
Author:
clem
Message:

rewrite iceistate.F90, add new variables (optional) for input fields, and prepare the code for a simplified ice restart

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icevar.F90

    r11362 r11397  
    4747   !!   ice_var_zapneg    : remove negative ice fields 
    4848   !!   ice_var_roundoff  : remove negative values arising from roundoff erros 
    49    !!   ice_var_itd       : convert 1-cat to jpl-cat 
    50    !!   ice_var_itd2      : convert N-cat to jpl-cat 
     49   !!   ice_var_itd       : convert N-cat to M-cat 
    5150   !!   ice_var_bv        : brine volume 
    5251   !!   ice_var_enthalpy  : compute ice and snow enthalpies from temperature 
     
    935934               IF ( i_fill == 1 ) THEN      !-- case very thin ice: fill only category 1 
    936935                  ph_i(ji,1) = phti(ji) 
    937                   pa_i (ji,1) = pati (ji) 
     936                  pa_i(ji,1) = pati (ji) 
    938937               ELSE                         !-- case ice is thicker: fill categories >1 
    939938                  ! thickness 
     
    978977               IF ( zconv < epsi06 )   itest(2) = 1                                        ! Test 2: volume conservation 
    979978               ! 
    980                IF ( ph_i(ji,i_fill) >= hi_max(i_fill-1) )   itest(3) = 1                  ! Test 3: thickness of the last category is in-bounds ? 
     979               IF ( ph_i(ji,i_fill) >= hi_max(i_fill-1) )   itest(3) = 1                   ! Test 3: thickness of the last category is in-bounds ? 
    981980               ! 
    982981               itest(4) = 1 
    983982               DO jl = 1, i_fill 
    984                   IF ( pa_i(ji,jl) < 0._wp ) itest(4) = 0                                ! Test 4: positivity of ice concentrations 
     983                  IF ( pa_i(ji,jl) < 0._wp ) itest(4) = 0                                  ! Test 4: positivity of ice concentrations 
    985984               END DO 
    986985               !                                         !---------------------------- 
Note: See TracChangeset for help on using the changeset viewer.