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 11893 for NEMO/branches/2019/dev_r11756_SI3restart_XIOS/src/ICE/icerst.F90 – NEMO

Ignore:
Timestamp:
2019-11-13T09:48:43+01:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2323 read/write SI3 restart file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11756_SI3restart_XIOS/src/ICE/icerst.F90

    r11870 r11893  
    138138      CALL iom_rstput( iter, nitrst, numriw, 'a_ip' , a_ip, ldxios = lwxios  ) 
    139139      CALL iom_rstput( iter, nitrst, numriw, 'v_ip' , v_ip, ldxios = lwxios  ) 
     140      IF(lwxios) THEN 
     141         CALL iom_rstput( iter, nitrst, numriw, 'e_s' , e_s, ldxios = lwxios ) 
     142         CALL iom_rstput( iter, nitrst, numriw, 'e_i' , e_i, ldxios = lwxios ) 
     143      ELSE 
    140144      ! Snow enthalpy 
    141       DO jk = 1, nlay_s  
    142          WRITE(zchar1,'(I2.2)') jk 
    143          znam = 'e_s'//'_l'//zchar1 
    144          z3d(:,:,:) = e_s(:,:,jk,:) 
    145          CALL iom_rstput( iter, nitrst, numriw, znam , z3d, ldxios = lwxios ) 
    146       END DO 
    147       ! Ice enthalpy 
    148       DO jk = 1, nlay_i  
    149          WRITE(zchar1,'(I2.2)') jk 
    150          znam = 'e_i'//'_l'//zchar1 
    151          z3d(:,:,:) = e_i(:,:,jk,:) 
    152          CALL iom_rstput( iter, nitrst, numriw, znam , z3d, ldxios = lwxios ) 
    153       END DO 
     145         DO jk = 1, nlay_s  
     146            WRITE(zchar1,'(I2.2)') jk 
     147            znam = 'e_s'//'_l'//zchar1 
     148            z3d(:,:,:) = e_s(:,:,jk,:) 
     149            CALL iom_rstput( iter, nitrst, numriw, znam , z3d, ldxios = lwxios ) 
     150         END DO 
     151         ! Ice enthalpy 
     152         DO jk = 1, nlay_i  
     153            WRITE(zchar1,'(I2.2)') jk 
     154            znam = 'e_i'//'_l'//zchar1 
     155            z3d(:,:,:) = e_i(:,:,jk,:) 
     156            CALL iom_rstput( iter, nitrst, numriw, znam , z3d, ldxios = lwxios ) 
     157         END DO 
     158      ENDIF 
    154159      ! fields needed for Met Office (Jules) coupling 
    155160      IF( ln_cpl ) THEN 
Note: See TracChangeset for help on using the changeset viewer.