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

Changeset 7107


Ignore:
Timestamp:
2016-10-26T15:21:30+02:00 (7 years ago)
Author:
flavoni
Message:

#1692 and ROBUST-3, add ISF computation for TOOLS, and remove .nc files

Location:
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/TOOLS/DOMAINcfg
Files:
5 deleted
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90

    r6951 r7107  
    8787   !!---------------------------------------------------------------------- 
    8888   !! NEMO/OPA 3.3.1 , NEMO Consortium (2011) 
    89    !! $Id: domzgr.F90 6492 2016-04-22 13:04:31Z mathiot $ 
     89   !! $Id: domzgr.F90 6827 2016-08-01 13:37:15Z flavoni $ 
    9090   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    9191   !!---------------------------------------------------------------------- 
     
    488488         risfdep(:,:)=0.e0 
    489489         misfdep(:,:)=1 
     490         ! 
     491         ! (ISF) TODO build ice draft netcdf file for isomip and build the corresponding part of code 
     492         IF( cp_cfg == "isomip" .AND. ln_isfcav ) THEN  
     493            risfdep(:,:)=200.e0  
     494            misfdep(:,:)=1  
     495            ij0 = 1 ; ij1 = 40  
     496            DO jj = mj0(ij0), mj1(ij1)  
     497               risfdep(:,jj)=700.0_wp-(gphit(:,jj)+80.0_wp)*125.0_wp  
     498            END DO  
     499            WHERE( bathy(:,:) <= 0._wp )  risfdep(:,:) = 0._wp  
     500         !  
     501         ELSEIF ( cp_cfg == "isomip2" .AND. ln_isfcav ) THEN 
     502         !  
     503            risfdep(:,:)=0.e0 
     504            misfdep(:,:)=1 
     505            ij0 = 1 ; ij1 = 40 
     506            DO jj = mj0(ij0), mj1(ij1) 
     507               risfdep(:,jj)=700.0_wp-(gphit(:,jj)+80.0_wp)*125.0_wp 
     508            END DO 
     509            WHERE( bathy(:,:) <= 0._wp )  risfdep(:,:) = 0._wp 
     510         END IF 
    490511         ! 
    491512         DEALLOCATE( idta, zdta ) 
Note: See TracChangeset for help on using the changeset viewer.