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 9123 – NEMO

Changeset 9123


Ignore:
Timestamp:
2017-12-18T23:34:28+01:00 (6 years ago)
Author:
clem
Message:

solve issue with Prather scheme and restart

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icedyn_adv_pra.F90

    r9019 r9123  
    708708            CALL iom_get( numrir, jpdom_autoglo, 'sxyage', sxyage ) 
    709709            !                                                        ! open water in sea ice 
    710             CALL iom_get( numrir, jpdom_autoglo, 'sxopw ', sxopw  ) 
    711             CALL iom_get( numrir, jpdom_autoglo, 'syopw ', syopw  ) 
     710            CALL iom_get( numrir, jpdom_autoglo, 'sxopw' , sxopw  ) 
     711            CALL iom_get( numrir, jpdom_autoglo, 'syopw' , syopw  ) 
    712712            CALL iom_get( numrir, jpdom_autoglo, 'sxxopw', sxxopw ) 
    713713            CALL iom_get( numrir, jpdom_autoglo, 'syyopw', syyopw ) 
     
    716716            DO jk = 1, nlay_i  
    717717               WRITE(zchar1,'(I2.2)') jk 
    718                znam = 'sxe'//'_il'//zchar1   ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sxe (:,:,jk,:) = z3d(:,:,:) 
    719                znam = 'sye'//'_il'//zchar1   ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sye (:,:,jk,:) = z3d(:,:,:) 
    720                znam = 'sxxe'//'_il'//zchar1  ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sxxe(:,:,jk,:) = z3d(:,:,:) 
    721                znam = 'syye'//'_il'//zchar1  ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   syye(:,:,jk,:) = z3d(:,:,:) 
    722                znam = 'sxye'//'_il'//zchar1  ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sxye(:,:,jk,:) = z3d(:,:,:) 
     718               znam = 'sxe'//'_l'//zchar1   ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sxe (:,:,jk,:) = z3d(:,:,:) 
     719               znam = 'sye'//'_l'//zchar1   ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sye (:,:,jk,:) = z3d(:,:,:) 
     720               znam = 'sxxe'//'_l'//zchar1  ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sxxe(:,:,jk,:) = z3d(:,:,:) 
     721               znam = 'syye'//'_l'//zchar1  ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   syye(:,:,jk,:) = z3d(:,:,:) 
     722               znam = 'sxye'//'_l'//zchar1  ;   CALL iom_get( numrir, jpdom_autoglo, znam , z3d )   ;   sxye(:,:,jk,:) = z3d(:,:,:) 
    723723            END DO 
    724724            ! 
     
    802802         CALL iom_rstput( iter, nitrst, numriw, 'sxyage', sxyage ) 
    803803         !                                                           ! open water in sea ice 
    804          CALL iom_rstput( iter, nitrst, numriw, 'sxopw ', sxopw  ) 
    805          CALL iom_rstput( iter, nitrst, numriw, 'syopw ', syopw  ) 
     804         CALL iom_rstput( iter, nitrst, numriw, 'sxopw' , sxopw  ) 
     805         CALL iom_rstput( iter, nitrst, numriw, 'syopw' , syopw  ) 
    806806         CALL iom_rstput( iter, nitrst, numriw, 'sxxopw', sxxopw ) 
    807807         CALL iom_rstput( iter, nitrst, numriw, 'syyopw', syyopw ) 
     
    810810         DO jk = 1, nlay_i  
    811811            WRITE(zchar1,'(I2.2)') jk 
    812             znam = 'sxe'//'_il'//zchar1   ;   z3d(:,:,:) = sxe (:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
    813             znam = 'sye'//'_il'//zchar1   ;   z3d(:,:,:) = sye (:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
    814             znam = 'sxxe'//'_il'//zchar1  ;   z3d(:,:,:) = sxxe(:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
    815             znam = 'syye'//'_il'//zchar1  ;   z3d(:,:,:) = syye(:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
    816             znam = 'sxye'//'_il'//zchar1  ;   z3d(:,:,:) = sxye(:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
     812            znam = 'sxe'//'_l'//zchar1   ;   z3d(:,:,:) = sxe (:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
     813            znam = 'sye'//'_l'//zchar1   ;   z3d(:,:,:) = sye (:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
     814            znam = 'sxxe'//'_l'//zchar1  ;   z3d(:,:,:) = sxxe(:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
     815            znam = 'syye'//'_l'//zchar1  ;   z3d(:,:,:) = syye(:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
     816            znam = 'sxye'//'_l'//zchar1  ;   z3d(:,:,:) = sxye(:,:,jk,:)   ;   CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
    817817         END DO 
    818818         ! 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icerst.F90

    r9122 r9123  
    135135      ! Snow enthalpy (1st snow layer only) 
    136136      z3d = e_s(:,:,1,:) 
    137       CALL iom_rstput( iter, nitrst, numriw, 'tempt_sl1' , z3d ) 
     137      CALL iom_rstput( iter, nitrst, numriw, 'e_s_l01' , z3d ) 
    138138      ! Ice enthalpy (all ice layers) 
    139139      DO jk = 1, nlay_i  
    140140         WRITE(zchar1,'(I2.2)') jk 
    141          znam = 'tempt'//'_il'//zchar1 
     141         znam = 'e_i'//'_l'//zchar1 
    142142         z3d(:,:,:) = e_i(:,:,jk,:) 
    143143         CALL iom_rstput( iter, nitrst, numriw, znam , z3d ) 
     
    221221!!gm dangerous !!!!!  ===>>>> better reading writing all snow layers ! 
    222222      ! Snow enthalpy (1st snow layer only) 
    223       CALL iom_get( numrir, jpdom_autoglo, 'tempt_sl1' , z3d ) 
     223      CALL iom_get( numrir, jpdom_autoglo, 'e_s_l01' , z3d ) 
    224224      e_s(:,:,1,:) = z3d 
    225225      ! Ice enthalpy (all ice layers) 
    226226      DO jk = 1, nlay_i  
    227227         WRITE(zchar1,'(I2.2)') jk 
    228          znam = 'tempt'//'_il'//zchar1 
     228         znam = 'e_i'//'_l'//zchar1 
    229229         CALL iom_get( numrir, jpdom_autoglo, znam , z3d ) 
    230230         e_i(:,:,jk,:) = z3d(:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.