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_extrapolation.f90

    r2143 r10383  
    4343       CALL get_detected_pts(G0%gphiu,G1%gphiu,G0%glamu,G1%glamu,   & 
    4444            G0%umask(:,:,n),G1%umask(:,:,n),detected(:,:))        
    45     ELSE IF( PRESENT(posvar) .AND. posvar == 'V' ) THEN 
    46        !       
     45    ELSEIF( PRESENT(posvar) .AND. posvar == 'V' ) THEN 
    4746       CALL get_detected_pts(G0%gphiv,G1%gphiv,G0%glamv,G1%glamv,   & 
    4847            G0%vmask(:,:,n),G1%vmask(:,:,n),detected(:,:))                                   
     48    ELSEIF( PRESENT(posvar) .AND. posvar == 'F' ) THEN 
     49       CALL get_detected_pts(G0%gphif,G1%gphif,G0%glamf,G1%glamf,   & 
     50            G0%fmask(:,:,n),G1%fmask(:,:,n),detected(:,:))                                   
    4951    ELSE 
    50        CALL get_detected_pts(G0%nav_lat,G1%nav_lat,G0%nav_lon,G1%nav_lon,   & 
    51             G0%tmask(:,:,n),G1%tmask(:,:,n),detected(:,:))         
     52       CALL get_detected_pts(G0%gphit,G1%gphit,G0%glamt,G1%glamt,   & 
     53            G0%tmask(:,:,n),G1%tmask(:,:,n),detected(:,:))                                   
     54!clem       CALL get_detected_pts(G0%nav_lat,G1%nav_lat,G0%nav_lon,G1%nav_lon,   & 
     55!clem            G0%tmask(:,:,n),G1%tmask(:,:,n),detected(:,:))         
    5256    ENDIF 
    5357    !       
     
    112116       lat = G0%gphiv  
    113117       lon = G0%glamv 
     118    ELSE IF ( PRESENT(posvar) .AND. posvar == 'F' ) THEN 
     119       WHERE(G0%fmask(:,:,k) == 1. ) 
     120          newmask(:,:) = .TRUE. 
     121       ELSEWHERE 
     122          newmask(:,:) = .FALSE. 
     123       END WHERE 
     124       ALLOCATE(mask(SIZE(G0%fmask,1),SIZE(G0%fmask,2),SIZE(G0%fmask,3))) 
     125       ALLOCATE(lat(SIZE(G0%fmask,1),SIZE(G0%fmask,2))) 
     126       ALLOCATE(lon(SIZE(G0%fmask,1),SIZE(G0%fmask,2))) 
     127       mask = G0%fmask         
     128       lat = G0%gphif  
     129       lon = G0%glamf 
    114130    ELSE       
    115131       WHERE(G0%tmask(:,:,k) == 1. ) 
     
    121137       ALLOCATE(lat(SIZE(G0%tmask,1),SIZE(G0%tmask,2))) 
    122138       ALLOCATE(lon(SIZE(G0%tmask,1),SIZE(G0%tmask,2))) 
    123        mask = G0%tmask 
    124        lon = G0%nav_lon 
    125        lat = G0%nav_lat 
     139       mask = G0%tmask         
     140       lat = G0%gphit  
     141       lon = G0%glamt 
     142!clem       mask = G0%tmask 
     143!clem       lon = G0%nav_lon 
     144!clem       lat = G0%nav_lat 
    126145    ENDIF 
    127146    ! 
Note: See TracChangeset for help on using the changeset viewer.