Changeset 391


Ignore:
Timestamp:
2006-03-10T11:35:43+01:00 (15 years ago)
Author:
opalod
Message:

RB:nemo_v1_update_038: first integration of AGRIF :

add agrif to ice files

Location:
trunk/NEMO/LIM_SRC
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC/iceini.F90

    r253 r391  
    6060      !!   8.5  !  02-08  (G. Madec)  F90: Free form and modules 
    6161      !!---------------------------------------------------------------------- 
    62  
     62       CHARACTER(len=80) :: namelist_icename 
     63        
    6364      ! Open the namelist file  
    64       OPEN( numnam_ice, FILE= 'namelist_ice', FORM='formatted', STATUS = 'old')  
     65      namelist_icename = 'namelist_ice' 
     66            
     67      CALL ctlopn(numnam_ice,namelist_icename,'OLD', 'FORMATTED', 'SEQUENTIAL',   & 
     68                     1,numout,.FALSE.,1)       
    6569 
    6670      CALL ice_run                    !  read in namelist some run parameters 
     
    114118      !!   2.0  !  03-08 (C. Ethe)  Original code 
    115119      !!------------------------------------------------------------------- 
     120 
    116121      NAMELIST/namicerun/ ln_limdyn, acrit, hsndif, hicdif 
    117122      !!------------------------------------------------------------------- 
     
    132137         WRITE(numout,*) '   computation of temp. in ice  (=0) or not (=9999) hicdif = ', hicdif 
    133138      ENDIF 
    134        
    135139   END SUBROUTINE ice_run 
    136140 
  • trunk/NEMO/LIM_SRC/limrhg.F90

    r258 r391  
    4646      !!  atmospheric (wind stress) and oceanic (water stress and surface 
    4747      !!  tilt) forcings. Ice-ice interaction is described by a non-linear 
    48       !!  viscous-plastic law including shear strength and a bulk rheology.    
     48      !!  viscous-plastic law including shear strength and a bulk rheology. 
    4949      !! 
    5050      !! ** Action  : - compute u_ice, v_ice the sea-ice velocity 
     
    501501            END DO 
    502502 
    503             ! Terms involving already up-dated velocities.      
     503            ! Terms involving already up-dated velocities. 
    504504            !-Using the arrays zu_ice and zv_ice in the computation of the terms ze leads to JACOBI's method;  
    505             ! Using arrays u and v in the computation of the terms ze leads to GAUSS-SEIDEL method.    
     505            ! Using arrays u and v in the computation of the terms ze leads to GAUSS-SEIDEL method. 
    506506              
    507507            DO jj = k_j1+1, k_jpj-1 
     
    530530                  zs21(ji,jj,1,1) =  zvis12 * ze21 + zvis21 * ze12 
    531531 
     532#if defined key_agrif 
     533             END DO 
     534          END DO 
     535 
     536          DO jj = k_j1+1, k_jpj-1 
     537             DO ji = 2, jpim1 
     538#endif 
    532539 
    533540                  iim1 = ji 
     
    571578                  zs21(ji,jj,1,2) =  zs21(ji,jj,1,2) + zvis12 * ze21 + zvis21 * ze12 
    572579 
    573 !i             END DO 
    574 !i          END DO 
    575  
    576 !i          DO jj = k_j1+1, k_jpj-1 
    577 !i             DO ji = 2, jpim1 
     580#if defined key_agrif 
     581             END DO 
     582          END DO 
     583 
     584          DO jj = k_j1+1, k_jpj-1 
     585             DO ji = 2, jpim1 
     586#endif 
    578587                  zd1(ji,jj) =   & 
    579588                     + alambd(ji,jj,2,2,2,1) * zs11(ji,jj,2,1) + alambd(ji,jj,2,2,2,2) * zs11(ji,jj,2,2)  & 
  • trunk/NEMO/LIM_SRC/limrst.F90

    r352 r391  
    7474         zsec, zdate0, zdt 
    7575 
    76       CHARACTER(len=45)  ::  & 
    77          ccfile = 'restart_ice_out.nc' 
     76      CHARACTER(len=45)  ::  ccfile 
     77 
     78      ccfile = 'restart_ice_out.nc' 
     79 
     80#if defined key_agrif 
     81      if ( .NOT. Agrif_Root() ) then 
     82         ccfile= TRIM(Agrif_CFixed())//'_'//TRIM(ccfile) 
     83      endif 
     84#endif 
    7885 
    7986      inumwrs  = 61 
     
    172179 
    173180      !- dummy variables : 
    174       CHARACTER(len=45)  ::  & 
    175          ccfile = 'restart_ice_in.nc' 
     181      CHARACTER(len=45)  ::  ccfile 
    176182      INTEGER :: & 
    177183        ji, jj 
     
    193199         clvnames(60)        
    194200 
     201 
     202       ccfile = 'restart_ice_in.nc' 
     203#if defined key_agrif 
     204      if ( .NOT. Agrif_Root() ) then 
     205         ccfile= TRIM(Agrif_CFixed())//'_'//TRIM(ccfile) 
     206      endif 
     207#endif 
    195208 
    196209      !Initialisations 
  • trunk/NEMO/LIM_SRC/limrst_dimg.h90

    r247 r391  
    3939      !! no additional check is necessary. 
    4040 
     41        
    4142       ! 0. Initializations 
    4243       ! ------------------ 
     
    5253       irecl8= jpi * jpj * wp 
    5354       WRITE(clres,'(a,".",i3.3)') TRIM(ccfile),narea 
    54        OPEN(inum,FILE=clres,FORM='UNFORMATTED', ACCESS='DIRECT', RECL=irecl8 ) 
     55       CALL ctlopn(inum,clres,'UNKNOWN','UNFORMATTED','DIRECT',irecl8,numout,lwp,0) 
     56!!       OPEN(inum,FILE=clres,FORM='UNFORMATTED', ACCESS='DIRECT', RECL=irecl8 ) 
    5557 
    5658 
     
    165167       ! Open direct access file, with reclength for 2D wp fields 
    166168       WRITE(clres,'(a,".",i3.3)') TRIM(ccfile),narea 
    167        OPEN(inum,FILE=clres,FORM='UNFORMATTED', ACCESS='DIRECT', RECL=8) 
     169       CALL ctlopn(inum,clres,'UNKNOWN','UNFORMATTED','DIRECT',8,numout,lwp,1) 
     170   !!!    OPEN(inum,FILE=clres,FORM='UNFORMATTED', ACCESS='DIRECT', RECL=8) 
    168171       READ(inum,REC=1) irecl8 
    169172       CLOSE(inum) 
Note: See TracChangeset for help on using the changeset viewer.