pkg/ptracers/ptracers_mnc_init.F

C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_mnc_init.F,v 1.2 2004/09/04 03:08:03 edhill Exp $
C $Name:  $

#include "PTRACERS_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP
C     !ROUTINE: 
      
C     !INTERFACE:
      SUBROUTINE PTRACERS_MNC_INIT( myThid )
      
C     !DESCRIPTION:
C     Initialize PTRACERS parameters, read in data.ptracers

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PTRACERS_SIZE.h"
#include "PTRACERS.h"
#include "PARAMS.h"
#ifdef ALLOW_MNC
#include "MNC_PARAMS.h"
#endif

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

#ifdef ALLOW_PTRACERS
#ifdef ALLOW_MNC

C     !LOCAL VARIABLES:
      INTEGER ii, jj, ilnb
      character*(80) name
      character*(MAX_LEN_MBUF) msgbuf

C     Functions
      integer ILNBLNK

      IF (useMNC) THEN
        DO ii = 1,PTRACERS_numInUse

          DO jj = 1,80
            name(jj:jj) = ' '
          ENDDO

C         Create the tracer
          ilnb = ILNBLNK(PTRACERS_names(ii))
          IF (ilnb .LT. 1) THEN
            write(msgbuf,'(2a,i3,2a)') 'PTRACERS_MNC_INIT ERROR: ',
     &           'PTRACERS_names(',ii,') is blank--please name it in',
     &           ' ''data.ptracers'''
            CALL print_error(msgbuf, mythid)
            stop 'ABNORMAL END: S/R PTRACERS_MNC_INIT'
          ENDIF
          CALL MNC_CW_ADD_VNAME(
     &         PTRACERS_names(ii),'Cen_xy_Hn__C__t',4,5,myThid)
          WRITE(name,'(a,a)') 'surf_',PTRACERS_names(ii)(1:ilnb)
          CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__-__t',3,4,myThid)

C         Add the long_name attribute if its defined
          ilnb = ILNBLNK(PTRACERS_long_names(ii))
          IF (ilnb .GT. 0) THEN
            CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),
     &           'long_name',PTRACERS_long_names(ii),myThid)
          ENDIF

C         Add the units attribute if its defined
          ilnb = ILNBLNK(PTRACERS_units(ii))
          IF (ilnb .GT. 0) THEN
            CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),'units',
     &           PTRACERS_units(ii),myThid)
          ENDIF
          
        ENDDO
        
      ENDIF
#endif /*  ALLOW_MNC  */
#endif /* ALLOW_PTRACERS */

      RETURN
      END

1	C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_mnc_init.F,v 1.2 2004/09/04 03:08:03 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	#ifdef ALLOW_MNC
24	#include "MNC_PARAMS.h"
25	#endif
26
27	C !INPUT PARAMETERS:
28	INTEGER myThid
29	CEOP
30
31	#ifdef ALLOW_PTRACERS
32	#ifdef ALLOW_MNC
33
34	C !LOCAL VARIABLES:
35	INTEGER ii, jj, ilnb
36	character*(80) name
37	character*(MAX_LEN_MBUF) msgbuf
38
39	C Functions
40	integer ILNBLNK
41
42	IF (useMNC) THEN
43	DO ii = 1,PTRACERS_numInUse
44
45	DO jj = 1,80
46	name(jj:jj) = ' '
47	ENDDO
48
49	C Create the tracer
50	ilnb = ILNBLNK(PTRACERS_names(ii))
51	IF (ilnb .LT. 1) THEN
52	write(msgbuf,'(2a,i3,2a)') 'PTRACERS_MNC_INIT ERROR: ',
53	& 'PTRACERS_names(',ii,') is blank--please name it in',
54	& ' ''data.ptracers'''
55	CALL print_error(msgbuf, mythid)
56	stop 'ABNORMAL END: S/R PTRACERS_MNC_INIT'
57	ENDIF
58	CALL MNC_CW_ADD_VNAME(
59	& PTRACERS_names(ii),'Cen_xy_Hn__C__t',4,5,myThid)
60	WRITE(name,'(a,a)') 'surf_',PTRACERS_names(ii)(1:ilnb)
61	CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__-__t',3,4,myThid)
62
63	C Add the long_name attribute if its defined
64	ilnb = ILNBLNK(PTRACERS_long_names(ii))
65	IF (ilnb .GT. 0) THEN
66	CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),
67	& 'long_name',PTRACERS_long_names(ii),myThid)
68	ENDIF
69
70	C Add the units attribute if its defined
71	ilnb = ILNBLNK(PTRACERS_units(ii))
72	IF (ilnb .GT. 0) THEN
73	CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),'units',
74	& PTRACERS_units(ii),myThid)
75	ENDIF
76
77	ENDDO
78
79	ENDIF
80	#endif /* ALLOW_MNC */
81	#endif /* ALLOW_PTRACERS */
82
83	RETURN
84	END
85