!-------------------------------------- LICENCE BEGIN ------------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
! version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
!without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
!See the above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
!
subroutine sumvogp 1,1
c
#if defined (DOC)
*
***s/r sumvogp - Get single profiles of background fields and store in
* in GOMOBSG for ONEOBS and effective variances experiments
*
*Author : C. Charette *ARMA/AES - 20 nov 98
*Revision:
*
* -------------------
** Purpose: Setup for Project Background statistics into observation space
* . using single profiles for background fields
*Arguments
*
#endif
IMPLICIT NONE
*implicits
#include "comdim.cdk"
#include "comdimo.cdk"
#include "comlun.cdk"
#include "commvog.cdk"
*
C
INTEGER IERR,ILEN,IDATA,JDATA,JK
REAL*8 PTOP(NOBTOT)
POINTER(PPTOP,PTOP)
C
C* 0.1 Memory allocation
C . -----------------
ILEN = NOBTOT
CALL HPALLOC(PPTOP,MAX(ILEN,1),IERR,8)
C
C* 0.2 Read background profiles
C
write(nulout,*)'SUMVOGP: ********UNDER CONSTRUCTION*****'
CALL ABORT3D(NULOUT,'SUMVOGP: UNDER CONSTRUCTION ')
c ***********************************************************
write(NULOUT,*) 'SETTING GOMOBSG TO FIXED PROFILES'
c
c FIXED PROFILES
c
do jdata=1,NOBTOT
c surface pressure
gompsg(1,jdata)=100231.4208984378
c log specific humidity
do jk=1,NFLEV
gomug(jk,jdata)=0.0
gomvg(jk,jdata)=0.0
enddo
gomqg(1 ,IDATA)= -12.89924640895085
gomqg(2 ,IDATA)= -12.89924640895085
gomqg(3 ,IDATA)= -12.89937593282709
gomqg(4 ,IDATA)= -12.90796218272754
gomqg(5 ,IDATA)= -12.58317461534159
gomqg(6 ,IDATA)= -11.50597894209372
gomqg(7 ,IDATA)= -11.38417325284509
gomqg(8 ,IDATA)= -11.08485494285742
gomqg(9 ,IDATA)= -10.93366738616775
gomqg(10 ,IDATA)= -11.00480369922830
gomqg(11 ,IDATA)= -10.91101578224992
gomqg(12 ,IDATA)= -10.00086176166695
gomqg(13 ,IDATA)= -9.061022603352350
gomqg(14 ,IDATA)= -8.347015549520336
gomqg(15 ,IDATA)= -7.641485977618439
gomqg(16 ,IDATA)= -7.103567127246686
gomqg(17 ,IDATA)= -6.754322343325619
gomqg(18 ,IDATA)= -6.619065452093341
gomqg(19 ,IDATA)= -6.699411220287870
gomqg(20 ,IDATA)= -6.603065971906735
gomqg(21 ,IDATA)= -6.305291476708035
gomqg(22 ,IDATA)= -6.097815281374210
gomqg(23 ,IDATA)= -5.874184595215093
gomqg(24 ,IDATA)= -5.806009184826604
gomqg(25 ,IDATA)= -5.774973255356879
gomqg(26 ,IDATA)= -5.769083987832715
gomqg(27 ,IDATA)= -5.733790495054194
gomqg(28 ,IDATA)= -5.579646889562327
c temperature
gomtg(1 ,IDATA)= 226.2727111816404
gomtg(2 ,IDATA)= 222.2787689208982
gomtg(3 ,IDATA)= 217.5619873046873
gomtg(4 ,IDATA)= 213.3916381835936
gomtg(5 ,IDATA)= 213.7779144287107
gomtg(6 ,IDATA)= 212.2317413330076
gomtg(7 ,IDATA)= 212.4781097412107
gomtg(8 ,IDATA)= 213.3258728027342
gomtg(9 ,IDATA)= 211.5532135009764
gomtg(10 ,IDATA)= 209.3011840820311
gomtg(11 ,IDATA)= 211.8721221923826
gomtg(12 ,IDATA)= 219.1873229980467
gomtg(13 ,IDATA)= 228.0132354736326
gomtg(14 ,IDATA)= 236.0305541992186
gomtg(15 ,IDATA)= 243.7344421386717
gomtg(16 ,IDATA)= 249.5092529296873
gomtg(17 ,IDATA)= 253.8461669921873
gomtg(18 ,IDATA)= 256.0853637695311
gomtg(19 ,IDATA)= 256.6787017822264
gomtg(20 ,IDATA)= 259.5489410400389
gomtg(21 ,IDATA)= 263.2389739990232
gomtg(22 ,IDATA)= 266.5882457733154
gomtg(23 ,IDATA)= 269.0130466461181
gomtg(24 ,IDATA)= 271.1701196670531
gomtg(25 ,IDATA)= 273.6114745974540
gomtg(26 ,IDATA)= 275.7961935043335
gomtg(27 ,IDATA)= 276.8473419189453
gomtg(28 ,IDATA)= 277.6344732284546
c gomtg(28,IDATA)= 288.2
c gomtg(27,IDATA)= 287.5
c gomtg(26,IDATA)= 286.0
c gomtg(25,IDATA)= 284.9
c gomtg(24,IDATA)= 283.5
c gomtg(23,IDATA)= 281.7
c gomtg(22,IDATA)= 278.4
c gomtg(21,IDATA)= 275.0
c gomtg(20,IDATA)= 271.0
c gomtg(19,IDATA)= 268.7
c gomtg(18,IDATA)= 265.4
c gomtg(17,IDATA)= 258.9
c gomtg(16,IDATA)= 255.7
c gomtg(15,IDATA)= 249.2
c gomtg(14,IDATA)= 242.7
c gomtg(13,IDATA)= 236.2
c gomtg(12,IDATA)= 229.7
c gomtg(11,IDATA)= 223.3
c gomtg(10,IDATA)= 216.8
c gomtg(9,IDATA)= 216.7
c gomtg(8,IDATA)= 216.7
c gomtg(7,IDATA)= 216.7
c gomtg(6,IDATA)= 216.7
c gomtg(5,IDATA)= 216.7
c gomtg(4,IDATA)= 216.7
c gomtg(3,IDATA)= 218.6
c gomtg(2,IDATA)= 222.5
c gomtg(1,IDATA)= 227.5
enddo
c
return
end