Ignore:
Timestamp:
2018-12-12T19:52:17+01:00 (23 months ago)
Author:
clem
Message:

ice restart should work in the nesting tools now. However ocean restart has been broken for some time

File:
1 edited

Legend:

Unmodified
Added
Removed
  • utils/tools/NESTING/src/agrif_readwrite.f90

    r10381 r10383  
    12231223    WRITE(*,*) 'Init masks in T points'     
    12241224    !      
    1225     ALLOCATE(Grid%tmask(nx,ny,N)) 
     1225    ALLOCATE( Grid%tmask(nx,ny,N) ) 
    12261226    ! 
    12271227    DO k = 1,N 
     
    12361236    !       
    12371237  END SUBROUTINE Init_Tmask 
     1238  ! 
     1239  !***************************************************** 
     1240  !   subroutine Init_ssmask(name,Grid) 
     1241  !***************************************************** 
     1242  ! 
     1243!  SUBROUTINE Init_ssmask(varname,filename,Grid,jpiglo,jpjglo) 
     1244!    ! 
     1245!    USE io_netcdf 
     1246!    !       
     1247!    CHARACTER(*) varname,filename 
     1248!    INTEGER :: nx,ny,k,i,j,jpiglo,jpjglo 
     1249!    TYPE(Coordinates) :: Grid 
     1250!    REAL*8, POINTER, DIMENSION(:,:) ::zwf => NULL() 
     1251!    ! 
     1252!    IF(jpiglo == 1 .AND. jpjglo == 1) THEN 
     1253!       CALL read_ncdf_var(varname,filename,Grid%bathy_level) 
     1254!    ELSE 
     1255!       CALL read_ncdf_var(varname,filename,Grid%bathy_level,(/jpizoom,jpjzoom/),(/jpiglo,jpjglo/) ) 
     1256!    ENDIF 
     1257!    ! 
     1258!    nx = SIZE(Grid%bathy_level,1) 
     1259!    ny = SIZE(Grid%bathy_level,2)  
     1260!    ! 
     1261!    WRITE(*,*) 'Init surface masks in T points'     
     1262!    !      
     1263!    ALLOCATE( Grid%ssmask(nx,ny), Grid%ssumask(nx,ny), Grid%ssvmask(nx,ny) ) 
     1264!    ! 
     1265!    WHERE(Grid%bathy_level(:,:) <= 0. )     
     1266!       Grid%ssmask(:,:) = 0. 
     1267!    ELSEWHERE 
     1268!       Grid%ssmask(:,:) = 1. 
     1269!    END WHERE 
     1270!    ! 
     1271!    Grid%ssumask(1:nx-1,:) = Grid%ssmask(1:nx-1,:)*Grid%ssmask(2:nx,:) 
     1272!    Grid%ssvmask(:,1:ny-1) = Grid%ssmask(:,1:ny-1)*Grid%ssmask(:,2:ny) 
     1273!    ! 
     1274!    Grid%ssumask(nx,:) = Grid%ssmask(nx,:) 
     1275!    Grid%ssvmask(:,ny) = Grid%ssmask(:,ny) 
     1276!    !       
     1277!  END SUBROUTINE Init_ssmask 
    12381278  ! 
    12391279  !***************************************************** 
Note: See TracChangeset for help on using the changeset viewer.