Ignore:
Timestamp:
2017-09-13T18:46:56+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part6 - commits of the day

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerhg_evp.F90

    r8517 r8518  
    4141   PRIVATE 
    4242 
    43    PUBLIC   ice_rhg_evp        ! routine called by icerhg.F90 
     43   PUBLIC   ice_rhg_evp   ! called by icerhg.F90 
     44   PUBLIC   rhg_evp_rst   ! called by icerhg.F90 
    4445 
    4546   !! * Substitutions 
     
    108109      !!------------------------------------------------------------------- 
    109110      INTEGER, INTENT(in) ::   kt      ! time step 
    110       REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) ::   pstress1_i, pstress2_i, pstress12_i 
    111       REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) ::   pu_ice, pv_ice, pshear_i, pdivu_i, pdelta_i  
     111      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   pstress1_i, pstress2_i, pstress12_i 
     112      REAL(wp), DIMENSION(:,:), INTENT(  out) ::   pu_ice, pv_ice, pshear_i, pdivu_i, pdelta_i  
    112113      !! 
    113114      INTEGER ::   ji, jj       ! dummy loop indices 
     
    825826   END SUBROUTINE ice_rhg_evp 
    826827 
     828   SUBROUTINE rhg_evp_rst( cdrw, kt ) 
     829      !!--------------------------------------------------------------------- 
     830      !!                   ***  ROUTINE rhg_evp_rst  *** 
     831      !!                      
     832      !! ** Purpose :   Read or write RHG file in restart file 
     833      !! 
     834      !! ** Method  :   use of IOM library 
     835      !!---------------------------------------------------------------------- 
     836      CHARACTER(len=*) , INTENT(in) ::   cdrw   ! "READ"/"WRITE" flag 
     837      INTEGER, OPTIONAL, INTENT(in) ::   kt     ! ice time-step 
     838      ! 
     839      INTEGER  ::   iter            ! local integer 
     840      INTEGER  ::   id1, id2, id3   ! local integers 
     841      !!---------------------------------------------------------------------- 
     842      ! 
     843      IF( TRIM(cdrw) == 'READ' ) THEN        ! Read/initialize 
     844         !                                   ! --------------- 
     845         IF( ln_rstart ) THEN                   !* Read the restart file 
     846            ! 
     847            id1 = iom_varid( numrir, 'stress1_i' , ldstop = .FALSE. ) 
     848            id2 = iom_varid( numrir, 'stress2_i' , ldstop = .FALSE. ) 
     849            id3 = iom_varid( numrir, 'stress12_i', ldstop = .FALSE. ) 
     850            ! 
     851            IF( MIN( id1, id2, id3 ) > 0 ) THEN      ! fields exist 
     852               CALL iom_get( numrir, jpdom_autoglo, 'stress1_i' , stress1_i  ) 
     853               CALL iom_get( numrir, jpdom_autoglo, 'stress2_i' , stress2_i  ) 
     854               CALL iom_get( numrir, jpdom_autoglo, 'stress12_i', stress12_i ) 
     855            ELSE                                     ! start rheology from rest 
     856               IF(lwp) WRITE(numout,*) '   ==>>   previous run without rheology, set stresses to 0' 
     857               stress1_i (:,:) = 0._wp 
     858               stress2_i (:,:) = 0._wp 
     859               stress12_i(:,:) = 0._wp 
     860            ENDIF 
     861         ELSE                                   !* Start from rest 
     862            IF(lwp) WRITE(numout,*) '   ==>>   start from rest: set stresses to 0' 
     863            stress1_i (:,:) = 0._wp 
     864            stress2_i (:,:) = 0._wp 
     865            stress12_i(:,:) = 0._wp 
     866         ENDIF 
     867         ! 
     868      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
     869         !                                   ! ------------------- 
     870         IF(lwp) WRITE(numout,*) '---- rhg-rst ----' 
     871         iter = kt + nn_fsbc - 1             ! ice restarts are written at kt == nitrst - nn_fsbc + 1 
     872         ! 
     873         CALL iom_rstput( iter, nitrst, numriw, 'stress1_i' , stress1_i  ) 
     874         CALL iom_rstput( iter, nitrst, numriw, 'stress2_i' , stress2_i  ) 
     875         CALL iom_rstput( iter, nitrst, numriw, 'stress12_i', stress12_i ) 
     876         ! 
     877      ENDIF 
     878      ! 
     879   END SUBROUTINE rhg_evp_rst 
     880 
    827881#else 
    828882   !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.