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

Changeset 4108


Ignore:
Timestamp:
2013-10-22T22:27:13+02:00 (10 years ago)
Author:
andrewryan
Message:

Fix reference to nemo_gcm in nemogcm.F90 in OOO_SRC. Updated write statements in ooo_read ooo_write modules.

Location:
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/nemogcm.F90

    r4106 r4108  
    8888   PRIVATE 
    8989 
    90    PUBLIC   nemo_gcm    ! called by model.F90 
    9190   PUBLIC   nemo_init   ! needed by AGRIF 
    9291   PUBLIC   nemo_alloc  ! needed by TAM 

  • branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/ooo_read.F90

    r4106 r4108  
    181181         sshn(:,:) = fbrmdi  
    182182      ELSE 
    183          WRITE(*,*) "Opening :", TRIM(filename) 
     183         WRITE(numout,*) "Opening :", TRIM(filename) 
    184184         ! Open Netcdf file to find dimension id 
    185185         istat = nf90_open(path=TRIM(filename), mode=nf90_nowrite, ncid=ncid) 
    186186         IF ( istat /= nf90_noerr ) THEN 
    187              WRITE(*,*) "WARNING: Could not open ", trim(filename) 
    188              WRITE(*,*) "ERROR: ", nf90_strerror(istat) 
     187             WRITE(numout,*) "WARNING: Could not open ", trim(filename) 
     188             WRITE(numout,*) "ERROR: ", nf90_strerror(istat) 
    189189         ENDIF 
    190190         istat = nf90_inq_dimid(ncid,'x',xdim) 
    … …  
    265265         ENDIF 
    266266         ! Close netcdf file 
    267          WRITE(*,*) "Closing :", TRIM(filename) 
     267         WRITE(numout,*) "Closing :", TRIM(filename) 
    268268         istat = nf90_close(ncid) 
    269269      END IF 

  • branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/ooo_write.F90

    r4106 r4108  
    33   !!                       ***  MODULE ooo_write  *** 
    44   !!====================================================================== 
     5   USE in_out_manager 
    56   USE netcdf 
    67   USE obs_utils, ONLY: chkerr 
    … …  
    194195         ! longitude 
    195196         cvar = 'longitude' 
    196          WRITE(*,*) TRIM(cvar) 
    197197         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    198198         CALL chkerr(nf90_put_var(ncid, varid, plam),cpname, __LINE__ ) 
    199199         ! latitude 
    200200         cvar = 'latitude' 
    201          WRITE(*,*) TRIM(cvar) 
    202201         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    203202         CALL chkerr(nf90_put_var(ncid, varid, pphi),cpname, __LINE__ ) 
    204203         ! depth 
    205204         cvar = 'depth' 
    206          WRITE(*,*) TRIM(cvar) 
    207205         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    208206         CALL chkerr(nf90_put_var(ncid, varid, pdep),cpname, __LINE__ ) 
    209207         ! varname 
    210208         cvar = 'varname' 
    211          WRITE(*,*) TRIM(cvar) 
    212209         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    213210         CALL chkerr(nf90_put_var(ncid, varid, cvnam),cpname, __LINE__ ) 
    214211         ! unitname 
    215212         cvar = 'unitname' 
    216          WRITE(*,*) TRIM(cvar) 
    217213         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    218214         CALL chkerr(nf90_put_var(ncid, varid, cunit),cpname, __LINE__ ) 
    219215         ! leadtime 
    220216         cvar = 'leadtime' 
    221          WRITE(*,*) TRIM(cvar) 
    222217         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    223218         CALL chkerr(nf90_put_var(ncid, varid, plead),cpname, __LINE__ ) 
    224219         ! observation 
    225220         cvar = 'observation' 
    226          WRITE(*,*) TRIM(cvar) 
    227221         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    228222         CALL chkerr(nf90_put_var(ncid, varid, pob),cpname, __LINE__ ) 
    229223         ! qc 
    230224         cvar = 'qc' 
    231          WRITE(*,*) TRIM(cvar) 
    232225         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    233226         CALL chkerr(nf90_put_var(ncid, varid, kqc),cpname, __LINE__ ) 
    234227         ! juld 
    235228         cvar = 'juld' 
    236          WRITE(*,*) TRIM(cvar) 
    237229         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    238230         CALL chkerr(nf90_put_var(ncid, varid, ptim),cpname, __LINE__ ) 
    239231         ! type 
    240232         cvar = 'type' 
    241          WRITE(*,*) TRIM(cvar) 
    242233         CALL inst_converter(cdtyp,nobs,cdpad) 
    243234         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    … …  
    245236         ! id 
    246237         cvar = 'id' 
    247          WRITE(*,*) TRIM(cvar) 
    248238         CALL chkerr(nf90_inq_varid(ncid,TRIM(cvar), varid),cpname, __LINE__ ) 
    249239         CALL chkerr(nf90_put_var(ncid, varid, cdwmo),cpname, __LINE__ ) 
    … …  
    346336         istat = nf90_open(TRIM(cdfilename), nf90_nowrite, ncid) 
    347337         IF (istat /= nf90_noerr) THEN 
    348             WRITE(*, *) TRIM(cdfilename) 
     338            IF (lwp) WRITE(numout, *) TRIM(cdfilename), ' opened successfully.' 
    349339            CALL chkerr( nf90_create(TRIM(cdfilename), nf90_clobber, ncid), cpname, __LINE__ ) 
    350340            !! Global Dimension section 
    … …  
    395385            !  Forecast day 
    396386            cvar = 'leadtime' 
    397             WRITE(*,*) TRIM(cvar) 
    398387            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dimid_f, varid), cpname, __LINE__ ) 
    399388            cattnam = 'long_name' 
    … …  
    409398            !  longitude 
    410399            cvar = 'longitude' 
    411             WRITE(*,*) TRIM(cvar) 
    412400            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dimid_o, varid), cpname, __LINE__ ) 
    413401            cattnam = 'long_name' 
    … …  
    420408            !  latitude 
    421409            cvar = 'latitude' 
    422             WRITE(*,*) TRIM(cvar) 
    423410            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dimid_o, varid), cpname, __LINE__ ) 
    424411            cattnam = 'long_name' 
    … …  
    430417            !  depth 
    431418            cvar = 'depth' 
    432             WRITE(*,*) TRIM(cvar) 
    433419            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim2a, varid), cpname, __LINE__ ) 
    434420            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    442428            !  varname 
    443429            cvar = 'varname' 
    444             WRITE(*,*) TRIM(cvar) 
    445430            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_char, dim2c, varid), cpname, __LINE__ ) 
    446431            cattnam = 'long_name' 
    … …  
    450435            !  unitname 
    451436            cvar = 'unitname' 
    452             WRITE(*,*) TRIM(cvar) 
    453437            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_char, dim2c, varid), cpname, __LINE__ ) 
    454438            cattnam = 'long_name' 
    … …  
    458442            !  observation 
    459443            cvar = 'observation' 
    460             WRITE(*,*) TRIM(cvar) 
    461444            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim3a, varid), cpname, __LINE__ ) 
    462445            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    467450            !  forecast 
    468451            cvar = 'forecast' 
    469             WRITE(*,*) TRIM(cvar) 
    470452            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim4a, varid), cpname, __LINE__ ) 
    471453            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    479461            !  persistence 
    480462            cvar = 'persistence' 
    481             WRITE(*,*) TRIM(cvar) 
    482463            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim4a, varid), cpname, __LINE__ ) 
    483464            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    491472            !  climatology 
    492473            cvar = 'climatology' 
    493             WRITE(*,*) TRIM(cvar) 
    494474            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim3a, varid), cpname, __LINE__ ) 
    495475            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    503483            !  daym2 
    504484            cvar = 'best_estimate' 
    505             WRITE(*,*) TRIM(cvar) 
    506485            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim3a, varid), cpname, __LINE__ ) 
    507486            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    515494            !  daym1 
    516495            cvar = 'nrt_analysis' 
    517             WRITE(*,*) TRIM(cvar) 
    518496            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim3a, varid), cpname, __LINE__ ) 
    519497            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    529507               !  mdt 
    530508               cvar = 'mdt_reference' 
    531                WRITE(*,*) TRIM(cvar) 
    532509               CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim3a, varid), cpname, __LINE__ ) 
    533510               CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    541518               !  altbias 
    542519               cvar = 'altimeter_bias' 
    543                WRITE(*,*) TRIM(cvar) 
    544520               CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dim3a, varid), cpname, __LINE__ ) 
    545521               CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    551527            !  qc 
    552528            cvar = 'qc' 
    553             WRITE(*,*) TRIM(cvar) 
    554529            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_short, dim3a, varid), cpname, __LINE__ ) 
    555530            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', NF90_FILL_SHORT), cpname, __LINE__ ) 
    … …  
    568543!            !  juld 
    569544            cvar = 'juld' 
    570             WRITE(*,*) TRIM(cvar) 
    571545            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_double, dimid_o, varid), cpname, __LINE__ ) 
    572546!            CALL chkerr(nf90_put_att(ncid, varid, '_FillValue', clrmdi), cpname, __LINE__ ) 
    … …  
    581555            !  modeljuld  
    582556            cvar = 'modeljuld' 
    583             WRITE(*,*) TRIM(cvar) 
    584557            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_float, dimid_f, varid), cpname, __LINE__ ) 
    585558            cattnam = 'long_name' 
    … …  
    592565            !  type 
    593566            cvar = 'type' 
    594             WRITE(*,*) TRIM(cvar) 
    595567            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_char, dim2d, varid), cpname, __LINE__ ) 
    596568            cattnam = 'long_name' 
    … …  
    600572            !  id 
    601573            cvar = 'id' 
    602             WRITE(*,*) TRIM(cvar) 
    603574            CALL chkerr(nf90_def_var(ncid, TRIM(cvar), nf90_char, dim2b, varid), cpname, __LINE__ ) 
    604575            cattnam = 'long_name' 
    … …  
    609580            CALL chkerr(nf90_close(ncid), cpname, __LINE__ ) 
    610581         ELSE 
    611             WRITE(*, *) cdfilename, 'already exists.' 
     582            IF (lwp) WRITE(numout, *) TRIM(cdfilename), 'already exists.' 
    612583            ! Close netcdf file 
    613584            CALL chkerr(nf90_close(ncid), cpname, __LINE__ ) 
Note: See TracChangeset for help on using the changeset viewer.