Ignore:
Timestamp:
Aug 21, 2015, 9:23:13 AM (9 years ago)
Author:
Ehouarn Millour
Message:

Physics/dynamics separation: get rid of all the 'include "temps.h"' in the physics; variables in module time_phylmdz_mod must be used instead. Also added JD_cur, JH_cur and JD_ref in module phys_cal_mod, in preparation for having physics handle its calendar internally.
EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/iophy.F90

    r2319 r2344  
    5454#endif
    5555  IMPLICIT NONE
    56   INCLUDE 'dimensions.h'
     56!  INCLUDE 'dimensions.h'
    5757    REAL,DIMENSION(klon),INTENT(IN) :: rlon
    5858    REAL,DIMENSION(klon),INTENT(IN) :: rlat
    … …  
    7777   
    7878!$OMP MASTER 
    79     ALLOCATE(io_lat(jjm+1-1/(iim*jjm)))
     79    ALLOCATE(io_lat(nbp_lat-1/(nbp_lon*(nbp_lat-1))))
    8080    io_lat(1)=rlat_glo(1)
    81     io_lat(jjm+1-1/(iim*jjm))=rlat_glo(klon_glo)
    82     IF ((iim*jjm) > 1) then
    83       DO i=2,jjm
    84         io_lat(i)=rlat_glo(2+(i-2)*iim)
     81    io_lat(nbp_lat-1/(nbp_lon*(nbp_lat-1)))=rlat_glo(klon_glo)
     82    IF ((nbp_lon*nbp_lat) > 1) then
     83      DO i=2,nbp_lat-1
     84        io_lat(i)=rlat_glo(2+(i-2)*nbp_lon)
    8585      ENDDO
    8686    ENDIF
    8787
    88     ALLOCATE(io_lon(iim))
    89     io_lon(:)=rlon_glo(2-1/(iim*jjm):iim+1-1/(iim*jjm))
     88    ALLOCATE(io_lon(nbp_lon))
     89    IF (klon_glo == 1) THEN
     90      io_lon(1)=rlon_glo(1)
     91    ELSE
     92      io_lon(1:nbp_lon)=rlon_glo(2:nbp_lon+1)
     93    ENDIF
     94
    9095!! (I) dtnb   : total number of domains
    9196!! (I) dnb    : domain number
    … …  
    104109
    105110    ddid=(/ 1,2 /)
    106     dsg=(/ iim, jjm+1-1/(iim*jjm) /)
    107     dsl=(/ iim, jj_nb /)
     111    dsg=(/ nbp_lon, nbp_lat /)
     112    dsl=(/ nbp_lon, jj_nb /)
    108113    dpf=(/ 1,jj_begin /)
    109     dpl=(/ iim, jj_end /)
     114    dpl=(/ nbp_lon, jj_end /)
    110115    dhs=(/ ii_begin-1,0 /)
    111116    IF (mpi_rank==mpi_size-1) THEN
    112117      dhe=(/0,0/)
    113118    ELSE
    114       dhe=(/ iim-ii_end,0 /) 
     119      dhe=(/ nbp_lon-ii_end,0 /) 
    115120    ENDIF
    116121
    … …  
    155160                                mpi_size, mpi_rank
    156161  USE ioipsl, only: flio_dom_set
     162  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    157163  IMPLICIT NONE
    158   INCLUDE 'dimensions.h'   
    159     REAL,DIMENSION(iim),INTENT(IN) :: lon
    160     REAL,DIMENSION(jjm+1-1/(iim*jjm)),INTENT(IN) :: lat
     164    REAL,DIMENSION(nbp_lon),INTENT(IN) :: lon
     165    REAL,DIMENSION(nbp_lat),INTENT(IN) :: lat
    161166
    162167    INTEGER,DIMENSION(2) :: ddid
    … …  
    169174
    170175!$OMP MASTER 
    171     allocate(io_lat(jjm+1-1/(iim*jjm)))
     176    allocate(io_lat(nbp_lat))
    172177    io_lat(:)=lat(:)
    173     allocate(io_lon(iim))
     178    allocate(io_lon(nbp_lon))
    174179    io_lon(:)=lon(:)
    175180   
    176181    ddid=(/ 1,2 /)
    177     dsg=(/ iim, jjm+1-1/(iim*jjm) /)
    178     dsl=(/ iim, jj_nb /)
     182    dsg=(/ nbp_lon, nbp_lat /)
     183    dsl=(/ nbp_lon, jj_nb /)
    179184    dpf=(/ 1,jj_begin /)
    180     dpl=(/ iim, jj_end /)
     185    dpl=(/ nbp_lon, jj_end /)
    181186    dhs=(/ ii_begin-1,0 /)
    182187    if (mpi_rank==mpi_size-1) then
    183188      dhe=(/0,0/)
    184189    else
    185       dhe=(/ iim-ii_end,0 /) 
     190      dhe=(/ nbp_lon-ii_end,0 /) 
    186191    endif
    187192   
    … …  
    198203  USE mod_phys_lmdz_para, only: is_sequential, is_using_mpi, is_mpi_root, &
    199204                                jj_begin, jj_end, jj_nb
     205  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    200206  use ioipsl, only: histbeg
    201207#ifdef CPP_XIOS
    … …  
    203209#endif
    204210  IMPLICIT NONE
    205   include 'dimensions.h'
    206211  include 'clesphys.h'
    207212   
    … …  
    217222!$OMP MASTER   
    218223    if (is_sequential) then
    219       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    220                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
     224      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     225                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
    221226    else
    222       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    223                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
     227      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     228                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
    224229    endif
    225230
    … …  
    240245
    241246  USE mod_phys_lmdz_para, only: jj_begin, jj_end, jj_nb, is_sequential
     247  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    242248  use ioipsl, only: histbeg
    243249
    244250  IMPLICIT NONE
    245   include 'dimensions.h'
    246251   
    247252    character*(*), INTENT(IN) :: name
    … …  
    255260#ifndef CPP_IOIPSL_NO_OUTPUT
    256261    if (is_sequential) then
    257       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    258                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
     262      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     263                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
    259264    else
    260       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    261                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
     265      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     266                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
    262267    endif
    263268#endif
    … …  
    274279                                is_sequential, klon_mpi_begin, klon_mpi_end, &
    275280                                mpi_rank
    276   USE mod_grid_phy_lmdz, only: klon_glo
     281  USE mod_grid_phy_lmdz, only: klon_glo, nbp_lon, nbp_lat
    277282  use ioipsl, only: histbeg
    278283
    279284  IMPLICIT NONE
    280   include 'dimensions.h'
    281285
    282286    REAL,DIMENSION(klon),INTENT(IN) :: rlon
    … …  
    304308    REAL, allocatable, DIMENSION(:) :: npplat, npplon
    305309    REAL, allocatable, DIMENSION(:,:) :: npplat_bounds, npplon_bounds
    306     INTEGER, PARAMETER :: jjmp1=jjm+1-1/jjm
    307     REAL, DIMENSION(iim,jjmp1) :: zx_lon, zx_lat
     310    REAL, DIMENSION(nbp_lon,nbp_lat) :: zx_lon, zx_lat
    308311
    309312    CALL gather(rlat,rlat_glo)
    … …  
    330333     endif
    331334!
    332      IF ( tabij(i).LE.iim) THEN
     335     IF ( tabij(i).LE.nbp_lon) THEN
    333336      plat_bounds(i,1)=rlat_glo(tabij(i))
    334337     ELSE
    335       plat_bounds(i,1)=rlat_glo(tabij(i)-iim)
     338      plat_bounds(i,1)=rlat_glo(tabij(i)-nbp_lon)
    336339     ENDIF
    337      plat_bounds(i,2)=rlat_glo(tabij(i)+iim)
     340     plat_bounds(i,2)=rlat_glo(tabij(i)+nbp_lon)
    338341!
    339342!    print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2)
    … …  
    351354     ENDDO
    352355
    353        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon_glo,zx_lon)
    354        if ((iim*jjm).gt.1) then
    355        DO i = 1, iim
     356       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon)
     357       if ((nbp_lon*nbp_lat).gt.1) then
     358       DO i = 1, nbp_lon
    356359         zx_lon(i,1) = rlon_glo(i+1)
    357          zx_lon(i,jjmp1) = rlon_glo(i+1)
     360         zx_lon(i,nbp_lat) = rlon_glo(i+1)
    358361       ENDDO
    359362       endif
    360        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat_glo,zx_lat)
     363       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat)
    361364
    362365    DO i=1,pim
    … …  
    367370
    368371     if (ipt(i).EQ.1) then
    369       plon_bounds(i,1)=zx_lon(iim,jpt(i))
     372      plon_bounds(i,1)=zx_lon(nbp_lon,jpt(i))
    370373      plon_bounds(i,2)=360.+zx_lon(ipt(i)+1,jpt(i))
    371374     endif
    372375 
    373      if (ipt(i).EQ.iim) then
     376     if (ipt(i).EQ.nbp_lon) then
    374377      plon_bounds(i,2)=360.+zx_lon(1,jpt(i))
    375378     endif
    … …  
    383386     endif
    384387 
    385      if (jpt(i).EQ.jjmp1) then
    386       plat_bounds(i,1)=zx_lat(ipt(i),jjmp1)+0.001
    387       plat_bounds(i,2)=zx_lat(ipt(i),jjmp1)-0.001
     388     if (jpt(i).EQ.nbp_lat) then
     389      plat_bounds(i,1)=zx_lat(ipt(i),nbp_lat)+0.001
     390      plat_bounds(i,2)=zx_lat(ipt(i),nbp_lat)-0.001
    388391     endif
    389392!
    … …  
    451454    use phys_output_var_mod, only: type_ecri, zoutm, zdtime_moy, lev_files, &
    452455                                   nid_files, nhorim, swaero_diag, nfiles
     456    USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    453457    IMPLICIT NONE
    454458
    455     INCLUDE "dimensions.h"
    456     INCLUDE "temps.h"
    457459    INCLUDE "clesphys.h"
    458460
    … …  
    478480       IF ( flag_var(iff)<=lev_files(iff) ) THEN
    479481          CALL histdef (nid_files(iff),nomvar,titrevar,unitvar, &
    480                iim,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
     482               nbp_lon,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
    481483               type_ecri(iff), zstophym,zoutm(iff))               
    482484       ENDIF
    … …  
    507509                                   nhorim, zdtime_moy, levmin, levmax, &
    508510                                   nvertm, nfiles
     511    USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    509512    IMPLICIT NONE
    510513
    511     INCLUDE "dimensions.h"
    512     INCLUDE "temps.h"
    513 !    INCLUDE "indicesol.h"
    514514    INCLUDE "clesphys.h"
    515515
    … …  
    535535       IF ( flag_var(iff)<=lev_files(iff) ) THEN
    536536          CALL histdef (nid_files(iff), nomvar, titrevar, unitvar, &
    537                iim, jj_nb, nhorim(iff), klev, levmin(iff), &
     537               nbp_lon, jj_nb, nhorim(iff), klev, levmin(iff), &
    538538               levmax(iff)-levmin(iff)+1, nvertm(iff), 32, type_ecri(iff), &
    539539               zstophym, zoutm(iff))
    … …  
    564564                                   nid_files, nhorim, swaero_diag
    565565    USE print_control_mod, ONLY: prt_level,lunout
     566    USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    566567#ifdef CPP_XIOS
    567568    use wxios, only: wxios_add_field_to_file
    … …  
    569570    IMPLICIT NONE
    570571
    571     INCLUDE "dimensions.h"
    572     INCLUDE "temps.h"
    573572    INCLUDE "clesphys.h"
    574573
    … …  
    621620       IF ( var%flag(iff)<=lev_files(iff) ) THEN
    622621          CALL histdef (nid_files(iff), var%name, var%description, var%unit, &
    623                iim,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
     622               nbp_lon,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
    624623               typeecrit, zstophym,zoutm(iff))               
    625624       ENDIF
    … …  
    651650                                   levmax, nvertm
    652651    USE print_control_mod, ONLY: prt_level,lunout
     652    USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    653653#ifdef CPP_XIOS
    654654    use wxios, only: wxios_add_field_to_file
    … …  
    656656    IMPLICIT NONE
    657657
    658     INCLUDE "dimensions.h"
    659     INCLUDE "temps.h"
    660658    INCLUDE "clesphys.h"
    661659
    … …  
    708706       IF ( var%flag(iff)<=lev_files(iff) ) THEN
    709707          CALL histdef (nid_files(iff), var%name, var%description, var%unit, &
    710                iim, jj_nb, nhorim(iff), klev, levmin(iff), &
     708               nbp_lon, jj_nb, nhorim(iff), klev, levmin(iff), &
    711709               levmax(iff)-levmin(iff)+1, nvertm(iff), 32, typeecrit, &
    712710               zstophym, zoutm(iff))
    … …  
    750748  USE ioipsl, only: histwrite
    751749  USE print_control_mod, ONLY: prt_level,lunout
     750  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    752751  IMPLICIT NONE
    753   include 'dimensions.h'
    754752   
    755753    integer,INTENT(IN) :: nid
    … …  
    760758    REAL,DIMENSION(klon_mpi) :: buffer_omp
    761759    INTEGER, allocatable, DIMENSION(:) :: index2d
    762     REAL :: Field2d(iim,jj_nb)
     760    REAL :: Field2d(nbp_lon,jj_nb)
    763761
    764762    integer :: ip
    … …  
    772770    CALL grid1Dto2D_mpi(buffer_omp,Field2d)
    773771    if(.NOT.lpoint) THEN
    774      ALLOCATE(index2d(iim*jj_nb))
    775      ALLOCATE(fieldok(iim*jj_nb))
     772     ALLOCATE(index2d(nbp_lon*jj_nb))
     773     ALLOCATE(fieldok(nbp_lon*jj_nb))
    776774     IF (prt_level >= 10) write(lunout,*)'Sending ',name,' to IOIPSL'
    777      CALL histwrite(nid,name,itau,Field2d,iim*jj_nb,index2d)
     775     CALL histwrite(nid,name,itau,Field2d,nbp_lon*jj_nb,index2d)
    778776     IF (prt_level >= 10) write(lunout,*)'Finished sending ',name,' to IOIPSL'
    779777    else
    … …  
    813811                                is_sequential, klon_mpi_begin, klon_mpi_end, &
    814812                                jj_nb, klon_mpi
     813  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    815814  use ioipsl, only: histwrite
    816815  USE print_control_mod, ONLY: prt_level,lunout
    817816  IMPLICIT NONE
    818   include 'dimensions.h'
    819817   
    820818    integer,INTENT(IN) :: nid
    … …  
    824822    REAL,DIMENSION(:,:),INTENT(IN) :: field  ! --> field(klon,:)
    825823    REAL,DIMENSION(klon_mpi,size(field,2)) :: buffer_omp
    826     REAL :: Field3d(iim,jj_nb,size(field,2))
     824    REAL :: Field3d(nbp_lon,jj_nb,size(field,2))
    827825    INTEGER :: ip, n, nlev
    828826    INTEGER, ALLOCATABLE, DIMENSION(:) :: index3d
    … …  
    837835    CALL grid1Dto2D_mpi(buffer_omp,field3d)
    838836    if(.NOT.lpoint) THEN
    839      ALLOCATE(index3d(iim*jj_nb*nlev))
    840      ALLOCATE(fieldok(iim*jj_nb,nlev))
     837     ALLOCATE(index3d(nbp_lon*jj_nb*nlev))
     838     ALLOCATE(fieldok(nbp_lon*jj_nb,nlev))
    841839     IF (prt_level >= 10) write(lunout,*)'Sending ',name,' to IOIPSL'
    842      CALL histwrite(nid,name,itau,Field3d,iim*jj_nb*nlev,index3d)
     840     CALL histwrite(nid,name,itau,Field3d,nbp_lon*jj_nb*nlev,index3d)
    843841     IF (prt_level >= 10) write(lunout,*)'Finished sending ',name,' to IOIPSL'
    844842   else
    … …  
    889887                                 nid_files
    890888  USE print_control_mod, ONLY: prt_level,lunout
     889  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    891890#ifdef CPP_XIOS
    892891  USE xios, only: xios_send_field
    … …  
    895894
    896895  IMPLICIT NONE
    897   INCLUDE 'dimensions.h'
    898896  include 'clesphys.h'
    899897
    … …  
    908906    REAL,DIMENSION(klon_mpi) :: buffer_omp
    909907    INTEGER, allocatable, DIMENSION(:) :: index2d
    910     REAL :: Field2d(iim,jj_nb)
     908    REAL :: Field2d(nbp_lon,jj_nb)
    911909
    912910    INTEGER :: ip
    … …  
    989987
    990988                  IF(.NOT.clef_stations(iff)) THEN
    991                         ALLOCATE(index2d(iim*jj_nb))
    992                         ALLOCATE(fieldok(iim*jj_nb))
     989                        ALLOCATE(index2d(nbp_lon*jj_nb))
     990                        ALLOCATE(fieldok(nbp_lon*jj_nb))
    993991#ifndef CPP_IOIPSL_NO_OUTPUT
    994                         CALL histwrite(nid_files(iff),var%name,itau_iophy,Field2d,iim*jj_nb,index2d)
     992                        CALL histwrite(nid_files(iff),var%name,itau_iophy,Field2d,nbp_lon*jj_nb,index2d)
    995993#endif
    996994!#ifdef CPP_XIOS
    … …  
    10481046                                 nfiles, vars_defined, clef_stations, &
    10491047                                 nid_files
     1048  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    10501049#ifdef CPP_XIOS
    10511050  USE xios, only: xios_send_field
    … …  
    10541053
    10551054  IMPLICIT NONE
    1056   INCLUDE 'dimensions.h'
    10571055  include 'clesphys.h'
    10581056
    … …  
    10651063!$OMP THREADPRIVATE(firstx)
    10661064    REAL,DIMENSION(klon_mpi,SIZE(field,2)) :: buffer_omp
    1067     REAL :: Field3d(iim,jj_nb,SIZE(field,2))
     1065    REAL :: Field3d(nbp_lon,jj_nb,SIZE(field,2))
    10681066    INTEGER :: ip, n, nlev, nlevx
    10691067    INTEGER, ALLOCATABLE, DIMENSION(:) :: index3d
    … …  
    11311129                  write(lunout,*)'histwrite3d_phy: call xios_send_field for ', &
    11321130                                  trim(var%name), ' with iim jjm nlevx = ', &
    1133                                   iim,jj_nb,nlevx
     1131                                  nbp_lon,jj_nb,nlevx
    11341132                endif
    11351133                CALL xios_send_field(var%name, Field3d(:,:,1:nlevx))
    … …  
    11381136#endif
    11391137                IF (.NOT.clef_stations(iff)) THEN
    1140                         ALLOCATE(index3d(iim*jj_nb*nlev))
    1141                         ALLOCATE(fieldok(iim*jj_nb,nlev))
     1138                        ALLOCATE(index3d(nbp_lon*jj_nb*nlev))
     1139                        ALLOCATE(fieldok(nbp_lon*jj_nb,nlev))
    11421140
    11431141#ifndef CPP_IOIPSL_NO_OUTPUT
    1144                         CALL histwrite(nid_files(iff),var%name,itau_iophy,Field3d,iim*jj_nb*nlev,index3d)
     1142                        CALL histwrite(nid_files(iff),var%name,itau_iophy,Field3d,nbp_lon*jj_nb*nlev,index3d)
    11451143#endif
    11461144
    … …  
    11941192                                is_sequential, klon_mpi_begin, klon_mpi_end, &
    11951193                                jj_nb, klon_mpi
     1194  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    11961195  USE xios, only: xios_send_field
    11971196  USE print_control_mod, ONLY: prt_level,lunout
    11981197
    11991198  IMPLICIT NONE
    1200   INCLUDE 'dimensions.h'
    12011199
    12021200    CHARACTER(LEN=*), INTENT(IN) :: field_name
    … …  
    12051203    REAL,DIMENSION(klon_mpi) :: buffer_omp
    12061204    INTEGER, allocatable, DIMENSION(:) :: index2d
    1207     REAL :: Field2d(iim,jj_nb)
     1205    REAL :: Field2d(nbp_lon,jj_nb)
    12081206
    12091207    INTEGER :: ip
    … …  
    12241222    !IF(.NOT.clef_stations(iff)) THEN
    12251223    IF (.TRUE.) THEN
    1226         ALLOCATE(index2d(iim*jj_nb))
    1227         ALLOCATE(fieldok(iim*jj_nb))
     1224        ALLOCATE(index2d(nbp_lon*jj_nb))
     1225        ALLOCATE(fieldok(nbp_lon*jj_nb))
    12281226
    12291227
    … …  
    12651263                                jj_nb, klon_mpi
    12661264  USE xios, only: xios_send_field
     1265  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
    12671266  USE print_control_mod, ONLY: prt_level,lunout
    12681267
    12691268  IMPLICIT NONE
    1270   INCLUDE 'dimensions.h'
    12711269
    12721270    CHARACTER(LEN=*), INTENT(IN) :: field_name
    … …  
    12741272
    12751273    REAL,DIMENSION(klon_mpi,SIZE(field,2)) :: buffer_omp
    1276     REAL :: Field3d(iim,jj_nb,SIZE(field,2))
     1274    REAL :: Field3d(nbp_lon,jj_nb,SIZE(field,2))
    12771275    INTEGER :: ip, n, nlev
    12781276    INTEGER, ALLOCATABLE, DIMENSION(:) :: index3d
    … …  
    12951293    !IF (.NOT.clef_stations(iff)) THEN
    12961294    IF(.TRUE.)THEN
    1297         ALLOCATE(index3d(iim*jj_nb*nlev))
    1298         ALLOCATE(fieldok(iim*jj_nb,nlev))
     1295        ALLOCATE(index3d(nbp_lon*jj_nb*nlev))
     1296        ALLOCATE(fieldok(nbp_lon*jj_nb,nlev))
    12991297        CALL xios_send_field(field_name, Field3d(:,:,1:nlev))
    13001298                       
Note: See TracChangeset for help on using the changeset viewer.