/[MITgcm]/MITgcm/pkg/ptracers/ptracers_mnc_init.F
ViewVC logotype

Contents of /MITgcm/pkg/ptracers/ptracers_mnc_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.2 - (show annotations) (download)
Sat Sep 4 03:08:03 2004 UTC (19 years, 9 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint55b_post, checkpoint55, checkpoint54f_post, checkpoint55a_post
Changes since 1.1: +10 -3 lines
 o add ability to output ptracers timeave fields to NetCDF
   - small cleanups to verification/dic_example

1 C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_mnc_init.F,v 1.1 2004/09/03 20:10:47 edhill Exp $
2 C $Name: $
3
4 #include "PTRACERS_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP
8 C !ROUTINE:
9
10 C !INTERFACE:
11 SUBROUTINE PTRACERS_MNC_INIT( myThid )
12
13 C !DESCRIPTION:
14 C Initialize PTRACERS parameters, read in data.ptracers
15
16 C !USES:
17 IMPLICIT NONE
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PTRACERS_SIZE.h"
21 #include "PTRACERS.h"
22 #include "PARAMS.h"
23
24 C !INPUT PARAMETERS:
25 INTEGER myThid
26 CEOP
27
28 #ifdef ALLOW_PTRACERS
29 #ifdef ALLOW_MNC
30
31 C !LOCAL VARIABLES:
32 INTEGER ii, jj, ilnb
33 character*(80) name
34 character*(MAX_LEN_MBUF) msgbuf
35
36 C Functions
37 integer ILNBLNK
38
39 DO ii = 1,PTRACERS_numInUse
40
41 DO jj = 1,80
42 name(jj:jj) = ' '
43 ENDDO
44
45 C Create the tracer
46 ilnb = ILNBLNK(PTRACERS_names(ii))
47 IF (ilnb .LT. 1) THEN
48 write(msgbuf,'(2a,i3,2a)') 'PTRACERS_MNC_INIT ERROR: ',
49 & 'PTRACERS_names(',ii,') is blank--please name it in',
50 & ' ''data.ptracers'''
51 CALL print_error(msgbuf, mythid)
52 stop 'ABNORMAL END: S/R PTRACERS_MNC_INIT'
53 ENDIF
54 CALL MNC_CW_ADD_VNAME(
55 & PTRACERS_names(ii),'Cen_xy_Hn__C__t',4,5,myThid)
56 WRITE(name,'(a,a)') 'surf_',PTRACERS_names(ii)(1:ilnb)
57 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__-__t',3,4,myThid)
58
59 C Add the long_name attribute if its defined
60 ilnb = ILNBLNK(PTRACERS_long_names(ii))
61 IF (ilnb .GT. 0) THEN
62 CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),1,'long_name',
63 & PTRACERS_long_names(ii),myThid)
64 ENDIF
65
66 C Add the units attribute if its defined
67 ilnb = ILNBLNK(PTRACERS_units(ii))
68 IF (ilnb .GT. 0) THEN
69 CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),1,'units',
70 & PTRACERS_units(ii),myThid)
71 ENDIF
72
73 ENDDO
74
75 #endif /* ALLOW_MNC */
76 #endif /* ALLOW_PTRACERS */
77
78 RETURN
79 END
80

  ViewVC Help
Powered by ViewVC 1.1.22