Ignore:
Timestamp:
Mar 25, 2021, 8:07:25 AM (4 years ago)
Author:
emillour
Message:

Generic GCM:
Small fix on the reading/initialisation of Non orographic GW related fields
loaded by lect_start_archive
EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F

    r2336 r2484  
    966966         PRINT*, "Setting it to zero"
    967967         du_nonoro_gwdold(:,:,:)=0
    968       ENDIF
    969 #ifdef NC_DOUBLE
    970       ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,count,du_nonoro_gwdold)
    971 #else
    972       ierr = NF_GET_VARA_REAL(nid,nvarid,start,count,du_nonoro_gwdold)
    973 #endif
    974       IF (ierr .NE. NF_NOERR) THEN
    975          PRINT*, "lect_start_archive: Failed loading <du_nonoro_gwd>"
    976          CALL abort
     968      ELSE
     969#ifdef NC_DOUBLE
     970        ierr=NF_GET_VARA_DOUBLE(nid,nvarid,start,count,du_nonoro_gwdold)
     971#else
     972        ierr=NF_GET_VARA_REAL(nid,nvarid,start,count,du_nonoro_gwdold)
     973#endif
     974        IF (ierr .NE. NF_NOERR) THEN
     975          PRINT*, "lect_start_archive: Failed loading <du_nonoro_gwd>"
     976          CALL abort
     977        ENDIF
    977978      ENDIF
    978979
     
    983984         PRINT*, "Setting it to zero"
    984985         dv_nonoro_gwdold(:,:,:)=0
    985       ENDIF
    986 #ifdef NC_DOUBLE
    987       ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,count,dv_nonoro_gwdold)
    988 #else
    989       ierr = NF_GET_VARA_REAL(nid,nvarid,start,count,dv_nonoro_gwdold)
    990 #endif
    991       IF (ierr .NE. NF_NOERR) THEN
    992          PRINT*, "lect_start_archive: Failed loading <dv_nonoro_gwd>"
    993          CALL abort
     986      ELSE
     987#ifdef NC_DOUBLE
     988        ierr=NF_GET_VARA_DOUBLE(nid,nvarid,start,count,dv_nonoro_gwdold)
     989#else
     990        ierr=NF_GET_VARA_REAL(nid,nvarid,start,count,dv_nonoro_gwdold)
     991#endif
     992        IF (ierr .NE. NF_NOERR) THEN
     993          PRINT*, "lect_start_archive: Failed loading <dv_nonoro_gwd>"
     994          CALL abort
     995        ENDIF
    994996      ENDIF
    995997
     
    10001002         PRINT*, "Setting it to zero"
    10011003         east_gwstressold(:,:,:)=0
    1002       ENDIF
    1003 #ifdef NC_DOUBLE
    1004       ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,count,east_gwstressold)
    1005 #else
    1006       ierr = NF_GET_VARA_REAL(nid,nvarid,start,count,east_gwstressold)
    1007 #endif
    1008       IF (ierr .NE. NF_NOERR) THEN
    1009          PRINT*, "lect_start_archive: Failed loading <east_gwstress>"
    1010          CALL abort
     1004      ELSE
     1005#ifdef NC_DOUBLE
     1006        ierr=NF_GET_VARA_DOUBLE(nid,nvarid,start,count,east_gwstressold)
     1007#else
     1008        ierr=NF_GET_VARA_REAL(nid,nvarid,start,count,east_gwstressold)
     1009#endif
     1010        IF (ierr .NE. NF_NOERR) THEN
     1011          PRINT*, "lect_start_archive: Failed loading <east_gwstress>"
     1012          CALL abort
     1013        ENDIF
    10111014      ENDIF
    10121015
     
    10171020         PRINT*, "Setting it to zero"
    10181021         west_gwstressold(:,:,:)=0
    1019       ENDIF
    1020 #ifdef NC_DOUBLE
    1021       ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,count,west_gwstressold)
    1022 #else
    1023       ierr = NF_GET_VARA_REAL(nid,nvarid,start,count,west_gwstressold)
    1024 #endif
    1025       IF (ierr .NE. NF_NOERR) THEN
    1026          PRINT*, "lect_start_archive: Failed loading <west_gwstress>"
    1027          CALL abort
     1022      ELSE
     1023#ifdef NC_DOUBLE
     1024        ierr=NF_GET_VARA_DOUBLE(nid,nvarid,start,count,west_gwstressold)
     1025#else
     1026        ierr=NF_GET_VARA_REAL(nid,nvarid,start,count,west_gwstressold)
     1027#endif
     1028        IF (ierr .NE. NF_NOERR) THEN
     1029          PRINT*, "lect_start_archive: Failed loading <west_gwstress>"
     1030          CALL abort
     1031        ENDIF
    10281032      ENDIF
    10291033
Note: See TracChangeset for help on using the changeset viewer.