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

Timestamp:

2015-10-22T17:18:18+02:00 (9 years ago)

Author:

mathiot

Message:

ice sheet coupling : add comments, rm PRINT statements, cosmetic changes

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

@@ -529 +529 @@
                WHERE( bathy(:,:) <= 0._wp )  risfdep(:,:) = 0._wp
 
-               ! set grounded point to 0 (treshold at 1cm, have to be update after first coupling experience)
-               WHERE (bathy(:,:) .LE. risfdep(:,:)+1e-2 )
+               ! 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(:,:) )
                   misfdep(:,:) = 0 ; risfdep(:,:) = 0._wp
                   mbathy (:,:) = 0 ; bathy  (:,:) = 0._wp
@@ -575 +576 @@
          ! patch to avoid case bathy = ice shelf draft and bathy between 0 and zhmin
          IF ( ln_isfcav ) THEN
-            WHERE (bathy == risfdep)
+            WHERE ( bathy == risfdep )
                bathy   = 0.0_wp ; risfdep = 0.0_wp
             END WHERE