Changeset 5920 for branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2015-11-25T17:58:51+01:00 (8 years ago)

Author:

mathiot

Message:

ice sheet coupling: add treshold to define grounded area, remove useless va
riable, change some variable name + add some namelist parameter

File:

• branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (2 diffs)

branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -531 +531 @@
                ! set grounded point to 0 
                ! (a treshold could be set here if needed, or set it offline based on the grounded fraction)
-               WHERE ( bathy(:,:) .LE. risfdep(:,:) )
+               WHERE ( bathy(:,:) .LE. risfdep(:,:) + rn_isfhmin )
                   misfdep(:,:) = 0 ; risfdep(:,:) = 0._wp
                   mbathy (:,:) = 0 ; bathy  (:,:) = 0._wp
@@ -1258 +1258 @@
       DO jl = 1, 10     
          ! check at each iteration if isf is grounded or not (1cm treshold have to be update after first coupling experiments)
-         WHERE (bathy(:,:) .LE. risfdep(:,:)+1e-2 )
+         WHERE (bathy(:,:) .LE. risfdep(:,:) + rn_isfhmin)
             misfdep(:,:) = 0 ; risfdep(:,:) = 0._wp
             mbathy (:,:) = 0 ; bathy  (:,:) = 0._wp